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1 6 Abstract 

A h i gh- t i p-speed, low-loading transonic fan stage was designed to deliver an overall pressure 
ratio of 1.5 with an adiabatic efficiency of 86 percent. The design flow per unit annulus area 
is 42.0 lbm/s-ft 2 (205.1 kgm/s-m 2 ). The fan features a hub/tip ratio of 0.46, a tip diameter 
of 28.74 in. (0.73 m) and operates at a design tip speed of 1600 fps (488.6 m/s). For these 
design conditions, -the rotor blade tip region operates with supersonic inlet and supersonic 
discharge relative velocities. A sophisticated quas i- three-d imens ional characteristic section 
design procedure was used for the all-supersonic sections and the inlet of the midspan transonic 
sections. For regions where the relative outlet velocities are supersonic, the blade operates 
with weak oblique shocks only. 
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SUMMARY 


A high-tip-speed, 1 ow-pressure-rat i o transonic fan stage was designed. 

The objective of this design was the provision of a rotor design that could 
deliver good efficiency at low work input via elimination of strong shock 
losses and shock- i nduced separation in the h igh-M^ch-number tip region. In 
the rotor blade tip region, the design incorporates near oblique shocks only, 
with supersonic relative outlet velocities. In order to accomplish thi$, a 
relatively sophisticated quas i-three-dimens ional characteristic section design 
procedure was used for the a 1 1 -superson i c sections ancj for the inlet of the 
midspan transonic sections. 

The subsonic hub sections were designed on the basis of conventional 
empirical cascade data and sections of the three regions (tip, pitch, and hub) 
were analytically joined via a special stage-stacking program to form an aero- 
dynamical ly and structurally compatible blade. Special attention was focused 
on inlet losses, the boundary layer growth rate, and the physical passage con- 
traction ratio to ensure passage "starting 11 below the rotor design speed. 

Although the rotor relative velocities were well into the supersonic range, 
the stator velocities were at modest, subsonic levels. Hence, the stator design 
is conventionally subsonic and has the primary role of turning the flow back to 
the axial d i rect i on . 


The final fan design calculations 
stage design parameters. 

Overall pressure ratio 

Adiabatic efficiency 

Equivalent total flow 

Flow per unit annulus area 

Equivalent tip speed 

Inlet hub-to-tip radius ratio 

Tip diameter 


resulted in the following values of 

1.5 

0.86 

148 lb/s (6.71 kgm/s) 

42.0 lb/s- ft 2 (205.1 kgm/s-m 2 ) 

1600 ft/s (488 . 6 m/s) 

0.46 

28.74 in. (0.73 m) 


It is anticipated that stages designed as described herein would permit 
the selection of h i gh-bypass- rat io fan speeds which minimize the number of fan 
drive turbine stages wi thout compromising fan or fan turbine efficiency. Other 
supplementary benefits, such as reduced hub loading and/or the ability to uti- 
lize lower radius ratios for a given stage pressure ratio, should also result. 
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INTRODUCTION 


The recent need for power plants to propel large transport aircraft capable 
of long-range flights has spurred the development of very-h i gh-bypass-rat i o 
turbofan engines (bypass ratios between 8 and 12 or even higher are being 
explored). The fan component of these high-bypass-ratio engines, in general, 
controls the engine diameter and requires large drive-turbine work. In addi- 
tion, the engine configurational constraints restrict the direct-drive turbine 
diameter to a fraction of that of the fan, leading to low turbine wheel speeds. 
Consequently, the resulting fan turbine usually consists of a large number of 
highly loaded stages and/or compromised turbine efficiency. Alternatively, the 
fan may be run up to tip speeds substantially above those for which good effi- 
ciency has been achieved for the fan pressure ratios commonly used in high- 
bypass-ratio engines. The ultimate compromise has usually meant acceptance of 
the penalties associated with a large number of fan turbine stages, compromised 
turbine efficiency, or compromised fan efficiency, or some combination of these. 
Such penalties could be reduced or conceivably eliminated by use of a fan 
capable of high efficiency at high rotative speeds and the low pressure ratios 
compatible with high bypass ratios. The purpose of the current investigation 
is to do precisely this, to design for, and demonstrate, good efficiency for a 
h i gh-t ip-speed , relatively low- total-pressure-ratio, axial-flow fan stage. A 
design tip speed of 1600 ft/s (A88.6 m/s) and a stage total pressure ratio of 
1.50 were specified for these purposes. The principles to be examined are 

(1) Whether low losses and therefore good efficiency could be obtained 
by designing a transonic rotor for weak oblique shock waves in the 
hi gh-Mach-number tip region at design-point operation (thereby 
avoiding normal shock losses and the potential blade suction surface 
and casing boundary layer separation) 

(2) Whether adequate stall margin would be obtained, particularly in the 
presence of distorted inlet flow, with this type operation 

Achievement of the objectives of this investigation would permit much 
improved fan and fan turbine compatibility, thereby minimizing the penalties 
previously mentioned. Still further, the concept could reduce substantially 
the problems associated with the inherent hub and tip mismatch in work poten- 
tial from which all 1 ow-rad i us-rat io rotors (whether fan or multistage compres- 
sor inlet stage) suffer. Indeed, the advantages of being able to select the 
appropriate hub speed for a given multistage compressor first-stage pressure 
ratio without penalty to the tip region performance may be equally as important 
as the benefits of matching the high-bypass fan to the fan turbine. Further, 
the current interest in turbojets with very high flow per unit of frontal area 
dictates use of low-i nlet^radi us-ratio, hi gh-i nlet-Mach-number rotors. The 
ability to use the minimum structurally feasible radius ratio reduces the 
required inlet Mach number (and hence the required total axial diffusion) for 
a given specific flow. In addition, there may be the potential for reduced 
noise through the minimization of upstream shock strengths for the selected 
rotor tip speed. 
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The rotor used for demonstration was selected to be approximately 30 in. 
(0.762 m) in diameter, a size which permits accurate measurement with conven- 
tional instrumentation. The rotor blade inlet radius ratio of 0.50, aspect 
ratio of 2.5 to 3, and airflow per square ft of annulus area of kO to k2 lb/s 
(195-3 to 205.1 kgm/s) were selected as being representative of current air- 
craft propulsion fan design practice. 

This report outlines the design procedures and presents the aerodynamic 
and mechanical design results. 
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AERODYNAMIC DESIGN 


DESIGN CONCEPTS AND OBJECTIVES 


Several new rotor aerodynamic concepts influenced the design of the 
current rotor. A design objective was to avoid tip normal shock losses by 
substitution of oblique shocks in the hi gh-Mach-number tip region. In addition, 
the strengths of the oblique shocks were restricted to values for which the 
static pressure rise would not exceed the separation criteria generally accepted 
for turbulent boundary layers on flat plates. The static pressure rise and the 
total pressure ratio values were kept below the maximum values obtainable from 
rotors of 1600 fps (1+88.6 m/s) tip speed to make them consistent with near- 
optimum hub potential. Lower losses and an improved stall margin were antici- 
pated to result from this design approach. Another feature essential to this 
approach was the calculation of the rotor design outlet velocity diagrams that 
satisfied radial equi 1 ibri urn wi th supersonic rotor-tip relative leaving veloci- 
ties. In spite of the excess tip speed (for the design pressure ratio), the 
design value of rotor-tip static density rise was small enough so that axial 
velocity did not incur a significant reduction. Accordingly, the use of nega- 
tive camber at the rotor tip proved to be unnecessary for the radially constant 
stage pressure ratio. 

To terminate the necessary supersonic waves at the blade trailing edges, 
a characteristic procedure allowing for change in radius and stream filament 
thickness was used for the wholly supersonic sections. The locations of the 
points for expansions and compressions and their points of impingement on oppo- 
site surfaces were carefully controlled to cancel all oblique shocks inside 
the blade passage. The shocks were also separated to prevent the accumulating 
effect from causing separation. The weak shocks that began at the blade lead- 
ing edges were prevented from increasing (in the high Mach number outboard 
region) by utilizing small leading-edge thicknesses and wedge angles. At the 
design condition, no tip waves emanated upstream of the cascade. An important 
design consideration involved the magnitude of the speed at which supersonic 
flow could be initiated in the tip region (i.e., the starting Mach number). 

The contraction ratio was prescribed to be sufficiently small to permit tip 
starting at a speed somewhere below design speed. 


The principle design objectives were set by contract. Others were selected 
for convenience of procurement and testing, but were consistent with the basic 
objectives. The contract-specified parameters are listed below: 


Overall total pressure ratio 
Flow per unit annulus area 
Equivalent tip speed 
Rotor and stator aspect ratios 
Rotor inlet hub-to-tip radius ratio 
Rotor tip solidity 


1.5 

42.0 lb/s-ft 2 (205.1 kgm/s-m 2 ) 
1600 ft/s (488.6 m/s) 

3.0 (originally) 

0.50 max. 

1.4 (approx.) 
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The following parameters were specified by the designer for the final design. 


Objective adiabatic efficiency 
Equivalent total flow 
Rotor inlet hub-to-tip radius 
Tip diameter 

FLOW PATH AND VELOCITY 


0.86 

148 lb/s ( 6.71 kgm/s) 
ratio 0.462 

28.74 in. (0.73 m) 

DIAGRAM DETERMINATION 


General Procedure 


The subject fan design was initiated utilizing the procedure of the usual 
quas i- three- d i mens iona 1 compressor design approach. That is, the flow path and 
the blade inlet and outlet velocity diagrams were determined under the assump- 
tion of axisymmetric, inviscid flow utilizing the AiResearch stream filament 
computer program. In this approach the assumed axisymmetric surfaces are repre- 
sented in the radial-axial plane by streamlines (fig. 1). The radial component 
of the equations of motion is solved to obtain the radial variations of the flow 
conditions between blade rows. The mass flow is then summed or integrated along 
radial or near-radial lines from hub to tip using the previously calculated 
radially varying flow conditions. When the calculated mass flow at one axial 
position equals the prescribed flow rate, both continuity and radial equilibrium 
are satisfied at that station and the same calculation is performed at the next 
station. As in most such programs, after satisfaction of continuity at every 
compressor calculation station, a curve-fit routine is used to connect points of 
constant mass flow. These connecting streamlines then provide new streamline 
curvatures and slopes in the meridional plane which permit a repeat solution of 
the radial equilibrium equation with the new streamline geometry. Subsequently, 
a new set of streamlines is defined. The process may be continued until a com- 
patible flow field is defined within the initial assumptions and further itera- 
tion leads to no appreciable changes in streamline curvatures. In addition, a 
new AiResearch axisymmetric computational procedure using a finite difference 
solution to the radial component of the equations of motion with interior grid 
points was used. In general, the agreement was good. 

Both procedures permit arbitrary definition of the radial variations of 
enthalpy and entropy (or losses). In general, the radial variation may be 
specified for any two of the three quantities: pressure ratio, temperature 

ratio, and efficiency (or other equivalent). The bases for these specifications 
are empirical with due consideration being given to the objective aerodynamic 
operating modes. 


Rotor Annulus Convergence 

Because a high rotor aspect ratio (which was originally specified to be 
3.0) was contractually required, the projected chord of the blade was short, 
leading to relatively steep inner and outer casing slopes as illustrated in 
fig. 1. These slopes and the accompanying concave inner- and outer-wall con- 
tours upstream of the rotor resulted in a parabolic inlet absolute velocity 
distribution with low inner- and outer-wall values at the rotor inlet and a 
high velocity near the pitch. Conversely, at the rotor exit the convex wall 
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Radius, in. 







contour and th.e decreasing wall slopes resulted in a rotor exit meridional 
velocity distribution with low velocity near the pitch radius and high velocity 
on the walls. These meridional velocity distributions have pronounced effects 
(some beneficial) on the rotor configuration. 

The meridional velocity ratios (V j) resulting from the rotor inlet 

and exit distributions varied from about unity at hub and tip to approximately 
0.85 in the central flow region. The low values in the central flow region 
were believed to be tolerable because the diffusion factors and static pressure 
rise coefficients were conservative and the region was nol; adjacent to wall 
boundary layers. A matter that could be of greater poncern was the effect of 
the resulting mid-blade camber and its structural implications. Due to the 
modest work input (for this wheel speed) required of this rotor, the change in 
tangential velocity produced is relatively low. With a low value of the meri- 
dional velocity ratio 1 ( 1 a rge diffusion), the required change in tangen- 

tial velocity could be accomplished at an almost constant relative flow direc- 
tion with essentially uncambered blade sections. Also, since the mid-passage 
flow is to be transonic, thin airfoil sections are desired for low losses. It 
is difficult to obtain the desired mechanical stiffness for such blades. Hence, 
in hopes of achieving a satisfactory rotor blade mechanical design with a single 
part-span damper, an aspect ratio change from 3-0 to 2.5 was permitted by NASA. 

A detrimental effect of the rotor meridional velocity ratios imposed by 
the flow path convergence results from the relocation of the maximum static 
pressure ratio region to the center sections instead of at the tip. Because 
the rotor inlet relative Mach number normally increases with radius and the 
allowable shock static pressure rise is proportional tp the preshock Mach 
number (e.g., see refs. 1 and 2), the problem of shock- i nduced separation on 
the pitch sections becomes critical. The use of multiple shocks provides a 
theoretical solution. However, little information exists regarding the dis- 
tance between shocks that is required before the boundary layer is restored to 
a healthy profile after the shock interaction. In spite of this uncertainty, 
spaced multiple shocks were incorporated in the mid-passage region for the 
current design to reduce each shock static pressure ratio below the value 
allowable for a single shock at the mid-passage Mach number. 

ROTOR DESIGN REGIONS 
Method of Classification 

The output from the axisymmetr ic, rad ial-equ i 1 i l?r i urn program yielded three 
types of rotor velocity diagrams based on rotor design-point inlet and outlet 
Mach numbers. Each of these three velocity diagram types characterizes a blade 
design region requiring a different section design procedure. These are sum- 
marized below. 



Inlet relative 
Mach number 


Outlet relative 
Mach number 


Approx flow division at 
design point, percent 


Outboard 
Centra 1 
I nboard 


Superson i c 
Superson i c 
Subson i c 


Supersonic 
Subson i c 
Subson i c 


35 

45 

20 


7 



/ 


To define a single, structurally feasible blade composed of the three 
separately designed blade types listed above, new computational programs were 
written. The al 1 -supersoni c outboard sections can be completely designed 
analytically with the aid of several assumptions. The subsonic-entry, 
subsonic-exit hub sections can be derived f rom wel 1 -establ ished empirical 
incidence and deviation angle rules and a selection of a meanline shape 
(circular arc, multiple circular arc, polynomial, etc.). The inlet region 
of the mid-span transonic sections can also be calculated. However, the 
terminal shock configuration and location must be determined by somewhat less 
formal methods for the conventional transonic sections. Still, the blade 
design program described "herei n did yield transonic sections that gave a 
smooth aerodynamic transition and a structurally satisfactory stacking between 
the supersonic tip and subsonic hub regions. 

Theoretically, a blade section is defined as the contour generated by 
the intersection of the rotor blade surfaces and the appropriate axisym- 
metric stream surface. From the picture of the final axisymmetric stream 
surfaces shown in fig. 1, they may be seen to be very closely approximated 
by conical surfaces through the blade leading- and t ra i 1 i ng-edge streamline 
intercepts. Therefore, for our purposes the blade sections referred to 
hereafter will be on these conical surfaces. 


Since the outboard section has supersonic inlet and outlet relative veloc- 
ities, an approximate method of characteristic was used. Because the principle 
is more easily understood for the two-dimensional problem, the application of 
the approximate method of characteristics is first illustrated for the two- 
dimensional flow field, and later generalized for the quasi-three-dimensional 
flow field. 

In the following presentation, the equat ions-' and formulas define relation- 
ship's among angles and other geometrical parameters, where the angles are cor- 
rectly expressed in radians. However, illustrative calculations and diagrams 
present angles in degrees for easier interpretation. 


' Outboard Section Design Theory (2-D) 

An idealized two-dimensional (2-D) section is shown in fig. 2. As shown 
in this figure, a simple 2-D section to meet the design conditions would have 
reverse camber. Let the leading- and tra i 1 ing-edge shock strengths (in terms 
of deflection angle) be denoted by and 6 2 , respectively, as illustrated in 
fig. 2. Then the relation between the relative flow angle P u upstream of the 
leading-edge shock and the flow angle g downstream of the tra i 1 ing-edge is 
given by the expression 

^2 = 8u " 61 + 6 2 0) 

Similarly, for weak shocks the relation between the upstream and downstream 
Prandtl -Meyer expansion angles v u and ^2 < s closely represented by the expression 



( 2 ) 
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S-72475 

Figure 2. --Sample Two-Dimensional Wave Pattern Illustrating Design Logic 
for the Compressions Contained in Leading-and T ra i 1 i ng-Edge 
Shocks . 


( 3 ) 


The sum of eqs . (l) and (2) yields the leading edge shock strength 

6 1 = ^u " ^2 

where 

v + P 
* “ 2 

As a consequence of the fact that the leading-edge shock has a greater slope 
than a corresponding upstream Mach (or characteristic) line, it follows that 
there is some point on the inlet line of the blade (denoted as point 1 in fig. 
2) whose negative characteristic line intersects the blade surface at point u. 
For a steady supersonic reversible adiabatic flow of a perfect gas with con- 
stant specific heats in the absence of body forces, the two-dimensional theory 
of characteristics states that the value of the characteristic parameter § 
remains constant along the negative characteristic, independent of the manner 
in which the characteristic line may weave through the flow field. Hence, it 
follows that 


? u - q < 4 > 

Since the flow conditions are uniform along the inlet line of the blade, it 
follows that the value of § is determined by the section inlet flow condition, 
i .e . , 


5 


u 


V 1 + P 1 


(5) 


In an analogous manner, the value of ' s determined directly from the discharge 
flow condi tion, i .e. , 



v 2 + e. 


( 6 ) 


It now follows from eqs. (3) > (5), and (6) that the leading-edge wedge angle of 
the simple section shown in fig. 2 is uniquely determined by the upstream and 
downstream flow conditions for the hypothesized blade and wave configurations. 
As an example, for the section shewn in fig. 2, = 14.86, = 70 deg, = 

3.56, and = 73«6 deg; hence, = 42.43 deg and §2 = 38.58 deg; and there- 
fore, the leading edge wedge angle is 3-85 deg from eq . (3). 

Development of the tra i 1 ing-edge shock-strength begins by subtracting 
eq. (2) from eq . (1) at the rotor blade trailing edge: 


where 




v - p 
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The value of ' s determined directly from the downstream flow conditions 
as 


5 

2 



( 8 ) 


p -J- — 

The value of 5 can be determined from the value of § and the blade surface 
u u 

angle at point u as 

§ + = - p 

u u u 


or from eq . (4) , 


5 U - 5, - S> u 


(9) 


Hence, it follows from eqs. (7), (8), and (9) that for the configurations 
assumed above, the tra i 1 i ng-edge shock strength is determined by the inlet and 
discharge flow conditions and by the blade surface angle just upstream of the 
leading-edge shock. 

For the section shown in fig. 2, = -35-02 deg, = 42.43 deg, and 

P u = 8 i = 70 deg, hence, 6 2 = 7.45 deg. 


Tra i 1 i ng-edge shock strength effect .--A1 though the leading-edge shock 
strength is uniquely determined by the inlet and discharge flow conditions, 
some degree of freedom is available to the designer relative to the choice of 
trai 1 ing-edge shock strength. The effect of the trai 1 ing-edge shock strength on 
the section aerodynamic and mechanical properties is discussed below. 


As a consequence of eqs. (7), (8), and (9), it follows that 


6 2 + 


( 8 . 




( 10 ) 


The left hand side of eq . (10) represents the total suction surface compres- 
sion from the first captured Mach wave to the blade trailing edge. Since the 
right hand side of eq . (10) is determined by the inlet and discharge flow condi- 
tions, it follows that the total net suction surface compression is fixed. 

The wave pattern for the nominal trai 1 ing-edge shock strength, defined as 
the shock strength corresponding to = P u (no suction surface compression or 

expansion forward of the final shock), is illustrated in fig. 2. For this case, 
the net required turning is simply the difference between the tra i 1 i ng-edge and 
» leading-edge shock strengths. For tra i 1 ing-edge shock strengths below nominal, 
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Figure 3."Sample Wave Pattern Showing 1.85 Degrees of the Trail ing- 
Edge Shock Shifted Part-way Forward. 
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Figure 6. --Blade Section Using Distributed Suction Surface Compressions 
to Weaken the Tra i 1 i ng-Edge Shock Below the Leading-Edge 
Shock Strength Resulting in an Imaginary Blade Aft Section. 



additional suction surface compression must be provided, either as a finite wave 
emanating from the suction surface (figs. 3 and 4), or as distributed turning 
(fig. 5) • For the type of blade section shown in fig. 2, where the leading-edge 
shock terminates at the adjacent-blade trailing edge, the minimum tra i 1 ing-edge 
shock strength is established by mechanical considerations. For the section of 
fig. 2, the trai 1 ing-edge wedge angle 6^.^ is 


6 


TE 



( 11 ) 


Accordingly, 6^.^ vanishes as 6^ approaches 6^ and becomes negative for < 
thereby yielding a nonphysical section (fig. 6). 

The trai 1 ing-edge shock can be made stronger than nominal by offsetting 
the excess compression with an equivalent distributed expansion on the suction 
surface between the points u and I. 


Aerodynamical ly , changing the trailing edge shock strength affects the 
blade surface Mach number distribution, and hence indirectly the blade viscous 
and profile losses. In addition, as a result of raising or lowering the mean 
Mach number level in the passage portion of the blade, the trai 1 i ng-edge shock 
affects the ability of the section to start. 

Weakening the t ra i 1 i ng-edge shock results in a reduction of the average 
surface Mach numbers, and hence should result in a reduction in viscous losses. 
In addition, distributing the overall blade diffusion and using the weaker 
t ra i 1 ? ng-edge shock should reduce the tendency of the boundary layer to 
separate and eliminate the profile losses associated with separation. On the 
other hand, section changes that weaken the t ra i 1 i ng-edge shock tend to reduce 
the mean passage Mach number and, hence, will make the section more difficult 
to start. This latter effect is partially offset by the reduced Mach number 
into the bow shock in the unstarted operating condition where upstream com- 
pressions usually precede the cascade. 

Using similar reasoning, it is anticipated that section suction surface 
changes that strengthen the tra i 1 ing-edge shock will increase the blade 
losses and, at the same time, make the resulting section easier to start. 

Solidity effects . --A section of higher solidity than is portrayed in 
figs. 2 through 6 will cause the leading-edge shock to impinge on the suction 
'surface. Cancellation of this wave by an expansion corner decreases 3 on the 
suction surface, downstream of the expansion kink, by the magnitude 6^. 

Denoting this value as 3 yields, in general, 


From eq. (1 ) , 


6 = 3-3 

TE 2 t 


B u >e 2 + 6, - 6 2 
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therefore , 


3 


t 



- 6 


2 


and, finally, for the high-solidity section of fig. 1 , 


6 


TE 


= 6 


2 


( 12 ) 


For this type of section, the trai 1 ing-edge shock strength 6 can evidently be 

less than for the section shown in fig. 2 before the trai 1 ing-edge region 
becomes excessively thin. 

It is apparent from a comparison of figs. 2 and 7 that the section • 
solidity can significantly influence the section mechanical properties. In 
particular, increasing the section solidity tends to increase the section 
thickness in the tra i 1 i ng-edge region. As will be seen later, three-dimensional 
sections have higher than minimum solidity in order to take advantage of these 
effects . 


From an aerodynamic point of view, the increased solidity (as a result 
of its effect on the trai 1 ing-edge wedge angle) allows a weaker t ra i 1 i ng-edge 
shock to be used than would otherwise be possible. In addition, it allows 
some distance for the suction surface boundary layer to energize after its 
I encounter with the leading-edge shock. This latter effect should tend to 

reduce profile losses, and may more than offset the increase in viscous losses 
due to the increase in section wetted surface. 

Outboard Section Design Theory (3-D) 

The three-dimensional (3-D) characteristic procedure used in the actual 
design of the rotor outboard sections is basically analogous to the 2-D method 
just described. The leading-edge shock strength may still be defined in terms of 
the characteristic parameters, eq. (3). However, because of the effects of the 
stream tube contraction, the radius change, and the blade row rotation, eq. (4) 
is no longer valid and . The procedure by which the actual outer blade 

sections are derived utilizes the simple two-dimensional character ist i c equat ions 
modified to reflect changes with axial distance in both stream tube thickness 
and radius. With this type of flow field, the two-dimensional variations in 
wave strength may be altered or even change signs as they move in the axial 
direction. That is, the variations are functions of stream tube contractions 
(or expansions) as well as blade turning. 

The actual 3-D sections contain modifications which make them consistent 

with: 

(1) The effects of blade row rotation, radius changes, and variations 
in stream tube height on the flow field 

(2) Certain necessary mechanical constraints such as leading- and 
tra i I ing-edge radii and blade section thicknesses 

(3) The calculated shock and viscous losses 
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Figure 7-~~Sample Wave Patterns Showing Leading- and Tra i 1 ing-Edge 
Shock Cancellation with Increased Solidity Which Allows 
Separation of Leading- and Tra i 1 i ng-Edge Shock. 


(4) Blade passage area variations which permit a reasonable passage 
"starting" speed (sufficiently below design speed) 

(5) Adjacent sections of the same and different types 

Characteristic equat ions . --As developed in Appendix A, the theory of 
characteristics gives the differential dC at any point on the negative character- 
istic joining the points 1 and u, fig. 2. 

dS' = e" (r, a, ou, v, p) (13) 

where r, a, and ^ ^ are all functions of the meridional location (m) of the 

point on the characteristic. Since uj is a design constant, eq . (13) may then be 
wr i tten as 


d§’ = e" (m, v , P) dm (14) 

where the e-function is defined in Appendix A by the expression of eq . (A-22) . 

In general, the v and 3 distributions on the negative characteristic 
joining the points 1 and u are not known a priori and depend upon the 
detailed contouring of the blade suction surface. H^nce, the exact relation- 
ship between and ^ cannot be determined without a detailed calculation of 

the intervening flow field. This calculation would be costly and time con- 
suming. However, examination of the function &2 indicates that under certain' 
conditions it may be closely approximated by substituting the inlet valyes v t 
and P^ for the local values of v and P. Then 

d|' « e~ (m, v 1 , pp 
or 

/ m 

u •; c. v 

m. * 


dm 


I . P,) dm 


(15 

( 16 ) 


In addition to the above approximations the conditions at the cascade 
exit (point 2) are assumed to be equal to t!:e core conditions as calculated 
from the trailing edge mixing analysis (Appendix B) . Further, the inlet con- 
ditions (point 1) are taken to be the inlet conditions after the leading edge 
blockage calculation described in Appendix C. The final equation for the 
leading edge shock strength in three-dimensional flow, analogous to eq. (3) 
for two-dimensional flow, is then 


6 

1 


1 ’ corr 


m 

P 

m. 


2 (nv 


1 corr 


P 1 corr^ 


dm 


_ l 


2, core 


(17) 
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where the first two terms on the right represent § u , the subscript corr denotes 
the cascade inlet conditions corrected for leading-edge blockage, and the sub- 
script core denotes inviscid flow conditions corrected for blockage. 

As a consequence of eq. (17), it follows that the leading-edge shock 
strength is uniquely determined in a manner which is essentially the same as 
the technique explained for the two-dimensional sections. For fixed inlet and 
outlet velocity diagrams (fixed CO rr an< ^ ^2> core)* t * 1e leading-edge shock 
strength can be varied without altering the cascade wave configuration only by 
changing the integral term in eq. (17) • This, in turn, can be done only by the 
blade annulus area distribution (stream tube height). This procedure allows 
some design freedom for altering the blade thickness. 

Using the values of corr anc l Pi corr as calculated in Appendix D and 
values of M core and 3 C ore from Appendix B, the leading-edge shock strength may 
be calculated as indicated in the following illustrative example. 


M 


1 


corr 


1.502 


P 


1 


corr 


= 62. 17 deg 


H = 1.143 

core 


P = 59-35 deg 
core 3 


The corresponding values of and § are: 

r 1 corr core 

§’ = 37.07 deg 

1 co r r 3 


= 30.79 deg 
core 3 


The value of the integral term in eq . (17) is 


m = 1 . 787 



m = 0 


P) 


dm = -1.94 deg 


The corresponding value of leading-edge shock strength as calculated by 
eq . (17) is therefore 

= 37.07 - 1-94 - 30.79 


or 


6 = 4. 34 deg 
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Accuracy evaluation of approximate e function . — I n the design of the 2-D 
section, the suction surface from the leading edge to tfip first captured Mach 
line was made straight in order to eliminate waves which would escape upstream 
of the rotor. The analogous procedure for 3 _ D sections consists of defining 
the suction surface from the leading edge to the first captured Mach line to 
be coincident with the 3 _ D free streamline. This streamline is the locus of 
the path which the flow would take in the absence of blades if it entered the 
blade annulus axi symmetr i cal ly with the blade row, inlet values of Mach number 
and flow angle and if it was restrained to flow within the same axisymmetric 
stream tube geometry for which the blade section is being designed. This free 
streamline may be referred to as the no-work streamline, and its flow field 
represents an exact solution of the fluid flow eqqations defined for the free 
streamline constraints. As a consequence it can be used as a partial check of 
the validity of the approximate eq. (15). 

i Using the actual rotor design defined herein, a comparison of the charac- 
teristic parameter ^"calculated for the free streamline flow field and from 
eq . (16) is shown in fig. 8. As can be seen, the comparison is excellent up 
to the point of free streamline choke. 

Description of 3~D methods . --The conical development of the passage 
between two adjacent blade sections associated wi th streaml i ne number k (a repre- 
sentative section) is shown in fig. 9* Similar developments for streamlines 1, 6 
9 and 12 are shown in figs. 10a through lOd. The grid shown in the figure is the 
conical development of lines of constant Z (lines A-A) and lines of constant 0 
(lines B-B) on the conical surface. For purposes of discussion, the significant 
points on the section are labeled f rom(j) to(J), as fpllows. 


Poi nt 
Number 

© 

© 

© 


Description 

Suction-surface, leading edge radius tangency point 

Intersection of first captured Mach line with suction surface 

Point just upstream of the intersection of the leading-edge shock 
with the shock surface 


© 


Point just downstream of the intersection of the leading-edge 
shock with the suction surface 


© 

© 

© 

© 

© 


Suction surface trai 1 i ng-edge radius tangency point 

Pressure surface t ra i 1 i ng-edge radius tangency point 

Point just downstream of the intersection of the trai 1 i ng-edge 
shock with the pressure surface 

Point just upstream of the intersection of the t rai 1 i ng-edge shock 
with the pressure surface 

Pressure surface leading-edge radius tangency point 


21 





Figure 9 • --Representat i ve Rotor Section (Streamline No. 4) Showing 
Pertinent Design Blade Segments and Points. 
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The line running across the passage from point (§) to point (J) represents 
the leading-edge shock, and the line running from point © to point (?) repre- 
sents the tra i 1 i ng-edge shock. The dashed line running from point © to point 
represents the first captured Mach line. The dashed lines running across 
the interblade passage are used during the calculation of the blade passage 
area distribution. The dashed line surrounding the blade section represents 
the section envelope which was calculated in the design process. The solid 
line within the envelope represents the blade surfaces which are obtained by 
subtracting the boundary layer displacement thickness from this envelope. 

The conical development of the section envelope can be conveniently divided 
into five steps dealing with the following segments. 

Step Segment of section 

number envelope developed 


1 

Poi nt 

© 

to 

po int 

© 

2 

Point 

© 

to 

point 

© 

3 

Point 

© 

to 

po int 

© 

4 

Po i nt 

© 

to 

point 

© 

5 

Point 

© 

to 

po i nt 

© 


The techniques associated with each step segment are presented below. 

The first step covering segment © - © comprises calculation of the 
free streamline or n no work streamline 11 . The initial values of Mach number 
and flow angle used for the no work streamline calculation differ from the 
axisymmetric program output by the blockage effects associated with the section 
leading-edge radius and the detached portion of the leading-edge shock. The 
equations used to calculate the magnitude of the changes are developed as indi- 
cated in Appendixes C and D. The net effect of these corrections is as follows: 


where I and i^- are given by eqs. D-8 and D-9 of Appendix D, respectively. 
I n add i t ion , 


M. = M(v.) 
1 1 ' 


where v. , the Prandtl -Meyer angle, is calculated from 

v i - v i + <s i - s i> 

The section envelope angle at point © is set equal to (3. 
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It was found to be most practical to assign i^. as an input item to the 

computer program developed to implement this design procedure. The program 
then calculates an exact value of i^. for use in the next iterative design pass. 

Since (3. is assigned to be the suction surface envelope angle 3^ it can 

be seen that the suction-surface incidence angle with respect to the cascade 
entrance conditions is zero. With respect to the free-stream inlet flow, the 
leading-edge suction surface incidence angle is i^ + i^. 

In addition to defining the forward portion of the section envelope, the 
free streamline calculation also defines the flow fi^ld in the forward region 
of the section. The location of point ( 2 ) on the section envelope is determined 
by constructing a Mach line from the pressure surface leading edge of the adja- 
cent section through the known free streamline flow field until it intersects 
the section envelope suction surface. The section envelope angle at point ( 2 ) , 

P/’-jkjis equal to the free streamline angle at that point. 


The portion of the section envelope suction surface developed in step 2 
consists of two smooth curves between points ( 2 ) - (|) and ® - ( 5 ), separated 
by a discontinuity between points (J) and (?) • For simplicity, the rate of 

change of 3” with respect to meridional distance on these two smooth surfaces 

,33 


was held constant (-r- = K) . This rate, defined as K, is 
* om 

- B%) - (»fo - B%) 


K = 


m 


® ‘ m © 


( 18 ) 


The angle is determined from step 1. The angle 3^ is calculated as 

the difference between the discharge flow angle of the inviscid core and the 
flow deflection across the trailing edge shock, ie. , 


_ ' wC 

P (5) S core 


09) 


The technique used to calculate S is described in Appendix B. The 

n core rr 

quantity 6 ^ is calculated from the expression 


= &i + A6 


( 20 ) 


where 6^ is the leading-edge shock strength, and A6 is the des igner- specif ied 

difference between the trailing- and leading-edge shock strengths. The tech- 
nique used to calculate 6^ was described in the preceding description of the 

outboard 3-D section design procedure. 
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For the current designs the magnitude of the leading-edge shock is assumed 
to be constant across the passage with complete cancellation at point ® ■ As 
a consequence 

% ' e © - 6 , < 2 ’> 

v v 

The value of p at any meridional location on the section envelope 
between points © and © is given by 


P ,W (m) = K( 


® + K.(m - m ) 

© a 

3 (m) = P g) + K(m - m@) - for m > 

and the corresponding angular location of the envelope is given by 


for m < 
6 


( 22 ) 



(23) 


The techniques used to determine the location of the leading-edge shock, and 
subsequently the value of nru , are presented in Appendix C. 

Kds 

In the blade layout, step 3 is simply an angular translation at the 
trailing edge of the section envelope to account for the blockage effects of 
the suction and pressure surface boundary layers and the physical blockage of 
the t ra i 1 i ng-edge radius. It is apparent from fig. 9 that a significant 
percentage of the section envelope thickness arises directly as a consequence 
of the tra i 1 i ng-edge blockage. The magnitude of the trai 1 i ng-edge blockage 
is evaluated by means of a mixing loss analysis which is described in Appendix B. 


For a 2-D section, the trailing edge region between poi nts © and© woul d 
be straight in order to eliminate wave disturbances downstream of the section. 
Correspondingly, in the 3 - D design, the pressure surface from points (6) to 
©is designed so as to be coincident with the free-streamline path correspond- 
ing to the cascade outlet inviscid core flow conditions and passing through 
point (6) so as to eliminate wave disturbances downstream of the cascade. In 
addition, the free streamline analysis provides a detail description of the 
core flow field in the vicinity of the trailing edge. This information, along 
with the knowledge of the strength of the t ra i 1 i ng-edge shock, is sufficient 
to determine the approximate location of the shock in the flow field. The 
intersection of the free streamline and the trai 1 ing-edge shock from the adja- 
cent blade determines the location of point©. 


In theory, knowledge of the pressure-surface flow conditions downstream 
of the leading-edge shock^nd upstream of the trai 1 i ng-edge shock, along with 
a reasonable assumption of the Mach number distribution (e.g., linear) between 
points © and ®, is sufficient to allow the pressure surface contour to be 
extended analytically from point (D to point ® . An assumed Mach number dis- 
tribution yields a corresponding 3-distribution, from which the polar coordinates 
can be computed by numerical integration from point (§) to point © . However, 
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the resulting value of 0^ would not necessarily agree with the previously 

calculated value. In view of this fact, this approach to completing the section 
envelope was not employed in the present design method. Instead, the angular 
location of each point on this surface was assumed to be a cubic function of 
the meridional distance, m, so that 


0 


3 2 

am + bm + cm + d 


(24) 


The location of the surface at points (§) and (§) is known as a result of 
the previous construction steps. The corresponding slopes at these points are 
also known since 


B ® = B © ' 6 i 

(25) 

B ® “ B © ‘ 8 k 

(26) 


where , the leading-edge envelope wedge angle, is equal to 6 • and 6 the 

I K 

pressure-surface kink angle, equals 6^, the tra i 1 ing-edge shock strength (eq. 20)). 
Hence, the four coefficients of eq. (24) can be easily evaluated. 


It might be noted at this point, that in the previous construction steps it 
was implicitly assumed that the flow deflection across the leading edge shock 
was the same at the suction and pressure surfaces. A similar assumption was 
made relative to the trailing edge shock. The validity of this assumption is 
related directly to the shape of the surface from point (8) to point (^) and any 
consequent waves from that surface. For the small changes in (3 that take place 
along the subject surfaces, the pressure surface compression waves are not 
expected to affect significantly the leading-edge and trai 1 ing-edge shock wave 
profiles. As stated previously, the suction-surface rate of change of 8** with 
respect to the meridional coordinate is constant. On the pressure surface, eq. 
(24) means that 


d(ta " km 

dm 

Accordingly, for small changes in g , the two surfaces are similar in contour 
and an appreciable degree of wave cancellation should result. 

In order to calculate the blade section from the section envelope it is 
necessary to choose a distribution of boundary layer displacement thickness on 
the section surfaces. In the present design procedure, the boundary layer dis- 
tribution on the suction and pressure surfaces is assumed to be a quadratic 
function of meridional distance. The three quadratic coefficients are defined 
by assuming that the displacement thickness at the leading edge is zero; and by 
specifying the rate of growth of the boundary layer at the leading edge, and the 
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magnitude of the displacement thickness at the trailing edge (see sketch below). 
The sum of the pressure and suction surfaces displacement thicknesses is evalu- 
ated using the mixing calculation described in Appendix B. The magnitude of 
the individual thickness was determined by assuming the suction surface thick- 
ness to be 20 percent larger than the pressure surface thickness. For simpli- 
city, the rate of change of boundary layer displacement thickness at the lead- 
ing edge is specified' in terms of the angle between the section envelope and 
the metal surface. 



Meridional distance 
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One approach to specifying the initial rate of boundary layer growth would 
be to evaluate the displacement thickness at the location of the first captured 
Mach line, and then to choose an initial growth rate consistent with this value. 
This approach was not used in the development of the present design; rather, an 
ini.tial growth angle of 0.2 deg was assumed for each blade section. 


Summary of outer-section design sequence . --The input to the outer section 
design procedure can be divided into two categories: (1) the input supplied 

by the axisymmetric program, and (2) the input supplied by the designer. The 
former includes the inlet and discharge flow conditions, the streamtube geome- 
try, and the axial chord of the section. The primary and secondary input sup- 
plied by the designer is listed in the following table. 


P r i ma ry 

i nput Description 


A6 


Secondary 
i nput 



L.E 


T.E. 


a 


O' 


H 

f 


Difference between trai 1 i ng-edge and leading-edge 
shock strength 

Pressure surface kink angle 


The profile loss factor to be used in the mixing 
cal culat ion 

Suction surface incidence factor (see Appendix D) 

The leading-edge radius 
The trai 1 ing-edge radius 

The suction surface initial boundary layer growth angle 

The pressure surface initial boundary layer growth 
angle 

The t rai 1 i ng-edge boundary layer shape factor 

The ratio of the suction surface to total boundary 
layer displacement thickness 


The secondary input parameters are those quantities which are generally 
held fixed during the design process. The only primary input variable with 
which the designer can influence the outboard section properties is A6. The 
choice of this variable directly affects the section maximum thickness and the 
minimum passage area. If the maximum section thickness is specified from mech- 
anical design considerations, then the value of A6 and consequently the minimum 
passage area are determined. If the value of minimum passage area is not 
acceptable then the solidity, i.e., the axial projection of the section must 
be changed. The design process can be speeded up considerably if the axial pro- 
jection is treated as a primary design parameter, with the condition that the 
axial projection value used in the final axisymmetric pass and the final blade 
design pass correspond. 
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Inboard Section Design Theory 


The inner rotor sections, by definition, have subsonic inlet and discharge 
Mach numbers. As a consequence, the analytical techniques used in the design of 
the outer sections are no longer valid and an alternate design approach must be 
used. One of the main considerations in choosing the inner section design tech- 
nique was that the method should provide, automatically, a simple means of fair- 
ing the inner and outer sections in a smooth, continuous, and structurally 
acceptable manner. For this reason, a design technique which was identical in 
structure to that used for the outboard design was selected rather than the 
simpler and more conventional mu 1 1 ip 1 e- ci rcula r-a rc approach. 

The outer section was developed in six basic steps, two steps to develop 
the suction surface, two steps to develop the pressure surface, and one step 
each to develop the section envelope blockage and the blade metal surface. 

The inner section is developed in four steps, one each for the pressure and 
suction envelope surfaces and one each for the blockage and metal section 
development. The suction and pressure surface steps are individually divided 
into two steps each so that the final section is consistent with the six-step 
outboard section. Consequently, construction steps 1 and 2, and steps k and 5 
are considered together in the following description. The manner in which the 
inner and outer sections are faired smoothly into one another is discussed in 
the description of the central section design technique. 

Description of Inboard section design methods . --Steps 1 and 2 encompass 
the generation of the section envelope suction surface. For the inner sections, 
the surface angle variations along this envelope are defined by 


dp** 


K(1 



(27) 


The parameter X allows the designer to control the distribution of suction 
surface curvature. For X = 0, the distribution of curvature is the same as 
was used for the outer sections. For X > 0~, the surface curvature tends to 
shift towards the trailing edge, and for X< 0, the surface curvature tends to 
shift towards the leading edge. Integrating eq. (27) from the leading edge 
(m =0) to a general point on the surface yields 

e»* - = Km 0 ♦ (28) 


where the subscripts identify locations designated in fig. 9. Eq. (28) can 
be used to evaluate K in terms of the change in section envelope surface angle 
from the known overall change in p**t 


. JL JU JL 


k, ®, to , 

' m (1 + X/2) 

© 
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Accordingly, eq . (28) becomes 


J. J. 


X m \ m 


P** = P** 




@ m Q 


(29) 


It will be noted that X = 0 results in 


0 + j) 




B ** , 9 »* * p (t) ' 6 <fD 


I 


which displays the same form as eq. (22). 




The suction surface angle f3^ at the leading edge is evaluated from 



2 r 


L.E. 


0 Y cos 3 00 


(30) 


The second and third terms on the right of eq . (30) are analogous to the 
leading-edge thickness and shock blockage terms used in the outboard section 
development. The leading-edge thickness blockage term has been chosen so as 
to just offset the loss in blade passage area for a two-dimensional cascade 
with a straight suction surface. The shock blockage term, now more correctly 
referred to simply as an incidence correction factor, is an input quantity in 
the inboard section design procedure. It is by means of this term that section 
incidence is controlled. 


from 


The section angle (e.g., 


f^,fig. 9) at the trailing edge is evaluated 


(3** = 8 - ( 31 ) 

© C ° re 5 ** 

where 3 core is the inviscid core flow angle, and is the section envelope 

suction surface deviation angle. The suction surface deviation factor is 
eval uated as 


6** 

s 





(32) 


where 6 and 6 YE are two parameters supplied by the designer, the former being 
the deviation angle and the latter the t ra i 1 i ng-edge wedge angle of the section 
envelope. The reason for specifying 6** and6y* instead of 6** directly is that 
these two quantities are closely related to the corresponding quantities of the 
metal section and it is convenient to be able to control these quantities 
independently during the design process. 

Once the parameters iy, ° , and X are specified, the section envelope 

at each point can be calculated from 
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m 


(33) 


-JB<- ** 

P = P 1 + 


dP 




dm 


dm 


and the corresponding section envelope angular location from 


•a-* 

e 



** 


tan P 
r 


dm 


(34) 


With the suction surface specified by eqs. (33) and (34), only the posi- 
tion of the point corresponding to the outboard section point (3) and (4) on the 
surface in a manner consistent with the outer section remains. In the outer^ 
section development, this position is determined by extending a Mach line (the 
first captured Mach line) from the leading edge of an adjacent blade to the 
suction surfaces. For sections toward the hub (lower relative Mach numbers), 
this Mach line tends to approximate a leading-edge passage normal. For the 
inner sections, where such a Mach line does not exist, the leading-edge normal 
is used in its place. 


The third step for the inner sections is identical to step 3 for the outer 
section. It specifies the location of point (6) on the section envelope via the 
angular translation of the coordinate to account for the total effective blade 
and displacement thickness blockage. Steps 4 and 5 encompass the definition 
of the section envelope pressure surface. The angle of this surface at the 

•»H<- ** 


leading edge was previously evaluated in eq. 

■iH* 


( 25 ) as - a, 


where 


3 

CD 


is the angle of the suction surface as calculated in step 1 and 6^ 


is the section envelope leading-edge wedge angle specified by the designer. 
The angle of the pressure surface at the trailing edge is evaluated as 


•iHf- 

* = & - 6 

© core P 


(35) 


where 6 is the section envelope pressure surface deviation factor. This 
P 

quantity is evaluated from the section envelope deviation angle, and the 
section envelope trailing edge wedge angle as 


6 =6 
P 



(36) 


Eqs. (25) and (35) serve to define the pressure surface angle at the leading- 
and tra i 1 ing-edge points. The location of the leading-edge point is calculated 
from the location of point (Y) and the designer's choice of leading-edge radius; 
the location of the tra i 1 ing-edge point is determined from step 3* With the 
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location and slope of the surface specified at the leading and trailing edge, 
the remainder of the surface is defined by a cubic curve, i.e., 

•iHt 3 2 

0 = am + bm + cm + d (37) 

Step 6 encompasses the generation of the metal section within the section 
envelope. This operation is performed in the same manner as for an outer 
sect ion. 

Summary of inboard section design sequence. --The design of an inner rotor 
section using the techniques described above requires two types of input: 

(1) input supplied from the axisymmetric program, and (2) input supplied by 
the designer. The input supplied by the axisymmetric program includes the 
inlet and discharge flow conditions, and the streamtube geometry, and the axial 
length of the section. The primary and secondary input supplied by the designer 
is listed in the table below. 


Descr ipt ion 

Suction surface curvature factor 
Suction surface incidence factor 

Envelope deviation angle with respect to envelope mean line 
Section envelope tra i 1 ing-edge wedge angle 

Section envelope leading-edge wedge angle 


Profile loss factor for use in mixing calculation 

Tra i 1 i ng-edge boundary layer shape factor 
Suction surface initial boundary layer growth angle 

Pressure surface initial boundary layer growth angle 

Leading-edge radius 

Tra i 1 i ng-edge radius 

The primary input values represent those quantities which are varied by 
the designer in the process of obtaining an acceptable metal section. The 
secondary input quantities, in general, are fixed during the iteration process. 
The manner in which the primary quantities tend to affect the section proper- 
ties are discussed in the following paragraphs. 


Pr imary 
input 


f 


TE 

6 ** 

1 

Secondary 

input 


H 

a 

c 

a 


L.E. 


T.E. 
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The leading- and t ra i 1 i ng-edge section envelope wedge angles control the 
magnitude and location of the section maximum thickness. Increasing either 
wedge angle tends to increase the maximum thickness. Increasing the leading- 
edge wedge angle tends to shift the maximum thickness forward, and increasing 
the trai 1 i ng-edge wedge angle tends to shift the maximum thickness aft. The 
section envelope deviation factor controls the metal section deviation angle. 
The suction surface curvature and incidence factors control the blade passage 
area distribution. In actuality, there is some cross -coupl i ng between each 
input parameter and all the section properties. This cross -coup 1 i ng is in 
general smal 1 . 

The section passage area distribution can be described quantitatively in 
terms of the value and location of the passage minimum area, and qualitatively 
in terms of how quickly the area distribution opens up downstream of the mini- 
mum value. In general, it was found that shifting the suction surface curva- 
ture aft tended to shift the passage minimum area in the same direction, and 
increasing the incidence factor tends to increase the value of the minimum 
area. In the development of the hub section, these two factors were chosen so 
as to achieve a minimum passage five percent larger than the critical area at 
the passage inlet, with a continuous increase in area downstream of this loca- 
tion. The hub deviation factor was chosen to achieve a section deviation angle 
of 10 deg, and the leading- and t ra i 1 i ng-edge wedge angles were chosen so as 
to achieve a t /C of 0.08 located slightly aft of the mid-chord position. 


Central Section Design Theory 

The central sections, by definition, have supersonic inlet and subsonic 
discharge relative Mach numbers. At present no theoretical technique is readily 
available for analyzing entire sections of this type. In developing the central 
sections, consideration was given to the following: 

(1) Achievement of a smooth transition from inboard to outboard sections 

(2) Achievement of adequate choke margin 

(3) Achievement of adequate section thickness 

(A) Avoidance of excessive suction surface Mach numbers in the leading- 
edge region 

A main requirement of the central section design technique is that it be 
structured in a manner which is consistent with item (l) above, and have suffi- 
cient flexibility to allow the designer to satisfy items (2) through (A). In 
order to satisfy item (1), the middle section is developed in six construction 
steps . 


Central section design methods . --In the following paragraphs, the 
techniques used in each step are discussed in detail. In certain instances 
where a nonnumerical description of a particular step is inconvenient, a finite 
difference description of the step is presented. The finite difference proce- 
dures are basically those that are used in a computer program that was developed 
to automate the blade generation technique. 
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Step 1 encompasses the generation of the suction surface of the section 
envelope from the leading edge to the first captured Mach line (point 2, fig 
9),. For all three types of sections the shape of this surface is defined by 
specifying the surface angle at a finite number of meridional locations. If 
conditions at the j th point are denoted by subscript j, then for the central 
sect i ons 






(V cos Pp (tan Pj 




,-iHr 


= F p, + 
s K fs . 

J 





(38) 


where the subscript f denotes the "no work" streamline. 

When F g = 1 , eq . ( 38 ) reduces to the definition used for the outer surface; 
where F = 0, it reduces to the definition used for the inner surface. The 

^ J 

parameter K in the above equation is defined as follows: 


ini' ## 



where the numbered subscripts denote points indicated on the reference blade 
sect i on (fig. 9 ) . 


With the surface angle defined at each discrete point j, the corresponding 
angular location of the surface 0 is calculated incrementally as 


e. 

j 


■iHf 

0 

J - ' 


■SS-rt 

'tan p. + tan P 
J 


r. + r. . 
J J-1 




m 


j -1 


) 


(40) 


The angular location of the first point (at the leading edge) is specified as zero. 


For the outboard sections, the termination point of the above surface is 
determined by its intersection with the first captured Mach line, the first 
captured Mach line being constructed from the free streamline angle, the Mach 
number distribution, and blade spacing. For the inboard sections, the termi- 
nation point is determined by the intersection of the surface with the passage 
inlet normal. For the central sections, the termination point is determined 
by its intersection with a pseudo-Mach- 1 i ne defined in terms of the surface 
angle and the free streamline Mach number distribution. The term " pseudo-Mach- 
line" is employed in recognition of the fact that the coordinates of this line 
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are calculated for a meridional distribution of Mach number that has not been 
corrected for adjustments to the suction-surface free streamline beta distri- 
bution imposed by eq. (38). This pseudo-Mach- 1 ine reduces to a true Mach line 
for the outboard sections where the surface and free streamline angles are 
identical. Similarly, it reduces to a passage normal close to the inboard sec- 
tions as the free streamline Mach number approaches unity. 


Step 2 encompasses the generation of the section envelope suction surface 
from point(2), fig. 9, to the trailing edge. The surface is defined in two sub- 
steps, the first referred to as the preliminary surface definition, and the 
second as the final surface definition. The preliminary surface definition is 
used in determining the incidence point of the leading edge shock. Once this 
point is located, appropriate modification of the preliminary surface down- 
stream of this location yields the final surface definition. 


The preliminary angle distribution for the critical sections is generated 
by the finite difference equation 

■JH'r -*Wr III 


where 


3. - P. , + K (m. - m. - ) 

J J-1 J J-1 


K" 1 = F K" + (1-F ) K 1 (1 + ^ -^-) 


(41) 

(42) 






K " . * f s 6 i : 

TD’ td 


m 


and K is given by eq . (39). 


To demonstrate that the structure of eqs. (41) through (43) applies also 
to outboard and inboard blade section design, eqs. (41) and (42) can be combined 


to display the finite-difference approximation dp**/dm in the region between 
points (l) and 


4 \ Q /V ** 

i - Pi i ii m. 

^ = F K + (1-F ) K 1 (1 + X -J- 
dm m.-m., s x s m/^\ 

j j - 1 © 

For outboard sections (F = 1), eq. (44) becomes 


) 


(44) 


pr - p: 




r\ + 6, - PI' 


m. - m. - 
j j- 1 


m © _ m © 

which is equivalent to eq. (18). (Note that the term (P^ - Pi*) of eq. (18) 
equals - P^.) For inboard sections (F = 0) eq. (44) reduces to 


(45), 


- P . - P/T\ ~ P/TN m • 

-1 = £) TO ,, + X _L) 

j j_1 m © (1 i) © 


(46) 
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which is identical to eq. (27)* 


The location of point®, the shock intercept point, is determined from 
the interception of the leading-edge shock and the preliminary suction surface. 
For cases in which the leading-edge wedge is too large to support an attached 
shock or where the designer wishes to control the location of the "kink" inde- 
pendently of the shock, the intersection location is specified manually. The 
expression used to determine the suction surface angle distribution downstream 
of the intersection is 


■iHS- 


3 

© 



F_ 6, 


(47) 


where is a specified parameter and 6^ is the leading-edge shock strength 

(see eq . (17))- For F g = 1 , the suction surface technique described above 

generates a surface which is consistent with the outer surface technique, and 
for F =0, generates a surface consistent with the inner surface definition. 

t- S 


The technique used in step 3 is identical to that used in the outer sec- 
tion procedure for defining the change in G"*** between points (5) and © . 


Step 4 encompasses the generation of the pressure surface from point® to 
point®- For the outboard sections, the surface angles between poi nts (6) and Q) 
are set equal to the corresponding tra i 1 ing-edge free streamline angles com- 
puted from the trailing edge RVg, and point (j) is determined by the intersection 
of the adjacent blade suction surface t rai 1 i ng-edge shock with the free stream- 
line surface. For the inboard section, the surface angle is defined by a cubic 
curve fit and the location of point© is specified at input. For a point j 
between po ints ©and®, let the section surface angle as calculated by the out- 
board section design method be denoted as 3j r “ c , and let the surface angle as cal- 


VwVf i 

culated by the inboard section method be . 
is defined as 


Then the central surface angle 


•SHf- 

3. 

J 



, % ## : 
+ c - Vb 


(48) 


where F^ is a parameter governing the transition from inboard 
sections. For F = 1 , the surface corresponds to an outboard 
and for F^ = 0 , the surface corresponds to an inboard section 


to output 
section surface 
surface . 


Step 5 defines the generation of the section envelope pressure surface 
from point® to point©- For both inboard and outboard sections, the surface 
is generated by a cubic equation whose parameters are defined by the location 
and the surface angles at poi nts (§) and® Since the location of point(8)is 
already known from the above steps and point® can be easily located in terms 
of the location of point©, only the angles need to be specified to complete 
the definition of the middle section. The angle at point® is defined in the 
same manner for all three types of sections, i.e., 


tt-tt -iHt 

b ‘ & 



(49) 
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( 50 ) 


For point(3), the central section angle is defined as 


P** = g** 
© ® 


F 6. 
P k 


where and 6^ are now parameters specified by the designer. For = 0, the 

angles at points (7)and(8) are the same and the surface definition is identical to 
the inboard section specification. Similarly with Fp = 1 , and the kink angle 6^ 

set equal to the tra i 1 ing-edge shock strength, the central section reduces to 
the outboard section specification. 


Step 6 consists of the development of the blade section within the section 
envelope. It is performed in the same manner as described for an outboard 
section by subtracting the boundary layer displacement thickness from the sec- 
tion envelope. 


Summary of central section design sequence .--The central section design 
procedure is the most complex of the three design procedures discussed in this 
report. This is so because the midsection procedure must be capable of pro- 
ducing (with appropriate input parameters) either an outboard or an inboard 
section or any section in between. In fact, in the computer program that has 
been written to automate the section design procedure, only the midsection 
procedure was actually programmed. The input to this program for each type of 
section is identical in format, although many of the input parameters are not 
used for the outboard sections since the parameters can be calculated inter- 
nally. The major aspects of the transition from inner to outer sections are 



section. The parameter F & controls the transition of the suction surface, and 


the parameter 


F controls 
P 


the transition of the pressure surface. 


The section using the mid-section design technique requires two types of 
input: (l) input supplied by the axisymmetric program, and (2) input supplied 

by the designer. The input supplied by the axisymmetric program includes the 
inlet and discharge flow conditions, the streamtube geometry, and the axial 
length of the section. The input parameters supplied by the designer are fur- 
ther divided into two categories: primary and secondary parameters. The 

latter refers to quantities that are, in general, held constant during the 
design of the section. The former refers to quantities that are varied during 
the design process and by means of which the designer controls the section 
mechanical and aerodynamic properties. The primary and secondary input quanti- 
ties required for development of the mid-sections are tabulated as follows. 


Primary input 

\ 
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Descr ipt ion 

Suction surface curvature factor 

Suction surface incidence factor 

Section deviation angle with respect to envelope 
mean 1 i ne 

Section envelope leading-edge wedge angle 


42 



6 te 

F 

S 

f p 

7D 

Secondary input 


Section envelope t ra i 1 i ng-edge wedge angle 

Suction surface transition factor 
Pressure surface transition factor 
Meridional location of point 3 
Meridional location of point 7 


Q Profile loss factor for use in mixing calculation 

P 

r 

H Tra i 1 i ng-edge boundary layer shape factor 

Suction surface initial boundary layer growth angle 

Pressure surface initial boundary layer growth angle 
r^ ^ Leading-edge radius 

r_p ^ Tra i 1 i ng-edge radius 

The manner in which the primary quantities tend to affect the section 
properties is discussed in the following paragraphs. 

The leading- and t ra i 1 i ng-edge section envelope wedge angles together 
control the magnitude and location of the section maximum thickness. Increas- 
ing either wedge angle tends to increase the maximum thickness. For small 
suction and pressure surface angle discontinuities, these two angles also tend 
to control the location of the section maximum thickness. Increasing the 
leading-edge wedge angle tends to shift the maximum thickness forward, and 
increasing the trai 1 i ng-edge wedge angle tends to shift the maximum thickness 
aft. For larger suction and pressure surface “kinks' 1 , the location of the 
kinks tends to control the location of the maximum thickness (unless the kinks 
are placed far forward or far aft on the section). This is because the sec- 
tions tend to become single or double wedge shapes with the maximum thickness 
located at or between the wedge vertexes. 

The section envelope deviation factor controls the metal section deviation 
angle. The suction surface curvature and incidence factors and the location 
and magnitude of the surface angle discontinuities control the blade passage 
area distribution as well as the suction surface Mach number distribution. 
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Overall Rotor Design Procedure 

In the previous sections of this report, the techniques used to design 
the outboard, central, and inboard rotor sections were discussed in detail. In 
the present section, the manner in which these techniques were used to arrive at 
the final geometry is discussed. 

The axisymmetric calculations indicated that the five outermost rotor 
sections (streamlines 1 through 5, fig- 1) were fully supersonic, and hence 
they were designed using the outboard section design techniques. The primary 
variables associated with the design of these sections are the axial projec- 
tion of the sections, and the difference between the leading- and trai ling- 
edge shock strengths. The major section parameters monitored during the design 

process were the maximum thi ckness-to-chord ratio t / C and the start marqin. 

m 

The start margin as calculated by the methods of Appendix D varied from approxi- 
mately 5.5 percent at the tip to approximately 2.5 percent at section 5. The 
t /C values were set equal to values previously deemed acceptable from mechani- 
cal design considerations. Using the program described herein and observing 
the constraints resulted in a set of sections having gradually decreasing chord 
lengths and increasing shock strength difference radially inward from the tip. 
The increasing shock strength difference was necessary in order to keep the 
section passages sufficiently open. That is, the channel Mach number in the 
exit region was increased by means of additional turning towards the axis down- 
stream of the expansion kink (eq. 19))* Turning the channel towards the axis 
increases its flow area. 

As a consequence of designing the five outboard sections first, some 
guidance in the selection of the variables required as input for the design of 
central and inboard sections could be gained. Specifically, this information 
included (l) the envelope leading- and trai 1 i ng-edge wedge angles, (2) the 
locations of the suction and pressure surface kinks, (3) the magnitudes of the 
respective kinks, i.e., the values of Fs and Fp, (4) the magnitudes of the 
incidence factors, and (5) the magnitudes of the deviation angles. 

In order to gain some idea as to the probable behavior of the above vari- 
ables from section 6 on down, section 12 was developed next. This section has 
subsonic inlet and discharge Mach numbers and hence was designed as an inboard 
section (no kinks on the pressure or suction surface). In choosing the input 
variables for this section, the major parameters monitored were the passage area 

distribution, the t /C, and the inlet incidence angle. The location of the 

m 

hypothetical suction surface kink was specified just aft of the passage inlet 
normal. The location of the hypothetical pressure surface kink was determined 
by ext rapol at ion f rom the outboard sections. 

With the central input variables specified in the outboard region and at 
the hub, the approximate values of the variables at the other sections were 
determined by curve fitting. The variables were modified somewhat in the design 
of the remaining sections in order to achieve acceptable distributions of pas- 
sage area suction surface Mach number, start margin, and t^/C. 


kk 


After the final smoothing of the input, r, 0 , and z coordinates of each 
section were input into the rotor "stacking program". This program was specifi- 
cally written as a part of the present design effort to handle sections with 
surface angle discontinuities on the pressure and suction surfaces. The program 
maintains the integrity of the kinks during the stacking and twist operations, 
and fairs the kinks over a specified surface distance after development of the 
section onto planar surfaces. In the present design, all kinks were faired over 
a 0.2 in. surface distance. In addition to fairing the kinks, the program also 
fairs a leading and trailing edge ellipse to each section. In the present design,, 
a 2:1 ellipse was used for all sections as a compromise between achievable manu- 
facturing tolerances and optimum aerodynamics. 


DESIGN RESULTS AND DISCUSSION 

The rather elaborate design procedure used in this program for the rotor 
design was felt to be justified on the basis that a conclusive evaluation of a 
new compressor concept was desired. The stator design procedure was standard 
inasmuch as the stator flow was entirely conventional with respect to both turn- 
ing angle and Mach number. The fan stage which resulted from these design 
efforts is summarized and discussed below. 

Flow Path and Velocity Diagrams 

The flow path was determined considering influences of blade blockage, 
boundary layer al lowance, shock and profile loss coefficients, and work input 
distribution within the blade row. The resulting iterated flow path showing 
the locus of the rotor and stator leading and trailing edges and design stream- 
lines is presented in fig. 1. The 12 streamlines shown divide the flow into 
11 actual flow annuluses. 

The final velocity diagram values for streamline Nos. 1, 6, and 12 at 
design speed and airflow are presented in fig. 11. The most noteworthy factors 
are the supersonic exit Mach number relative to the rotor tip and the small 
amount of fluid flow turning required at streamline No. 6. The tip inlet Mach 
number of 1.647 is higher than present commercial practice, but not beyond pub- 
lished data for transonic and supersonic compressors (refs. 3 and 4). The maxi- 
mum absolute Mach number leaving the rotor is only 0.732 at an angle of 39*3 
deg from the axial direction (at the hub). Providing that the flow leaving the 
rotor is as uniform as the design concept predicts and the flow parameters have 
approximately the design values, a low-loss stator design appears routine. 

Rotor Design Summary 

The overall summary of the final rotor design is presented in tabular form 
in table 1. Because the axisymmetri cal flow calculations are carried out along 
stream surfaces, all of the design values are presented as functions of the 
streamline numbers rather than the dimensional radius. However, for orientation 
purposes, plots of streamline location at rotor blade inlet and outlet are pre- 
sented in fig. 12. In addition, 4 rotor sections defined along the conical sur- 
faces representing streamlines 1, 6, 9, and 12 are presented in figs. 10a 
through 1 Od . The presentation of the aerodynamic design results is completed 
by the following discussion and the accompanying figures. 
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Rotat ion 


St reaml i ne No. 

1 

6 

r 
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61.2 

51.1 
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Figure 11. --Rotor Velocity Diagram Values for Tip, Mid-Span 
and Hub Stream Surfaces of Design. 











TABLE I 


ROTOR AERODYNAMIC SUMMARY 

Pjl = 14 696 psia = 10 133 N/cm^ 

Tjl = 518.69 °R = 288 I6°K 
N//T = 12 781 rpm - 1338 4 rad/s 
W/T/6 = 147 91 Ibm/s = 6 709 kgn/s 
* 1^,= 1603 0 ft/sec = 488 594 m/s 
r T | = 14 370 in = 36 50 cm 


Blockage factor = 0 980 Inlet 
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3 
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0 

0 

0 

0 

0 

0 

0 

0 

0 
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3 
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1 1394 

1 1407 

1 1442 

1 1488 

1 1498 

1 1444 

1 1400 

1 1375 

1 1388 

U/U T 
1 1 

0 9701 

0 9388 

0 9068 

0 8747 

0 8392 

0 8010 

0 7613 

0 

7189 

0 6722 

0 6228 

1 0 5713 

' 0.5150 

ve/u T| 

0 7854 

0 7573 

0 7168 

0 6818 

0 6376 

0 5858 

0 5253 

0 

4660 

0 4117 

0 3512 

0 2789 

0 1884 

V^t, 

0 4735 

0 4317 

0 4086 

0 3893 

0 3774 

0 3699 

0 3656 

0 

3637 

0 3643 

0 3755 

0 3930 

0 3993 

0 

-11 5 

' -8 2 

-6 U 

-4 3 

-2.5 

-0 7 

0 8 

2 

7 

4 0 

7 6 

12.0 

17 7 

V/Uti 

0 5072 

0.4654 

0 4510 

0 4354 

0 4298 

0 4292 

0.4361 

0 

4429 

0 4479 

0 4641 

0 4910 

0 5170 

v’/u T( 

0 9170 

0.8715 

0. 8260 

0 7835 

0 7405 

0 6925 

0 6400 

0 

5908 

0 5496 

0 5140 

0 4822 

0 4417 

&' 
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8 
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Design flow regimes . --The Mach numbers relative to the rotor at inlet and 
exit are presented in fig. 13. At rotor inlet the outer half of the flow may be 
considered supersonic and the inner half transonic. Flow exiting from the rotor 
is supersonic from the tip to streamline no. 5, and subsonic inboard from that 
locat i on . 

Relative flow angles . --The distribution of the airflow angles relative to 
the rotor blade at the rotor inlet and outlet are presented in fig. 14. The 
relatively high values of the inlet flow angles, varying from about 50 to 65 deg, 
> reflect the high rotational speed of this rotor. The distribution of the rotor 
blade outlet flow angles is somewhat unusual, showing a maximum value at 
streamline no. 4 with the minimum required turning angle occurring in this zone. 
As mentioned previously, this variation of rotor outlet flow angle and of the 
required turning angle, g j minus g ^ > results from the wall curvature effects on 

the radial variation of rotor inlet and outlet meridional velocity. 

Section solidity and th i ckness .--The rotor blade solidity and section thick- 
ness to chord ratio (t^/c), figs. i~5 and 16, respectively, result from compro- 
mises between the aerodynamic and mechanical design considerations. The rotor 
tip solidity of about 1.62 was chosen solely from aerodynamic design considera- 
tions. The number of rotor blades, their taper, and average chord were funda- 
mentally fixed by the contractually specified aspect ratio and structurally 
allowable chord and cross-sectional area taper. Initially a linear variation 
with radius of maximum thickness in percent of chord length was assumed. The 
final distribution shown in fig. 16 does not vary greatly from this selection. 
Fig. 17 gives the chordwise location of the maximum thickness as a function of 
streamline number. 


Total pressure ratio and loss coefficients .--The objective radial distri- 
bution of stage pressure ratio for this transonic fan is constant. The varia- 
tion in rotor outlet total pressure ratio (fig. 1 8) therefore reflects the 
expected higher total pressure loss in the stator end wall regions. The rotor 
outlet total pressure ratio is highest at the hub because the Mach number, turn- 
ing requirements, and hence loss for the stator vane sections are greatest at 
the hub. The rotor diffusion factor calculated using the axisymmetric program 
is presented in fig. 19* The profile component of the rotor total loss coeffi- 
cient, u) , is considered to be a function of the diffusion factor and radius. 


The variation of this profile loss coefficient, u)^, is presented in fig. 20 and 

shows the trend with streamline number as would be expected from diffusion factor 
variations. The component of rotor total loss coefficient at design operation 
due to supersonic flow and shock waves is expected to be different in this rotor 
from any previous design because of the avoidance of strong tip shock. 


There is no direct precedent upon which to estimate the total pressure loss 
coefficients for this design wherein the total is composed of shock and profile 
losses. In the al 1 -supersoni c outboard region only weak oblique shocks are 
expected so the shock losses are anticipated to be very low. This should also 
be true near the hub contour where the relative flow has either subsonic or low 
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transonic values. In the intermediate zone between streamlines nos. 6 and 10, 
bow waves are anticipated. The strong shocks are expected to occur in the 
region around streamline no. 8 (see fig. 20). The additional total pressure 
loss due to rapid boundary layer growth and probable separation is expected to 
be maximum inthis region. For this reason, the additional component of loss 
coefficient, u) s , attributable to shock losses, has been given a maximum value 

of twice the normal shock loss in this zone (fig. 20). From this maximum, addi- 
tive values of ou s have been tapered off smoothly to low values in the subsonic 

and in the al 1 -supersonic regions. At both walls, the loss coefficient has been 
increased to reflect wall boundary layer and clearance effects, but the increase 
is anticipated to be small due to the design condition of low diffusion across 
this rotor and particularly across the hub and tip elements where high calculated 
meridional velocity ratios exist at the walls. These considerations result in 
the fig. 20 variation of shock loss coefficient with radius. The resulting com- 
bined total -pressure- loss coefficient is also indicated in fig. 20. 

Although a localized loss associated with the midspan damper will occur, 
the decision was taken to distribute a l-percent loss in relative total pressure 
uniformly across the rotor exit. This distributed loss is included in fig. 20. 
The reason for not attempting to model the loss locally is that for an essen- 
tially constant total pressure ratio the computer program would tend to give a 
locally increased camber yielding small local blade irregularities which could 
compromise mechanical integrity in the high stress region of the midspan damper. 
Due to the high sweep angle of the damper (71 deg) and the low normal Mach num- 
ber (0.67), there is reason to expect minimum flow disturbance by the damper even 
though the relative Mach number of 1.43 is higher than used in current practice. 

Element temperature rise and efficiency . — The normalized rotor total 
temperature rise presented in fig. 21 reflects the work input required in order 
to produce the desired stage total pressure ratio while considering the combined 
total loss coefficient anticipated to occur. The calculated element efficiency 
of the rotor for each streamline is shown in fig. 22, presenting a near mirror 
image of the fig. 21 temperature rise distribution. In practice, it might be 
expected that lower efficiencies than the 85 percent calculated at the tip may 
occur in that region and that higher values than the 82 percent minimum calcu- 
lated between the streamlines nos. 7 and 8 may also be realized. Hcwever, the 
calculated mass-averaged rotor efficiency of approximately 86.7 percent (includ- 
ing the distributed one percent total pressure loss allowance for the midspan 
damper) appears to be a valid expectation. 

Incidence, deviation, and blockage . — The incidence was primarily influenced 
by the section incidence correction as defined in Appendix D. The resulting 
incidence angles with respect to the suction surface and mean camber line are 
shown in fig. 23. 

In the all-supersonic outboard region of the rotor, the deviation angles, 
fig. 2k, were determined by aerodynamic analysis during the calculation of the 
airfoil sections. For the inboard zone of the rotor having subsonic exit Mach 
number, the deviation angle was , cal cul ated using Carter's rule with an adjusted 
additive factor shown in fig. 2k. 
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Figure 23. --Rotor Incidence Angles 



and deviation angle, 



Figure 2k . --Rotor Deviation and Carter's Rule Adjustment 






An inlet total blockage factor of 0.98 and an exit total blockage factor of 
0.955 were assumed for the design pass and held constant span-wise. Realistic 
blockages based on a preliminary design were used through the blade. Typical 
rotor blockages for streamlines 1, 6, and 12 are listed below. 


Streaml ine 

Percent Axial Chord 

0 

12.5 

25.0 

37.5 

50 

62.5 

75.0 

87.5 

100 

I 

0.980 

0.945 

0.928 

0.919 

0.910 

0.903 

0.894 

0.898 

0.955 

6 

0.980 

0.960 

0.916 

0.896 

0.870 

0.834 

0.831 

0.890 

0.955 

12 

0.980 

0.870 

0.820 

0.794 

O.78I 

0.784 

0.810 

0.860 

0.955 


Rotor start margin . --The rotor starting margin for streamlines 1 through 8, 
which include the all-supersonic outboard sections and sections similar in shape 
to these outboard sections, was calculated by the procedures described in Appen- 
dix D. Start margin is defined as the percentage by which the channel flow area 
exceeds the minimum area needed to pass the normal shock. These characteristics 
are presented in fig. D-1 , Appendix D. As shown therein, the start margin varies 
from approximately 5 percent for section 1 to about 2.5 percent for section 8. 

Additional rotor calculations . — In the process of the current design, three 
separate versions of the axisymmetric design program were used. 

(1) The streamfi lament program with calculat ional stations at the inlet 
and outlet of the blade rows (total of 14 computational stations) 

(2) The streamfi lament program with calculat ional stations at seven axial 
positions inside the blade row (total of 28 computational stations) 

(3) A new axisymmetric program using a finite difference solution to the 
radial component of the equations of motion with blade interior grid 
points 

All aerodynamic results presented prior to this point were obtained from 
program (3) above. The calculated flow conditions including relative angle 
agreed well for all cal cul at ional procedures except for relative angle at the 
rotor trai 1 ing-edge station. The results are indicated in fig. 25, showing a 
significant difference, particularly at the inner radii. The finite difference 
results were selected for the final design because it was felt that the particu- 
lar calculation appeared to handle the blade blockage just upstream of the rotor 
trailing edge somewhat better than the stream filament program. 

Stator Design Summary 

The stator design is based on inlet and outlet velocity diagrams compatible 
with the rotor final aerodynamic design and zero stator outlet swirl. The axi- 
symmetric computer pass on which this design was based utilized calculating 
stations at the blade row leading and trailing edges. A blockage factor of 
0.955 which was compatible with the rotor tra i 1 i ng-edge value was used at the 
inlet to the stator. On the basis of available empirical data this blockage 
was assumed to decrease to 0.940 at the stator exit plane. A tabulated summary 
of aerodynamic parameters at the leading- and trail ing-edge stations is pre- 
sented in table 2 with some of the more pertinent parameters plotted in fiqs. 

26 through 28. 
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TABLE 2 


STATOR AERODYNAMIC SUMMARY 

p Tl = 14.696 psta = 10,133 N/ cm 2 

' TTI = 518.69 °R = 288. I6°K 

N/ AIT = 12781 rpm = 1338.4 rad/s 
w/e 7 fi= 147.91 Ibm/s = 6.709 kgm/s 
U T | = 1603.0 ft/s = *<88 59*4 m/s 
rjj ss 14.370 in. = 0.3650 m 


Blockage factor = 0.955 Inlet 


Stream! i ne 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

'2 /r T, 

0.9605 

0.9313 

0. 9008 

0 8690 

0. 8354 

0 7999 

0.7623 

0 7221 

0.6797 

0 635) 

0.5879 

0 5377 

Flow 

0 

9 09 

18.18 

27 27 

36 36 

45.45 

54 55 


72 73 

81.82 

91 91 

100.0 

wm 

1.512 

1.510 

1.509 

1 509 

1.509 

1 510 

1 .51 1 

1 512 

1.513 

1 514 

1 516 

1 527 

V t t, 

1.1474 

1 1421 

1 1399 

1.1396 

1.1406 

1. 1441 

1.1488 

1 1498 

1 1443 

1 1397 

1 1374 

1 1388 

ve/un 

0. 1878 

0.1871 

0. I 903 

0 1965 

0 2065 

0.2208 

0.2389 

0.2539 

0 2601 

0.2695 

0.2863 

0 3163 

vm/ut. 

0.4423 

0.4292 

0.4242 

0 4180 

0.4117 

0.4080 

0.4042 

0.4049 

0.4124 

0.4230 ; 

0.4410 


V/UTI 

0.4810 

0.4685 

0.4654 

0.4616 

0.4604 

0.4635 

0 4697 

0.4779 


0.5016 

0 5259 

0. 5689 

3 

23.0 

23.5 

24.1 

25.2 

26.6 

28.4 

30.6 

32. 1 


32.5 

32.4 

33.7 

M 

0.669 

0.653 


0 643 

0.641 

0 645 

0.653 

0.665 

0.681 

0 704 

, 0.743 

0 81 1 

‘MCL 

-2,90 

-2,62 

Eg| 

-1.87 

m 

-0.62 

| 0.25 


1.05 

1.28 

1.70 

2.30 

'ss 

-9,05 

-8.98 

-8.88 

-8.88 

-7.60 


-5.80 



-4.89 

-5.20 

-4.10 


Overal 1 parameters 


r 

0.9605 

0 9308 

0.9001 

0 8681 

0.8345 

0 7994 

0 7623 

0 7230 

0 6815 

0 6377 

0.5912 

0.5415 

AT/T, 

T l 

0.1473 

0 1417 

0.1398 

0. 1 396 

0 1407 

0.1442 

0 1488 

0. 1498 

0 1446 

0 1400 

0.1377 

0.1388 

U) 

0.030 

0.025 

0.024 

0 025 

0.026 

0.027 

0.0285 

0 030 

0.0315 

0 033 

0 035 

0 050 

D 

0.344 

0.321 

0.314 

0.309 

0.308 

0.316 

0 330 

0.342 

0.348 

0.358 

0 379 

0 424 

^ad 

0.834 

0.866 

0.878 

0.879 

0.873 

0 852 

0.826 

L 

0.820 

0 851 

0 879 

0 894 

0 885 

a 

1.010 

' 1.043 

1.080 

1. 1 19 

1.165 

1.220 

1.283 

1 359 

1.445 

1.542 

1 .643 

1.825 


Blockage factor = 0.940 Exit 


r 3 /r T, 

0.9605 

0. 9303 

0.8993 

0 8672 

0. 8337 

0.7988 

0 7623 

0 7238 

0.6832 

0.6403 

0 5945 

0. 5452 

m 

1.500 

1.500 

1 500 

1.500 

1.500 

1 500 

1 500 

1.500 

1 .500 

1 500 

1.500 

1 500 

asm 

1.1474 

1. 1421 

1 1400 

1 1396 

1.1407 

1.1441 

1.1489 

1 1498 

1.1444 

1.1398 

1 . 1374 

1 1388 

ve/Uj, 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Vm/I->t i 

0.4080 

0.4074 

0 4074 

0.4074 

0 0474 


0.4080 

0 4080 

0 4080 

0 4092 

0 4124 

0 4148 

v/Uti 

0.4080 

0.4074 

0.4074 

0.4074 

0.0474 

0.0474 

0 4080 

0 4080 

0 4080 

0. 4092 

0 4124 

0 4148 

e 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

M 

0.561 

0.562 

0.562 

0. 562 i 

0 562 

0.561 

0.561 

0 561 

0.562 

0.565 

0 570 

0 573 

6° 

7.50 

7.50 

7.60 

7.85 

8.39 

9.00 

9.95 

10.40 

9.86 

9.50 

9.19 

8.95 
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Figure 27* — Stator Inlet Absolute Mach Number. 
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Vane se lect ion .--Both a special uncambered leading edge section and a 
double-circular-arc (DCA) section were considered for this design. The inlet 
hub Mach number of only 0.811 (maximum)and corresponding diffusion factor of 
less than 0.43 (sol id i ty of 1.82) seemed near ideal conditions for a DCA stator. 
However, a general survey was made of the published literature on high and 
moderately high Mach number stators, and recent AiResearch experience on 
stators at inlet Mach numbers up to 0.94 was examined. The study showed strong 
evidence that DCA sections would perform at a level equal to or better than any 
other section type at the design Mach number and turning level (see table 2 and 
fig. 2 % To further verify the suitability of DCA stator sections, a suction 
surface velocity option in the axisymmetric computer program incorporating 7 
intra- blade- row stations (for both the rotor and stator) was used to approxi- 
mate the stator pressure and suction surface velocity and Mach number levels. 
The results are shown in figs. 29 and 30 for the hub and tip streamlines. 

A linear variation (with axial chord) of turning and loss was assumed and the 
calculations based on moment of momentum change between the designated stations. 
Realistic blockages based on a preliminary design were used through the vane. 
Typical stator blockages used are tabulated below for streamlines 1, 6, and 12. 
For the above assumptions, the resulting peak suction surface velocity for the 
critical streamline (hub) is less than sonic at the design point. It also can 
be seen from fig. 29 that the ratio of peak suction surface velocity to trail- 
ing edge velocity is less than 1.7. 


Percent axial chord 

Streaml ine 

0 

12.5 

25.0 

37.5 

50.0 

62.5 

75.0 

87-5 

100 

1 

0.955 

0.9135 

0.8948 

0.8833 

0.8804 

0.8847 

0.8947 

-d" 

0 

c\ 

• 

0 

0.9400 

6 

0.955 

0.9054 

0.8835 

0.8699 

0.8670 

CM 

00 

• 

0 

0.8848 

0.9041 

0.9400 

12 

0.955 

0.8829 

0.8513 

0.8323 

0.8288 

0.8376 

0.8569 

0.8862 

0.9400 


Losses . The assumed spanwise loss distribution for the stator shown in 
fig. 31 is compatible with that of ref. 5 for similar type blading, inlet con- 
ditions, and loading. 


Incidence and deviat ion . --The stator vane parameters of blade angle and 
camber angle corresponding to a given set of velocity triangles were defined 
utilizing the procedures given in Chapter VII of ref. 6. Some of the more perti- 
nent parameters are plotted as figs. 32 through 35. The design mean-cambei — line 
incidence is shown in fig. 35. The deviation based on Carter's Rule plus an 
additive factor is also indicated by fig. 36 with the additive factor shown 
separately in fig. 37- Small adjustments to original calculated blading param- 
eters were made to obtain smooth leading and trailing edges. All the above 
tables and figures reflect the final design. 

Choking inci depce . --The designed stator was evaluated for the incidence 
angle corresponding to sonic flow at the throat. This incidence angle, denoted 
the ideal choking incidence angle, imposes the minimum permissible design angle. 
This choking incidence angle is shown in fig. 38. The procedure for determin- 
ing the ideal choking incidence angle is thoroughly explained in ref. 7. 
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Figure 30. — Stator Surface Mach Number Distributions. 
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Figure 32.— Stator Solidity.. 







Streamline number 


s- 7<t060 


Figure 34.— Stator Stagger Angle. 
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Figure 37. — Stator Deviation and Carter's Rule Adjustment. 
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MECHANICAL DESIGN 


A relatively elaborate mechanical design was carried out on the 
subject transonic fan rotor because of the combination of relatively high 
tip speed and high aspect ratio and the aerodynamic necessity for a thin 
b 1 ade . 

This section presents the details of that analysis. Design criteria, 
applicable steady-state and vibratory stresses, as well as margins of 
safety are summarized herein. The study encompassed the rotor blades, fan 
disk, disk/blade attachments, stator vanes and rotor critical speeds. 

DESIGN CONCEPTS AND OBJECTIVES 

The design criteria cover the allowable stress levels, anticipated 
operating conditions, and material properties used to size the critical 
components . 


Material Properties 

The materials selected for the critical components are as follows: 


Component 



Materi al 

(Speci f icat i on) 

Fan blade 



Ti -6A1 

-1*V (AMS 1*928) 

Fan disk 



T i - 6 A 1 

-1*V (AMS 1*928) 

Stator vane 



17-** PH (AMS 561*3) 

The material 

strengths that form 

the bases for the al 

lowable stresses 

from Ml L-HDBK 

-5B 

(ref. 8) and are 

tabulated below: 




Fan blade 

Fan disk 

Stator vane 



(AMS 1*928) 

(AMS 1*928) 

(AMS 561*3) 

Temperature , 

°F 

250 

200 

300 


(°K) 

(395) 

(367) 

(1*22) 

Ft jj , ksi 


113 

117 

135 

(N/cm 2 ) 


(779 x 10 2 ) 

(806 x 10 2 ) 

(930 x 10 2 ) 

Ft v , ksi _ 
(N/cm ) 


102 

107 , 

115 7 


(702 x 10 ) 

(737 x 10*) 

(792 x 10 z ) 


80 


Allowable Design Parameter Values 

The following allowable design stresses refer to minimum material proper- 
ties based on MIL-HDBK-5B definition (ref. 8). 

Rotor blade . --The blade design limits were set on the basis of the follow 
ing critical safety margin requirements: 

(1) The maximum steady-state airfoil stress allowed is 82 ks i 

2 2 

(566 x 10 N/cm ) which is 80 percent of the minimum 0.2-percent 
yield stress. 

(2) The maximum midspan damper stress allowed is 102 ks i (703 x 

2 2 

10 N/cm ) which is 100 percent of the minimum 0.2-percent yield 
st ress . 

(3) Airfoil vibratory stresses shall not exceed 50 percent of the 
smooth specimen endurance strength based on a combined stress 

Goodman diagram at 10^ cycles. 

Disk . --The disk was sized based on the following allowable limits: 

2 2 

(1) The maximum allowable bore stress is 91 ks i (626 x 10 N/cm ) which 
is 85 percent of the minimum 0.2-percent yield stress (100 percent 
of the minimum 0.2-percent yield stress locally where bending exists 

2 2 

(2) The maximum web stress is 85.5 ksi (590 x 10 N/cm ) which is 
80 percent of the minimum 0.2-percent yield stress (100 percent 
of the minimum 0.2-percent yield stress locally where bending 
exi sts) . 

(3) Disk burst speed is 125 percent of the design mechanical speed 

at the maximum disk metal temperature, based on 90 percent of the 

2 2 

minimum ultimate tensile stress of 105 ksi (723 x 10 N/cm ). 

Blade/disk attachment . --The maximum allowable attachment stresses listed 
in table 3 were based on the following limits: 

(1) Peak neck tension or tang bending stresses do not exceed 90 
percent of the minimum 0.2-percent yield stress of the material. 

(2) Peak combined fillet stresses do not exceed 125 percent of the 
minimum 0.2-percent yield stress. 

(3) Peak maximum fillet stresses do not exceed the low-cycle-fatigue 
allowable stress for 10 000 cycles from rest to 110 percent speed. 

(A) Peak bearing stresses do not exceed 125 percent of the minimum 
0.2-percent yield stress. 
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(5) Peak tang shear stresses do not exceed 58 percent of the neck 
tension allowable stress. 

Stator vane . --The principal stationary aerodynamic component is the stator 
vane row. The stator vane design is based on an allowable limit of airfoil 
maximum stress of less than 80 percent of the minimum 0.2-percent yield stress 
and a margin of 2 on smooth specimen endurance strength based on combined 
stresses. 


uperdLing tuna 




From the mechanical design viewpoint, the most pertinent operating speeds 
were the following: 

Aerodynamic design speed (100%) 12 800 rpm (1340.4 rad/s) 

Mechanical design speed (110%) 14 100 rpm (1476.6 rad/s) 


Maximum blade temperature 
Maximum disk temperature 


250°F (395°K) 
200°F (367°K) 


DESIGN RESULTS AND DISCUSSION 
Rotor Blade 

Geometric def in i t ion . --The geometric summary of the final rotor blade 
design is presented in tabular form in table 3* The rotor airfoil coordinates 
of the resulting 12 rotor blade sections are presented in table 4. The coor- 
dinates are L-plane blade sections defined by the intersection of the blade 
with planes tangent to the cylindrical surfaces. A typical L-section is pre- 
sented in fig. 39* 

Steady stress . --The blade nominal centrifugal stress versus radius is 
shown in fig. 40. This stress does not include any bending stress and, there- 
fore, is useful as a basis for optimizing the blade tilt. 

The airfoil combined stress distributions due to aerodynamic and centri- 
fugal loads including blade untwist are shown in fig. 41a. These stresses were 
calculated at the aerodynamic design speed (100 percent) with the blade tilted 
tangentially 0.0175 in. /in. (cm/cm) to optimize the stress distribution. The 
airfoil stresses caused by design aerodynamic loads only are shown in fig. 41b. 
The effect of blade tilt on the stress distribution is 5 percent per 0.010 
in. /in. (cm/cm) tangential tilt. The effect of dimensional tolerance 
is less than ±2 percent of the blade stresses. 
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table 3 


ROTOR GEOMETRIC SUMMARY 

[l-P lane Blade Sections, defined by the intersection of the blade with planes tangent to cylindrical surfaces 


Blade height =* 7>728 in. Number of blades = 40 

(19.63 cm) 

Hub/tip ratio = 0.462 Aspect ratio = 2.76 



14.372 

13.800 

13. 100 

12.399 

11.700 

11.700 

10.300 

9.600 

8. 899 

8.201 

7 500 

6.641 

l\uO l US _ ' 

cm 

_3_6-50 

35.05 _.. 

33.27 

31 .4<3 

26.72 

27.64 

26.16 

24 38 

22.60 

20.8S 

19.06 

16.87 

design 

63. 17 

62.21 

61.20 

60.06 

58.86 

57.64 

56.06 

54.20 

52.61 



51.88 

51.24 

50 21 

* 

^2 design 

52.74 

55.47 

56.93 

56.25 

54.38 

51 .63 

48. 14 

44. 14 

37.52 

28 07 

15. 13 

2.42 

design 

59.82 

59.76 

59.36 

58. 15 

56.45 

54.25 

51 .35 

48.07 

44. 10 

39. 50 

34.25 

27. 87 

ft * 

Pj static 

63.17 

62.43 

61.93 

61.13 

60.06 

58.84 

57. 1 1 

54.98 

53.02 

52.00 

51.26 

50.21 

P 2 static 

52.74 

55.69 

57.66 

67.32 

55.58 

52.83 

49.19 

44.92 

37.93 

28. 19 

15. 15 

47.79 

ft * 

static 

59.82 

59.98 

60.09 

59.22 

57.65 

55.45 

52.40 

48.85 

44.51 

39.62 

34.27 

27.87 

Pretwi st 

0.00 

0.22 

0.73 

1.07 

1.20 

1.20 

1.05 

0.78 

0.41 

0.12 

0.02 

0.00 

* 

A3 

10.43 

6.74 

4.27 

3.81 

4.48 

6.01 

7.92 

10.06 1 

15.09 

23.81 i 

36. 1 1 

47 79 


1.792 

1.740 

1.702 

1.675 

1 . 669 

1.687 

1.702 

1.726 

1 761 

1 791 

i — 825 1 

! 866 

rA/\ 1 d 1 UlUl U ™ " 

cm 

L.55 • 

l*.42 

4.32 

4.25 

4.24 

4.28 

4.32 

4.38 

h.47 

4.55 

4.64 

4rf4 I 

Chord ,n * 

3 714 

3.565 

3.357 

3.194 

3.065 

2.939 

2.795 

2.664 

2.526 

2 394 

2.285 

2. 184 

cm 

9.43 ' 

9.06 

8.53 

8.11 

7.79 

7.47 

7.10 

~T7n 

6.42 

6.08 

5.80 

5.5Q 

T m a*. /C 

0.0177 

0.0248 

0. 0306 

0. 0370 

0.0445 

0.0493 

0 0558 

0.0612 

0.0673 

0.0722 

0 0782 

0.0858 

v s .« <*> 

0.0015 

0.0016 

0.0019 

0.0022 

0.0025 

0.0029 

0. 0033 

0 0038 

0.0045 

0. 0052 

0.0059 

0 0065 

r T.E./C (a) 

0.0014 

0.0016 

0.0019 

0.0022 

0.0026 

0.0030 

0.0036 

0.0043 

0.0052 

0.0064 

0.0072 

0.0080 


( a ) Leading and trailing edge radii are defined to be the semi-minor axes of 2-to-l ellipses 

Oo 
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TABLE 4 


ROTOR AIRFOIL COORDINATES 
[Section no. 1; radius 14.372 in. (36.505 cm)] 


Coordinates, in. a Coordi nates, cm 


X 

Y 

X 

Y 

- . y 6 2 4 

-1.0271 

-2.4 4 4 5 

- 4 ♦ 13 2 6 

- • 9 c 0 5 

-1.6316 

-2.4 393 

-4.1442 

- . 9 5 5 9 

-1.0313 

-2.4 3 05 

- 4 . 14 3 4 

- . y 5 2 y 

-1.0277 

-2.4 20 3 

-4.1343 

-.9314 

-1.5890 

- 2 . 3 o 5 6 

o 

0 

O 

• 

1 

-.7472 

-1. 25 2 1 

-1. 89 76 

- 3 . 1 0 0 2 

- . s> 1 0 3 

- . 82 00 

- 1 . 2962 

- 2 . 0 o 2 9 

- . 2 6 9 7 

- . 3o 32 

- . 6 d 5 1 

-.9733 

-.0233 

. 05 3 7 

-. 0 59 1 

.15 6 3 

• 2 3 1 3 

. 4 9 0 1 

• 3 8 7 o 

1.2449 

. 4 9 71 

. 9 2 3 6 

1 . 2o2o 

2.3459 

. 5 7 y 5 

1.0505 

1.4716 

2 . 6 6 o 1 

. 6024 

1 . 08 5 7 

1. 3 300 

2 . 7 p 7 7 

.0231 

1.2112 

1 . 3 8 7 6 

2.6477 

• o 4 7 4 

1.1568 

1 . 6 44 3 

2 . 9 3 8 3 

. 6 69 2 

1.1926 

1 . 69 9 7 

3 . U 2 9 7 

.6B6c 

1.2220 

1 . 7 4 4 0 

3.1037 

. 7741 

1.3675 

1.9662 

3.4734 

• o 3 7 3 

1.4964 

2.1775 

3 . 8 008 

. 9 06 3 

1.5687 

2.3020 

5.9645 

. y C 9 1 

1.5763 

2.309 1 

4.0037 

. 907 9 

1.5794 

2 . 3 0 6 0 

4 , U 1 1 c 

. 9 0 4 3 

1 .5796 

2.2975 

4.0125 

. 9002 

1.5773 

2.2865 

4.0062 

• o 9 3 5 

1.5721 

2.2 74 5 

3 . 99 3 0 

• o 4 5 4 

1.5128 

2.1472 

3 . o4 26 

.7919 ' 

1.4484 

2.0114 

3.0786 

. 7 3 7 7 

1.3818 

1.8738 

3 . 5097 

.0629 

1.3132 

1 . 734 6 

3.335o 

. 6 6 U 4 

1 . 2845 

1 . o 7 7 4 

3.2627 

. 6337 

1.2489 

1. 6095 

3 . 1722 

. 0 0 7 o 

1.2127 

1 .5438 

3.0801 

.5627 

1 . 1 75o 

1. 4 801 

2 . 966 6 

.3584 

1.1385 

1.4182 

2.6916 

.4586 

.9a22 

1.1&53 

2.4947 

.1392 

. 455 1 

. 35 36 

1.1556 

-.4717 

- . 6420 

-1.19 8 1 

-1.6307 

-.9415 

-1.5800 

-2.3914 

-4.0132 

-.9618 

-1.6232 

-2.44 30 

-4 . 1226 

Coord i nates 

on L-plane blade sections 

defined by the 

intersection 


of the blade with planes tangent to cylindrical surfaces. 



TABLE 4.--Cont inued 



ROTOR AIRFOIL 

coordinates 


[Sect Jon no. 2; radius 13 

-800 in. ( 35-052 

cm) ] 

Coordinates, in. 

a Coord i nates 

, cm 

X 

Y 

X 

Y 

-.9567 

- 1 . 5 9 2 9 

-2.4 29 9 

4 . u 4 5 o 

- * 9 5 4 t> 

-1.5976 

- 2.4 24 6 

4 • u 5 7 o 

-.9507 

-1.5973 

-2.414b 

4.05? 1 

- ♦ > 4 0 4 

- 1 . 5 9 3 5 

- 2 . 4 u 5 7 

4.0474 

-.9208 

-1.5570 

-2.35)4 

3.9561 

- . 7 5 8 6 

- 1 . 2 o 0 1 

- 1 . 9 2 o 7 

3 . 2 0 0 6 

- • 5 4 4 9 

- • o 7 9 o 

- 1 . 3 o 3 9 

2 . 2 3 4 1 

-.5305 

-,499b 

-•ob94 - 

1 . 2 6 9 4 

-.1128 

-.1197 

- • 2 c 6 5 

- • 3 o 4 0 

.1115 

. 2o 1 4 

• 2 c 3 0 

. 06 3 9 

. 3 4 3 0 

. 6 397 

,o712 

1 . o £. 4 7 

. 4 07 8 

.7414 

1 . 0 3 5 o 

1 . 8 a 3 1 

.4316 

. 7 7 9 4 

1.0901 

1 . 9 7 9 7 

.4540 

.8179 

1 • 1 5 4 6 

2.07 7 3 

. 4 709 

. 8 56 7 

1.2114 

2.1759 

• 4 9 8 7 

. b 9 5 9 

1 *26oo 

2.2 7 5 4 

.5306 

. 955 7 

1.347b 

2.4274 

• 6 4 7 7 

1.1067 

1 • 0 4 5 0 

2.9634 

. 76 0 9 

1 . 364 0 

1. 9 3 25 

3 . 4 6 4 6 

• 8 2 b 3 

1.4777 

2 . 1 0 3 o 

3.7532 

.6310 

1 • 4 b 5 o 

2 . 1 1 o 7 

3.7740 

♦ 8 3 0 0 

1 • 4 o B 3 

2 . 10 9 5 

3 • 7 6 o 2 

. 82 7a 

1.4901 

2 . 1 U2o 

3.7849 

.o235 

1 • 4 b 8 7 

2.0916 

3 . 7 a 1 2 

.8190 

1 . 4 8 4 4 

2 . U8 o 1 

3 • 7 7 o 4 

.7895 

1.4462 

2.0053 

3 • o 7 3 4 

. 73 00 

1 . 36b4 

1 .o543 

3.4757 

. 6099 

1 . 288 7 

1.7016 

3 . 2 7 3 2 

. 609 3 

1.2072 

1 . 347 5 

3.0063 

. 5597 

1.1400 

1 . 4 2 1 5 

2 . b 9 5 4 

. 53 3 7 

1 . 1038 

1 . 3 5 5 5 

2 . a 0 3 5 

. 5084 

1.0671 

1.2914 

2.7103 

. 4839 

1.0299 

1.2290 

2,6 160 

. 4600 

. 99 24 

1 . 1 oB3 

2.5205 

• 4 3 6 5 

• 9 5 4 2 

1.1087 

2.4237 

.1769 

. 52 1 0 

. 4 4 9 2 

1 . 323 3 

-.2185 

- . 17 4 2 

- . 55 4 9 

-. 4 4 24 

-.7739 

-1.2170 

- 1 . 9 6 5 o 

3 . U9 2 7 

-.9547 

-1.5853 

-2.4249 

4.0266 

- . 9 5 o 6 

-1.5916 

-2.4297 

4.04 25 


a Coord i nates on l-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 



TABLE 4.--Cont inued 


ROTOR AIRFOIL COORDINATES 
[Section no. 3; radius 13*100 in. (33-274 cm)] 


Coordinates, in. Coord i nates, cm 


X 

Y 

X 

Y 

- .909 2 

-1.5035 

-2.5092 

-3.818a 

- • 9069 

-1.5086 

-2.3036 

-3.831b 

- ■ . ■ } ' L .21 

-1.5085 

- 2 . 2 9 2 1 

- 5 . 8 3 1 1 

- • o 9 7 o 

-1.5042 

-2.2o 0 4 

- 5 . o t 0 5 

- • o 7 o 9 

-1.4713 

-4.231? 

-3.750/ 

- . i 2 b 4 

- 1 . 207 4 

-1.o4b1 

- 5 . j 6 0 6 

-.6520 

- . c70 1 

- 1 . 3 5 13 

-2.4099 

- . 5516 

-.5350 

- . o5 7 3 

-1.5530 

- . 14 2!? 

- . 20 50 

- . 5 0 2 u 

-.51^7 

.05 3 / 

.1244 

• 1 5 o <* 

• 3160 

.25 17 

• 4 4 5 9 

. o 5 9 1 

1 . 1 ; 2 3 

.2911 

. 5 1 Bd 

• 7 5 6 U 

1.517c 

.54 10 

. 556 7 

.0154 

1.4139 

• 5 4 3 o 

. 59 5 0 

.872/ 

1.5113 

• 3 6 5 4 

. 6 3 3 8 

. 9 e. o 1 

1 .6090 

• 5 c 0 9 

. 67 30 

. 9 8 1 o 

1.7093 

. 40 7 0 

.7128 

1. 0 3 5 7 

1 . c 10 6 

. 4 9 2 5 

. 8 7 7 6 

1 • 2 9 0 a 

2 .2 29 0 

. 6 229 

1.1243 

1 . ?d 2 1 

2.8556 

. 7 o 3 2 

1.3824 

1 . 9 5 o 5 

> . 5 1 1 2 

. 76 9 1 

1.5952 

1 . 9 5 5 3 

3 • 5 4 3 o 

. 76 9 2 

1 . 39 7 5 

1.953d 

5 • 5 4 y 0 

• 7 6 o 6 

1. 4003 

1.9522 

3.5967 

. 7 o 5 5 

1.4022 

1.9445 

5 • 5 o 1 6 

• 7 6 0 o 

1.4005 

1 . 9 32 3 

5.5572 

.16b 9 

1 . 39 5 6 

1.9200 

3.9447 

. 7 23d 

1.3509 

1 . o 5 8 5 

3.4312 

. 059 1 

1.2603 

1.6740 

3 . 20 1 2 

. 59 3o 

1 . 1682 

1 . 5u77 

2.9 cl 2 

. 5274 

1 . 07 45 

1 .3396 

2.129 5 

. 4 7 91 

1.0057 

1.2168 

2 . bb kb 

. 4 5 3 7 

. 96 93 

1.1523 

2.4620 

. 4 2 8 8 

. 9 3 24 

1 . 0o9 1 

2.3683 

. 4 04 6 

. 695 1 

1 .0277 

2.2735 

. 3o 1 1 

. 8 5 7 3 

. 96 79 

2 . 1776 

. 3 5 8 2 

.8192 

. 9 09 7 

2 • 0 o 0 7 

. 2 85 5 

. 695 6 

. 725 0 

1. 766 7 

. 0 24 5 

. 2432 

. 06 22 

.6176 

-.4876 

- . 6825 

-1.2335 

- 1 . 7 3 3 5 

- . d 9 U 4 

- 1 . 4620 

-2.2616 

-3.7134 

- . 9 0 d 3 

-1. 4986 

-2.3071 

- 3 • d 0 6 3 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 



TABLE 4.--Cont i nued 


ROTOR AIRFOIL COORDINATES 
[Section no* 4; radius 12.400 in. (31*496 cm)] 


a Coord i nates , in. 

a Coord i nates, cm 

X 

Y 

X 

Y 

- . 688 9 

-1.4215 

-2 . 2578 

-3.6106 

• • 8 8 6 5 

-1.4270 

-2.2517 

-3.6244 

- • c d 1 y 

- 1 . 4 26 7 

-2.2399 

- 5 . o 2 3 o 

- . b 7 6 5 

-1.4223 

- 2 . 2 2 6 2 

.>•6125 

* • b b u 0 

-1.3920 

-2.1793 

- 5 . 8 5 5 6 

-.7130 

-1.1498 

- 1 . a 1 1 0 

2 . 9 0 P 

- . 5' 2 8 5 

- . a4 30 

-1.3418 

-2.1411 

-.3437 

- . 5 39 7 

- . o 7 2 9 

- 1 . 5 7 0 0 

- . 1 5 d 9 

-.2410 

—.4036 

-,o120 

. 0262 

. 05 2 2 

. 0664 

.15 2o 

.2 IIS 

. 3 388 

. 5 5 7 o 

. o 6 0 5 

. 24 92 

. 395 3 

. 6 329 

1 . 004 1 

. 2 7 3 3 

. 4 3 2 b 

. 69 4 0 

1 . U 9 o 6 

. 296 5 

. 4 70 3 

. 7 5 3 0 

1.1944 

• 3 1 8 8 

. 5085 

. d 0 9 8 

1.2914 

. 34 0 3 

. 54 72 

• 8 6 4 3 

1 . 3d9d 

.3611 

. 586 5 

. 9 1 7 1 

1.4897 

.4543 

. 7624 

1.1538 

1.9365 

. 59 3 2 

1.0208 

1 . 5066 

2.5927 

. 7 4 3 2 

1.2952 

1 . 88 76 

3.2898 

. 749 9 

1.3098 

1 . 90 46 

3.3268 

. 750 1 

1.3120 

1 . 9052 

3.3323 

. 74 9 4 

1.3154 

1.9034 

3.3408 

. 746 0 

1.3176 

1.8947 

5 . 54 6 5 

. 7 40 7 

1.5157 

1 . 8 d 1 4 

3.3417 

. 7 3 5 3 

1.3103 

1 .8o77 

5 . 3 2 o 1 

. 69 7 5 

1.2590 

1.7715 

3 . 1 978 

. 6209 

1.1548 

1.5770 

2.93 52 

. 543 2 

1.0486 

1 . 3 797 

2 • 6 6 3 4 

• 4 6 4 5 

. 9404 

1.1799 

2 . 3885 

. 3 8 5 0 

. 8 30 3 

. 9 7 79 

2 . 1090 

. 3b90 

. 8079 

. 9 3 73 

2.0521 

. 3437 

.7713 

. 87 29 

1.9591 

.3191 

. 7343 

.8104 

1.8651 

. 295 1 

. 6968 

. 74 95 

1 . 7 6 9 9 

.2718 

. 6589 

. 6 903 

1.6736 

.2104 

. 557 1 

. 5 344 

1.4149 

- . 0289 

.1522 

-.0733 

. 3866 

- . 5000 

-.6803 

-1.2699 

-1.7279 

-.87 1 1 

-1.3825 

-2.2125 

-3.5 1 16 

- . 8879 

-1.4160 

-2.2553 

-3.5965 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces.^ 


TABLE 4.--Cont inued 


ROTOR AIRFOIL COORDINATES 
[Section no. 5; radius 11. 700 in. (29-718 cm)] 


a Coordi nates, in. a Coord i nates, cm 


X 

Y 

X 

Y 

- . a 8 0 a 

- 1 . 3 4 3 1 

-2.2371 

-3.4115 

- . o 7 8 3 

- 1 . 3 4 8 9 

- 2 . 2 3 j 8 

-3.4262 

- • t) 7 / i 

- 1 . 3488 

-2.2179 

-3.4258 

- . 0 G 7 2 

-1. 3 4 4 0 

-2.2026 

-3.4 137 

- . 0 4 y 7 

-1.3167 

- 2 . 15 8 1 

-3.5445 

-.7057 

- 1 . Udo 1 

- 1 . 79 25 

1 

l\; 

• 

-O 

a 

V*4 

—0 

-.52 3 b 

- . 802 4 

-1.3299 

- 2 • O 5 6 1 

- • 5 * 1 6 

-.5214 

- • o 6 7 6 

-1.3242 

- . 15 9 6 

-. 245 1 

- . 405 3 

- ♦ o 2 2 5 

. 02 2 4 

. 0260 

. O 5 6 y 

.0659 

,12 1b 

.17 11 

. 50 8 6 

• 4 3 4 0 

• 2 0 y C 

.3014 

. 5 308 

. 7 656 

.29 39 

. 4 3 3 4 

. 74 65 

1.1007 

.3765 

. 5 fc 6 7 

. 956 1 

1.4394 

.4509 

.7014 

1 • 1 6 0 6 

1 . 7 u 1 o 

. b 2 5 6 

.8199 

1.3355 

2.0824 

.0714 

1.0712 

1.7052 

2.7207 

. 7 6 0 3 

1.2245 

1 . y 5 i i 

5 . 1 1 U 3 

. 7 o 3y 

1.235o 

1.9404 

3 . 1 3 8 3 

. 76 3 3 

1.2390 

1.9387 

5.1470 

. 7 5 5 6 

1.2415 

1.9293 

3.1534 

, 7 5 3b 

1 .2397 

1.9145 

3.1487 

, 7 4 7 o 

1.2340 

1 . 898 9 

5.1344 

. 6988 

1.1729 

1.774b 

2 . 9 7 9 ^ 

. 59 9 8 

1.0481 

1.5233 

2 . 66 2 2 

. 4 99 2 

. 9 200 

1 . 2678 

2 . 3 3 o 7 

. 3 9 7 2 

. 7 8 86 

1 . 0089 

2 . 0029 

. 3 385 

.7113 

. 8 5 9 7 

1 . 8 0 o 7 

. 2 5 68 

. 5 987 

. 652 3 

1. 5207 

.17 8 1 

. 4 8 4 1 

. 4 52 3 

1 .2296 

.10 19 

. 3678 

. 258 8 

. 9 34 1 

. 02 7 9 

. 24 9 9 

. 070 7 

• 6 3 4 o 

-.1024 

. C37 1 

-.2601 

. 09 41 

- . 52 75 

- . o 9 0 o 

-1.3398 

-1.7541 

—.8640 

- 1 . 3068 

-2. 1944 

-3.3192 

-.8796 

-1.3369 

-2.234^ 

-3.3955 

a Coordi nates 

on L-plane blade sections 

defined by the 

intersection 


of the blade with planes tangent to cylindrical surfaces. 



TABLE 4. --Con tinned 


ROTOR AIRFOIL COORDINATES 
[Section no* 6; rad i us 11.000 in. (27-940 cm)] 


a Coord i nates, in. 

a Coordi nates, cm 

X 

Y 

X 

Y 

• • o 7 4 6 

-1.251o 

-2.2213 

- 3 . 17 9 4 

-.8720 

-1.257b 

-2.2148 

- 3 • 1 9 4 8 

- • d o o 4 

- 1 . 2 5 7 7 

- 2 . d 0 0 7 

- 3 . 19 4 7 

• • b b 9 o 

- 1 . 2 3 2 o 

- 2 . 1 o 3 9 

- 3 . 1 a 2 u 

-.8433 

- 1 . 2 2 o 5 

-2.1419 

- 5 . 1 d U 3 

- . o 9 7 7 

-1.0104 

-1.7721 

-2.3664 

- . 5 1 5 4 

-.7440 

-1.5089 

- 1 • d a 9 d 

-•3331 

- • 4 o 4 0 

- • b 4 b 0 

-1. 22 9 3 

- . i 509 

-.2291 

. 38 3 1 

-.381/ 

. 0 3 1 3 

. 0220 

. 07 9 5 

. 055 7 

.1357 

• 16 4b 

• 3 4 4 o 

. 4 1 a 1 

.2155 

. ^ 74b 

. 3 4 7 5 

• D 9 8 0 

. 29 41 

. 3 o 5 9 

. 7470 

• 9 o 0 1 

.3714 

• 4 9 7 6 

. 9 4 3 4 

1 • 2 o 4 4 

• 4 4 7 4 

.6107 

1 . 1 3 b 2 

1.3512. 

. 5 4 5 d 

. 7o 1 C 

1 • 3 a b 3 

1 . 9 330 

» o 9 b 6 

. 9 9 a 4 

1 . 7 7 4 3 

2.3 33d 

. 7920 

1.145b 

2 . U 1 17 

2.90 9 7 

. 7 9 65 

1.1560 

2.0230 

2 . 9 4 1 2 

. 7 9 5 5 

1.1612 

2.02 J 6 

2 . 9 4 y 4 

. 79 20 

1 . 1 g42 

2.01 17 

2 . 9 56 * 

• 7 b 5 g 

1.1625 

1.9953 

2 . 9 528 

. 7 7 6 o 

1.136b 

1 . 9 7 7 5 

2 . 9 3 fa 2 

.7151 

1 . 0a62 

1 . d 1 b 4 

2 . 7 5 a 9 

. 5 b 7 1 

.9416 

1.4911 

2 . 3y 1 7 

. 4 57 b 

. 7924 

1 . 1 o 2 2 

2 . o 1 d 7 

. 32 7 0 

• 6 3 3 6 

• 6 3 0 o 

1 . 6 22 1 

. 2 54 2 

.5513 

. 04 57 

1 . 4 00 3 

.17 36 

. 4509 

.4415 

1.14 5 3 

. 0567 

. 5462 

. 24 55 

• d b 4 4 

. 022 3 

. 24 34 

. 0567 

. 6 1 a 2 

-. 04 96 

.1308 

-.1258 

. 3 4 7 5 

-.1119 

.0411 

- • 2 o 4 3 

.1043 

-.3714 

- . 3745 

-.9434 

- . 9 5 1 2 

- . 7 4 7 5 

-1.0188 

- 1. 8987 

- 2 • 5 d 7 d 

- . d7 1 1 

-1.2403 

-2.2126 

- 3 . 15 0 2 

-.674 3 

-1.2450 

-2.2207 

-3.1725 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 


TABLE 4.--Cont inued 


ROTOR AIRFOIL COORDINATES 
[Section no. 7; radius 10.300 in. (26.162 cm)] 

a Coordi nates, in. a Coordi nates , cm 


X 

Y 

X 

Y 

-.0767 

- 1 . 1427 

-2.22 bo 

- 2 • 9 0 2 5 

- • o 7 4 1 

-1.1490 

- 2.2 20 1 

-2.9 164 

— # o o o 0 

-1.1492 

— 2 . 2 u 4 d 

-2.9 loo 

* . b 6 0 b 

-1.1439 

-2.1060 

-2.903d 

-.64 52 

-1.1230 

- 2 . 14 6 7 

-2.6523 

- . 6 y 7 o 

-.9193 

-1.7719 

-2.3349 

- . 5 1 y 

- . 6 7 7 5 

-1.3077 

-1.7209 

-.3321 

-.4444 

-.6435 

- 1 . 1 2 o d 

- . 14 9 4 

-.2169 

- . 3 7 9 3 

- . 5 5 0 fa 

.6334 

. 008 2 

. 0648 

. U 2 U 9 

.14 6 7 

.14 6 1 

. 3 7 20 

. 3 7 6 1 

.2193 

. 2 3 8 3 

. 5 56c 

. 6 0 ~j 2 

.2914 

. 3286 

. 7 4 0 0 

. c 3 5 1 

. 3o29 

.4197 

. y 2 1 o 

1.0661 

• 4 3 3 6 

. 5 1 1 2 

1.1019 

1 • 2 y b 4 

.4826 

. 5 7 5 2 

1 . 2 2 6 2 

1.461 1 

• 6 4 6 0 

. 794 5 

1 .6407 

2 . U 1 7 9 

. o 1 b 9 

1.0360 

2 . 0 7 9 9 

2.6314 

. 6 302 

1 . 052b 

2.1066 

2 . 6 7 3 b 

. d 3 2 b 

1 . 0b 1 3 

2.1147 

2.6957 

. 6 3 06 

1 . U 6 6 6 

2.1095 

2. 7Qy 1 

. 6 24 7 

1 . 067 6 

2.0947 

2.7120 

• d 1 6 6 

1.0642 

2 . 0 7 4 b 

2 • 7 0 3 0 

. d047 

1 . 0 53 2 

2.0439 

2.6750 

• b 5 0 9 

. 9 04 9 

1 . b 3 3 3 

2 . 2 9 d 3 

. 4 880 

.7416 

1 . 2 5 9 5 

1 . d 8 4 0 

. 3 2 4 3 

.3713 

. 8 2 3 7 

1.4511 

. 1 o 0 3 

. 3 9 3 5 

. 407 2 

. 9 99 3 

. 0678 

. 2889 

.17 2 1 

. 7 3 3 8 

-.0079 

.1979 

- . 0200 

. 5023 

- • 0 a 0 6 

.1045 

-.2048 

. 26 5 5 

- . 1 6 0 o 

. 0092 

- . 36 29 

. 02 3 3 

-.2166 

- . 0d79 

-.5353 

-.2232 

- • 2 d 3 1 

-.1831 

-.7190 

-.4650 

-. 604 3 

- • 6 6 4 6 

- 1 . 3 3 53 

-1.7393 

-.8620 

-1.1120 

- 2 • 1 o 9 4 

-2.8243 

- . 6750 

-1.1347 

-2.2224 

-2.6820 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 



TABLE 4. — Cont inued 


ROTOR AIPFOIL COORDINATES 
[Section no. 8; radius 9.600 in. (24.384 cm)] 


a Coord i nates, in. 

a Coordi nates, cm 

X 

Y 

X 

Y 

- « o 9 6 7 

-1.0490 

-2.2775 

- 2 . 6 o 4 4 

-.696 7 

-1.0494 

-2.2775 

-2.66 5 4 

- . c 9 4 u 

-1.0554 

-2.2706 

- 2 . Oof! o 

-.691 1 

- 1 . U 5 6 4 

- 2 . ^ o 3 4 

- 2 • o a 3 ; 

-.6674 

-1.0559 

-2.2640 

- 8 . c o 1 9 

— • D 6 3 3 

- 1 • 0 5 3 o 

-2 . 24 3 6 

- 2 . 6 7 o 7 

- . o ? 9 1 

-1.0505 

- 2 . 2 3 3 0 

- 2 . 6 6 6 3 

-.6735 

-1.0440 

— 2.2 1 6 t> 

- 2 • o p 1 t> 

- . d 6 4 2 

-1.0321 

-2.1951 

- 2 • o 2 1 P 

- • 6 0 o 5 

- . 9 5 o 4 

-2.0466 

- 2 . 4 3 4 3 

- . M5U 

- • 0 4 4 9 

- 1 . 8 1 6 1 

- . 14 6 9 

- . 6635 

-, 7 3 4 8 

- 1 . 5 8 3 5 

- 1 . o 6 6 a 

-.5319 

- . 627 9 

- 1 . 3 6 1 U 

-1.5949 

-.4403 

-.5237 

-1.1163 

-1.3301 

-.3467 

-.4217 

-.6656 

-1.0709 

-. 25 7 1 

-.3214 

-.6530 

- . 8 1 o 3 

- • 16 5 5 

- . 2224 

-.4202 

-.5649 

-. 0 7 3 8 

- . 124 3 

- • 1 o 7 ? 

-.3156 

.0176 

- . 0266 

. 04 6 2 

-.067t 

.1094 

.0712 

. 27 79 

.18 0 7 

.14 4 0 

.1089 

. 36 7 2 

. 2 7 6 5 

. 1 7 79 

.1447 

. *•; 1 7 

. 3 6 7 4 

.2111 

.1605 

. 5 36 1 

• 4 5 6 P 

. 2 4 4 3 

.2166 

.6^04 

. 5 4 9 8 

. 2 77 3 

. 2525 

. 704 4 

.6412 

.3103 

. 2886 

. 766 1 

. 7 32 9 

. 3 4 31 

. 3248 

. d 7 1 5 

. o 2 4 9 

. 3 7 5 8 

.3611 

. 95 4 6 

.9172 

• 4 U d 4 

. 39 76 

1.0373 

1 . 009 6 

. 4 4 06 

• 4 3 4 2 

1.1196 

1.1028 

.4813 

. 4602 

1.2224 

1.2196 

. 5 6 9 3 

. 5825 

1.4459 

1.4795 

. 65 7 2 

• 6 6 7 5 

1 . 669 3 

1.7462 

. 74 5 1 

. 7952 

1.8924 

2.0196 

. 8432 

.9188 

2 . 1 4 1 8 

2.3337 

. 65 26 

. 9 309 

2.1661 

2.36 4 5 

• o 5 6 8 

. 9366 

2.1762 

2 . 3 7 8 d 

.0599 

. 94 55 

2 . 16 4 0 

2.4015 

. d599 

. 9465 

2 . 1 b4 1 

2.404 1 

.8596 

. 94 9 1 

2.1833 

2.4107 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 



TABLE 4.--Cont inued 


ROTOR AIRFOIL COORDINATES 

[Section no. 8--cont i nued ; radius 9*600 in. (24.384 cm) 1 


a Coord i nates 

, in. 



a Coord 1 nates. 

cm 

X 


Y 



X 


Y 

• o 5 8 Q 


. 9b 1 8 

2 


17 9 3 

2 

.4175 

• c s 5 5 4 


.9333 

2 


1726 

2 

• 4 2 1 4 

. o 5 1 d 


. 9 b 36 

2 


1 o 3 5 

2 

. 4 2 2 2 

. o 4 7 ? 


. 9 5 27 

2 


15 2 6 

2 

. 4 1 y d 

. 3 4 2 3 


. 9 5 0 o 

2 


14 0/ 

2 

• 4 1 4 4 

. a 3 o 1 


. 9 4 74 

2 


1 2 c 7 

2 

• 4 0 6 4 

• o 2 8 5 


. 9396 

:> 


10 4 4 

2 

• 3 oV 1 

.7 4 2 3 


.6716 

1 


8 c 5 3 

2 

.2 13/ 

• Ci 4 6 3 


. 79 33 

1 


6 4 17 


. u 1 4 y 

.5^04 


.7123 

1 


39 7 d 

1 

. o G y 1 

. 4 54 3 


. o 2 8 5 

1 


15 3c 

1 

.5/6) 

. 3 5 8 2 


.5418 



9 o 9 d 

1 

. 3 7 o 1 

. 2622 


. 4 5 22 



6 t P u 

1 

• 1 4 o 6 

• 1 o 6 1 


. 359 6 



4 2 19 


. y 1 3 8 

. u702 


. 26 39 



178 4 


• o 7 0 5 

- . U 2 2 7 


. 16 7 9 

- 


0 5 7 7 


. 4 265 

- . Cl 5 9 9 


.1383 

- 


152 0 


. 325 9 

-.0963 


. 08 d 1 

- 


24 4 6 


.*257 

-.1321 


. 04 73 

- 


3 35 4 . 


.1202 

- . 16 7 2 


. 0060 

- 


4 2 4 7 


. j 1 5 3 

-.2017 

- 

. U356 

- 


5 12 3 

- 

. 0906 

-.23b7 

- 

. 0 761 

- 


5 9 5 o 

- 

.19 8 2 

-.2691 

- 

• 12 0 8 

- 


6 b 3 3 

- 

• 3 0 6 o 

-.3619 

* 

.1640 

- 


7 6 o 8 

- 

• 4165 

-.3343 

- 

. 2076 

- 


d 4 9 0 

- 

. 52 7 2 

-.3772 

- 

• 2 6 6 4 

- 


9 5 7 9 

- 

• 6 76 6 

- . b 7 7 3 

- 

• 4 6 5 8 

- 1 


3 1 39 

- 1 

. 1 o 3 1 

- . o b 7 a 

- 

. 6742 

• 1 


o 7 09 

- 1 

.7125 

- . 7 9 9 1 

- 

. 6897 

*“ il 


0^96 

-2 

. 2597 

- . b a 2 b 

- 1 

. 020 1 

*“ 2 


2 4 14 

-2 

.5910 

- . b 9 1 d 

- 1 

• 0 346 

- 2 


265 0 

-2 

.6*76 

- . d 9 4 5 

- 1 

. 0 39 o 

-“\ 


2 72 0 

2 

.o409 

-.6962 

- 1 

. 0448 

- 2 


2 7 6 3 

-2 

. 6 538 




-2 


2 7 7 5 

-2 

• 6644 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 



TABLE- 4.--Cont inued 


ROTOR AIRFOIL COORDINATES 
[Section no. 9; radius 8.900 in. (22.606 cm)l 


a Coord i nates, in. 

Coordinates, cm 

X 

Y 

X 


Y 

- • y 16 4 

-.9565 

-2.3326 

-2 

. 4 29 4 

-.>177 

-.9605 

-2.3306 

- 2 

. 4 3 9 6 

-.y 124 

-.9644 

-2.3175 

-2 

. 44 9 6 

- . y 0 a 4 

-.9640 

- 2 . 3 07 2 

- 2 

• 44 8 4 

- . V 0 3 o 

- . 9 6 1 9 

-2.2955 

— ') 

. 4 4 3 2 

- . a y 9 u 

- . 9 5 o 4 

-2.2834 

-2 

. 4 3 4 4 

- . d y 2 i? 

-.9514 

- 2 • 2 6 6 6 

* 

• 4 1 o 6 

“•663b 

-.9413 

-2.2449 


. 3 9 0 a 

-.o253 

— . o 7 2 o 

- 2 . 0 9 o 3 

-2 

.2163 

-.733b 

- . 7 6 d 7 

-1.6633 

- 1 

• 9 8 .. 4 

- • 6 4 1 y 

-.6669 

- 1 . 630 5 

- 1 

. 69 90 

-.5503 

-.5729 

- 1 . 3 9 7 7 

- 1 

. 4 55 1 

- . 4 5 d 6 

- . 4 o 0 1 

-1.1649 

- 1 

.2193 

-.3670 

-.3901 

-.9322 

- 

• 9 9 0 b 

-.2754 

- . 3 024 

-.6996 

- 

. 7 6 6 1 

-. 1639 

-.2167 

- • 4 6 7 U 

- 

• 5 5 u 2 

- . U924 

-.1324 

- . 2346 

- 

. 5 5 6 2 

- . U009 

-.0492 

-.0023 

- 

.12 4 8 

. 09 05 

. 0 3 34 

. 22 9 6 


• 0 d 4 8 

.1361 

. 07 64 

. 35 0a 


• 1 y 4 0 

.1663 

.10 37 

. 4 2 7 5 


. 26 3 3 

.1986 

.13 10 

. 504 3 


. 3 3 2 6 

. 2267 

.1583 

• 5 o 0 9 


. 4 0 2 1 

. 258 9 

. 1 o57 

. 6 5 7 5 


.4716 

. 289 0 

.2131 

. 7 3 4 0 


.5412 

.3190 

. 2406 

.6103 


.6110 

. 349 0 

. 2 662 

• 66 64 


• 6 a 1 1 

. 3 769 

. 29 58 

. 96 2 3 


.7513 

. 4 08 7 

. 32 36 

1.0360 


.8216 

• 4 7 5 a 

. 3666 

1.2084 


. 9a 1 9 

. 5705 

. 4 7 77 

1.4490 

1 

.2133 

. 665 4 

.571 1 

1 . 69 00 

1 

. 4 506 

. 760 3 

. 667 0 

1.9312 

1 

. 69 4 1 

. 6666 

. 7 76 7 

2.201 1 

1 

. 9 7 27 

. 8 77 1 

. 7 6 7 5 

2.2278 

2 

. 0003 

. 6822 

, 79 3 5 

2.2407 

2 

.0154 

. 88 6 6 

. 603 1 

2.2520 

2 

. 0 397 

• 6 d 6 6 

. 805 4 

2.2525 

2 

. 04 5 7 

. 666 7 

. 607 1 

2.2522 

2 

. 0500 

. 8853 

.6103 

2.2466 

2 

• 0 5 a 1 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 


TABLE 4. — Continued 


ROTOR AIRFOIL COORDINATES 

[Section no. S-- cont inued ; radius 8.900 ' n * (22.606 cm) 1 


Coordinates, in. 

a Coordi 

mates, cm 

X 

Y 

X 

Y 

. d o 2 o 

.0124 

2.2417 

2 • 0 o 3 5 

. a 7 & 6 

. 6 1 33 

2.231c 

2 . 0 6 5 o 

. o 7 3 7 

. o 1 30 

2.2192 

2 . 0 c 4 9 

« o c c 2 

.6114 

2 . 2 0 6 3 

2 . U 6 0 o 

• o o 2 6 

• d 0 o 6 

2.1907 

2 . U 5 3 o 

» d ? 2 0 

. o 0 2 5 

2.1641 

2 . 0 5 b 2 

.76 7c 

. 7 4 7 1 

1. 9 2 4 1 

1 . 6 9 7 c 

.6626 

• C 6 1 b 

1.667c 

1 . 7 3 1 o 

.64 77 

.6126 

1 . 39 1 2 

1 • 5 5 c u 

. 4 4 29 

. 5 3 9 3 

1 . 1 4. 4 9 

1 . 509 0 

. 3 3 b 1 

. 4 6 1 b 

■ b 5 b 7 

1.1729 

.2334 

. 3 7 9 9 

. 592b 

. 9 c 5 0 

.12 6 0 

. 29 35 

.3271 

. 7 4 5 5 

. 0 2 4 3 

. 2025 

.0617 

.5143 

- . 0 o 0 0 

. 10 5 7 

-.2030 

. 2 7 0 9 

-. 096 3 

. 090b 

-.244c 

. 2 306 

-.1310 

.0573 

- . 3 3 2 7 

.14 5 4 

- . 1652 

. U*33 

-.4195 

. 0 590 

- . 1 9 b a 

-.0113 

- . 5050 

-.0286 

- . 2 3 1 9 

- . 04 63 

-.5891 

- . 1 174 

-. 2646 

- . U a 1 7 

- • o 7 2 0 

-.2075 

- . 2 9 6 7 

-.1176 

- . 7 6 5 6 

- . 2 987 

- . 3 264 

- . 1540 

- • o 3 4 U 

-.3910 

-» 3 5 9 u 

- . 19 07 

-.9132 

- • 4 d 4 4 

-. 590 3 

- . 2279 

- . 99 1 3 

-.6788 

-.4315 

- . 2 2 b 6 

-1.095c 

-.7077 

-. 5625 

- . 44 85 

-1.428c 

-1.1391 

- . 6 9 4 U 

- • 6 2 7 4 

- 1 . 76 2o 

-1.6936 

-.6260 

-.o139 

- 2 . 09 a 1 

- 2 . U6 7 1 

-.903b 

-.9277 

-2.295c 

- 2 . 3 5 6 3 

-.9124 

-.9405 

-2.3175 

-2.3884 

-.916o 

- . 94 59 

-2.325c 

-2.4026 

-.9177 

-.9515 

-2.330o 

-2.4 1 c 7 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces.. 



TABLE 4.--Cont i nued 



ROTOR 

AIRFOIL COORDINATES 


[ Sect 

ion no. 10; 

radius 8.200 in. (20.828 cm)] 

a Coord i nates, in. 

a Coordi nates, cm 

X 

Y 

X 

Y 

-.9364 

- . b 7 1 0 

- 2 . 3 7 o 3 

- 2 . 2 1 2 3 

- • y 3 5 7 

-.6753 

- 2 . 3 76 5 

-2.2231 

-.930J 

— . a 7 9 7 

-2.3622 

-2.23 43 

- • 9 C 5 o 

-.6793 

-2.3509 

- 2 . 2 3 3 3 

- . 9 2 U 5 

-.0772 

-2.3379 

- 2 • 2 7 9 

-.9132 

-.6733 

-2 . 3 24 5 

- 2 • ^ 1 b 6 

-.9060 

-.0662 

-2 . 306 2 

-2.2000 

- . b y y 7 

-. 8 569 

- 2 . 2 b 3 1 

-2.1763 

1 

• 

U 

4^ 

-.7920 

- 2.1367 

-2.01 Jo 

- . 7 5 0 

- . 69 5 0 

-1.9054 

- 1 . 7 b 5 3 

-.6392 

-.6029 

-1.0743 

-1.5312 

-.3683 

-.5152 

-1.4435 

-1.3065 

-.4773 

- • 4 3 1 6 

-1.2129 

- 1 . 0 9 6 2 

-. 3b6S 

- . 3 ? 1 7 

- . * 6 2 4 

-•o931 

- . 29 62 

-.2750 

- . 7 5 2 2 

- . o y d 5 

- . 203 6 

-.2013 

-.3222 

-.5112 

- . 11 5^ 

- . 1300 

- . 29 2 5 

-.3302 

-.0248 

-.0608 

-.0626 

- . 15 4 4 

• 0 b 5 3 

. 00o7 

.16 6 4 

. 0 1 7 Ci 

.12 24 

. 0465 

• 3 1 0 o 

.12 3 2 

.14 9b 

. 0686 

. 3 d 0 5 

.17 4 2 

.17 7 3 

. 0 d d 7 

. 4 502 

. 225 1 

. 204 7 

.10 6 7 

. 5200 

. 2 7 59 

. 2 322 

.1266 

• 5 6 9 b‘ 

. 3 2 6 b 

. 23 9 7 

.1486 

. 6 596 

. 3 7 7 3 

. 2o72 

. 1 6 8 o 

. 7 295 

. 4 2 61 

.3147 

. 1 o 6 6 

. 7 99 3 

. 4 7 69 

. 342 2 

. 208 6 

. 6690 

. 5297 

• 3 6 9 6 

. 2287 

• 9 3 8 6 

. 5607 

• 4 6 3 4 

. 296 1 

1.1769 

. 7 57 1 

. 5 b 5 3 

. 3 7 4 5 

1 . 4 364 

.9513 

. 666 2 

. 4 52 1 

1.6972 

1.1483 

.7714 

. 5 308 

1.9593 

1.3461 

• 667 6 

.6196 

2.2544 

1.5736 

. b 9 6 7 

. 6262 

2.2827 

1.5955 

. 9049 

. 6 335 

2. 2983 

1.6090 

.9114 

. 64 3 3 

2 . 3 1 49 

1.6338 

.9122 

. 64 7 5 

2.3169 

1.6447 

. y 1 22 

. 64 7 7 

2.3169 

1.6451 

.9112 

.6515 

2.3145 

1 . 654 8 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces.. 


TABLE 4. — Continued 


ROTOR AIRFOIL COORDINATES 

[Section no. 1 0--cont inued ; radius 8.200 in (20.828 cm)] 


Coordinates 

, in. 

Coordinates, cm 

X 

Y 

X 

Y 

• y 0 8 6 

• o 3 4 4 

2 • 3 0 7 d 

1.6621 

• y 0 4 4 

. 6 562 

2 . 2 9 7 2 

1 • 6 6 6 7 

• o y 9 u 

• O 3 b b 

2 • <l c 3 5 

1 . 6 6 c 2 

• b y 2 7 

• 0562 

2 . 2 o 7 4 

1 . 6666 

• o c 5 u 

• 6 5 4 3 

2 . 2500 

1 . 06 20 

. b 7 4 7 

. 6502 

2 . 2 2 1 6 

1 . 6 0 1 O 

. 7 7 4 1 

.6134 

1 . 9 6 6 2 

1 . 5 5 7 9 

. 6 0 2 6 

.6671 

1 • o b 2 o 

1 . 4 4 0 5 

.>311 

.3153 

1.399b 

1 . 30b 9 

, 4 3 y / 

. 4 6 b 0 

1.1172 

1.1653 

. 3 2 6 b 

. 3 9 5 0 

. b 35 1 

1 . Gu 3 2 

• £.161 

. 3 2 6 2 

. 5 5 4 0 

• u 2 o 6 

.1077 

.2514 

. 273 4 

. O 0 O 4 

*»uU2b 

.17 0 1 

-.0066 

. 4 3 2 0 

-.112o 

. 06 21 

- . 2 b 0 0 

. 2 G d 4 

-.1434 

. 0 5 54 

-. 3 6 4 1 

.14 0 6 

- . 17 5 9 

. U 2 7 2 

-.4466 

• 0 6 9 1 

- . 20o0 

-.0014 

- . 5 2'o 3 

-.0033 

-.^597 

-. 0305 

- • 6 0 b o 

-.0773 

-.2710 

-.0599 

-.6862 

- . 1322 

-.3019 

— . 0 B 9 6 

- . 766 7 

- . 22o 1 

CvJ 

• 

1 

-.1201 

- • 6 4 4 1 

-.3051 

- . 3 o 2 5 

- . 150b 

-.9206 

- • 3 b 5 0 

- . 3 y 2 2 

-.Idly 

-.9961 

-.4620 

-.4216 

- . 2 1 3 3 

- 1 . 0 70 7 

- .5 4 1 6 

- • 4 6 6 0 

- . 2 o 2 0 

-1.1835 

-.6654 

- . 5 9 2 o 

-. 4 08 3 

-1.5051 

-1.037 1 

-.7195 

-. 565 9 

-1.6273 

- 1 . 4 3 76 

- • o 4 6 6 

-.7354 

-2 . 1504 

- 1 . 6 o 7 9 

-.9214 

-.6515 

-2.3403 

-2.1374 

- . 9297 

- • b 5 3 3 

-2.3614 

-2.1673 

-.9333 

- . b 5 9 4 

-2 . 3 7 0 4' 

-2.1827 

-.9335 

- • b 6 5 4 

-2.3762 

- 2 • 1 9 b 0 

- • 9 3 6 4 

-.8708 

-2.3 783 

-2.21 19 


Coordinates on L-pIane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 



TABLE 4.--Cont i nued 



ROTOR 

AIRFOIL COORDINATES 


[ Sect 

ion no. 11; 

radius 7*500 in. (19 

050 cm)] 

a Coordi nates, in. 

Coordinates, cm 

X 

Y 

X 

Y 

- • y 4 6 0 

-.7918 

- 2 . 4 0 2 o 

- 2.01 10 

-.'>453 

- • 7 9 6 2 

-2.4010 

- 2.0222 

-.9393 

-.8010 

— 2 . 3 c 8 0 

-2.0344 

-.9346 

- . 8007 

-2.3 737 

- 2 • 0 3 3 0 

-.9291 

-.7985 

- 2 . 3 3 9 9 

- 2 . 0 2 b 2 

-.9234 

-.7947 

-2.3434 

- 2 • 0 1 b 0 

-.915o 

-.7873 

-2.3260 

- 1 . 9 > 9 7 

- . 9 o 7 7 

-.7787 

-2.3054 

-1. 9 7 7 9 

-.0499 

-.7175 

-2 . 1 3 0 7 

- 1 . 6 2 2 ? 

-.7601 

-.6273 

-1.9307 

-1.2934 

-.0707 

-.3427 

-1.7034 

- 1 . 576 5 

- • 5 c 1 5 

—. 468 I 

-1. 4 7 6 9 

-1.1762 

-.4926 

- . 388 3 

-1.251 1 

-.9661 

- . 4039 

-.3177 

- 1 . 0 2 60 

-.8070 

-.3136 

-.251 1 

- • 0 0 1 5 

-.6577 

-.2273 

- . 1880 

- . 5 7 7 7 

- • 4774 

-.1396 

- . 1 2o 0 

-.3546 

- . 3 2 2 

-.0520 

-.0709 

-.1321 

-.1600 

. 0 35 2 

-.0162 

. 089 5 

-.041 1 

. 09 3 4 

.0190 

. 2 3 71 

. 04 61 

• 1 1 6 O 

. 034 2 

.3017 

• 0 6 6 8 

.1443 

. 049 3 

. 3 6 b 5 

.1251 

.1699 

. 064 3 

.4314 

.1o33 

.1955 

. 07 92 

. 496 4 

.2011 

.2211 

. 0940 

.3616 

• 2 3 d b 

. 2 4 60 

• 1 0 8 8 

. 6269 

. 27 6 3 

. 2 7 2o 

.12 3 5 

. 69 2 3 

. 3 1 55 

• 29 6 4 

.1381 

. 7 5 7 8 

. 3 507 

. 3 2 4 2 

.1527 

. 0 2 3 4 

. 38 7 7 

. 3 350 

.15 8 8 

. 6509 

. 4 0 3 3 

• 4 4 6 5 

. 2 209 

1 . 13 4 0 

. 3 6 1 0 

. 5591 

. 282 1 

1.4199 

.7166 

.o728 

. 34 22 

1 . 70o5 

.8690 

. 78 7 9 

. 4 006 

2.0012 

1.0175 

.9192 

. 4 6 3 9 

2.3 3 4 7 

1.1782 

.9311 

• 4 6 9 6 

2.3649 

1.1927 

. 9 37 5 

. 4 7 32 

2.381 1 

1 . 20 1 8 

. 9427 

. 4 7 7 3 

2.3945 

1.2123 

. 948 3 

. 4862 

2 . 4 0 B 6 

1.2349 

. 94 85 

• 4 6 8 1 

2.409 1 

1.2398 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 



TABLE 4.-- Cont i nued 


ROTOR AIRFOIL COORDINATES 

[Section no. 1 1 --cont i nued ; radius 7*500 in. (19*050 cm)] 


a Coord i nates 
X 

, i n. 

Y 


a Coord inates, cm 
X Y 

• y 4 a 2 

.4.9 0 4 

2 

• 4 0 o.4 

1 . 2 4 5 5 

. 9 4 0 2 

. 494 0 

2 

. 4 0 32 

1. 2 54 7 

.9423 

. 4 96 9 

2 

. 5 9 3 3 

1 . 2 c 2 1 

• y 3 6 9 

. 49d9 

2 

. 5 7 9 7 

1 . 2 o 7 1 

. y 3 0 t 

. 499d 

2 

• 5 o £ b 

1 . 2 o 9 4 

. 9 22 7 

. 4 9 9 6 

2 

• 3 4 5 6 

1 .2o90 

. 9 1 0 d 

. 4 96 1 

2 

• 51^4 

1 . 2 6 3 2 

. o 0 3 7 

• 4 c 4 7 

2 

. u 4 1 3 

1 . 2 5 1 u 

.eO/2 

• 4 o 1 7 

1 

• 7 4 U 4 

1 • 1 7 2 o 

.5673 

.4510 

1 

• 4 4 0 6 

1 .1*9 4 7 

• 4 4 9 o 

. 3 9 2 6 

1 

.14 2 4 

. 9 9 7 1 

• 3 3 2 d 

. 3 4 6 4 


• o431 

. e 7 9 a 

.2161 

. 292 3 


. 5 4 6 9 

. 7 4 2 4 

. 09 9 9 

. 2 302 


. 2 3 3 6 

. 564 7 

-.0160 

.159b 

- 

. 04 0 5 

• 4 0 5 9 

-.1314 

. U 807 

- 

. 3 3 5 d 

. 2050 

-.1759 

. 0466 

- 

. 4 4 6 7 

.1184 

-.2067 

. 0229 

- 

.5249 

• 05 8 ^ 

- . 2 37 1 

-.0012 

- 

. 602 3 

-. 0 0 3 0 

- . 26 7 2 

- . 0256 

- 

. 6 7 d 7 

- . 0 6 5 4 

-.2970 

- . 0507 

- 

. 7 3 4 2 

- . 1 2 o 7 

- . 3 2 6 3 

-. 076 1 

- 

. 82 69 

- . 1 y 3 2 

-.3554 

-.1013 

- 

. 9 0 2 6 

-.2586 

-.3841 

-.1260 

- 

. 9 7 5 5 

-.5250 

-.4124 

-.1545 

- 1 

. 04 7 5 

- . 5 9 2 5 

-.4405 

-.1615 

- 1 

. 1 1 o7 

-.4609 

-.4410 

- . 1620 

- 1 

.1200 

-.4621 

— . 4 c 7 4 

- . 228 0 

- 1 

. 2 3 6 0 

- . 5 7 9 U 

-.6105 

- . 5 5 79 

- 1 

• 55 06 

-. 909 1 

-.7340 

-.501 1 

- 1 

• o 6 4 3 

- 1 . 2 7 2 6 

-.8579 

-.6596 

-2 

. 1 7 o 9 

- 1 . 67 5 7 

- • y 3 0 o 

-.7615 

-2 

• 36 41 

- 1 . 934 3 

- . 9 3 d 9 

- . 7729 

-2 

. 3 d 4 6 

- 1 . 96 30 

- . 94 27 

-.7793 

-2 

. 39 4 3 

- 1 . 9 79 4 

-.9451 

-.7657 

-2 

. 40 03 

- 1 . 99 56 

-.9460 

-.7915 

-2 

. 40 28 

- 2.0103 

a Coordi nates on 

L-plane blade sections 

defined by the 

intersection 


of the blade with planes tangent to cylindrical surfaces. 



TABLE 4.--Cont i nued 


ROTOR AIRFOIL COORDINATES 
[Section no*. 12; radius 6-641 in. ( 1 6 . 868 cm)] 


Coordinates, in. 

a Coord i nates, cm 

X 

Y 

X 

Y 

- . y 4 4 y 

- .69d7 

-2 . 4000 

- 1 . 7 7 0 J 

-. 94 4 2 

-.7013 

-2 . 398 2 

-1.7612 

-.9362 

-. 706 3 

- 2 . 3 8 2 9 

- 1 . 7 9 4 0 

-.9333 

- . 70o 1 

-2. 3706 

-1.793? 

-.9277 

-. 704 1 

-2.5563 

- 1 . 1 O i V 

- • 9 2 1 o 

- . 7005 

-2.3413 

- 1 . / 7 9 2 

-.9142 

-.6936 

-2.3220 

-1.7618 

-.9064 

-.6859 

-2.5022 

- 1 . 7 4 2 2 

-.opliti 

-.6310 

- 2 . 1 6 1 0 

- 1 • 6 0 2 6 

-.7644 

-.5504 

-1.9416 

-1.5981 

-.07 6'? 

-.4755 

-1.7233 

-1.2077 

-.592 9 

-.4058 

-1.5058 

- 1 • U 3 0 o 

- • ? 0 7 o 

-.3410 

- 1 . 289 3 

- • a 6 e G 

- . 4 2 2 a 

- . 2807 

-1.0737 

-.7130 

-. 3 3 o 2 

-.2249 

-.8591 

-.5711 

-.2541 

-.1730 

-.6453 

-.459; 

-.1703 

-.1249 

-.4525 

- . 5 1 7 2 

- • 0 d 6 9 

- . 0803 

- . 220b 

-.2039 

- . 004 0 

- . 0 390 

- . U 1 0 0 

-.0990 

. 05 2 3 

-.0129 

.1329 

-.0327 

. 0 7 6 2 

-.0021 

.1935 

- . 005 2 

.1002 

• 0 0 O 6 

• 2 5 4 4 

.0218 

.1242 

.0191 

.3154 

• 0 4 a 5 

• 1 4 d 3 

. 02 9 ? 

. 376 5 

. 07 4 8 

.1724 

. 0 39 7 

. 4 3 7 a 

.1008 

.19 6 6 

. 0499 

. 4992 

. 1 i. 6 6 

. 2206 

. 0599 

. 56 0 7 

.1521 

. 2 4 5 0 

. 0698 

. 622 3 

.1773 

. 269 3 

. 07 9 7 

. 6840 

. 2024 

. 3005 

. 0926 

. 76 3 1 

. 235 1 

. 4 24 6 

.1403 

1.0784 

. 356 4 

. 5505 

.1850 

1 . 398 1 

. 469 7 

• 6 7 a 3 

. 2 260 

1.7229 

. 57 39 

• b 0 8 3 

. 2627 

2.0530 

. 66 7 3 

. 9 56 3 

. 2988 

2.4366 

. 7 59 0 

.9711 

.3019 

2 • 4 6 6 4 

. 7667 

. 9 7 7b 

. 3039 

2.4634 

.7716 

• 9 a 4 3 

. 3070 

2.5000 

. 7 79 7 

. 99 2 1 

.3149 

2.5200 

. 7999 

.99 30 

. 3 1 69 

2.5223 

.6100 


Coordinates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 



TABLE --Continued 



ROTOR AIRFOIL 

COORDINATES 


[ Sect ion 

no. 1 2--cont inued ; radius 6.641 in. 

(16.868 cm) 1 

a Coord i 

nates, in. 

a Coord i nates, cm 

X 

Y 

X 

Y 

.9930 

.3193 

2.5222 

• o 1 09 

.991b 

. 3 2 3 6 

2 . 5 1 90 

. d 2 1 d 

• 9 8 8 5 

. 32 75 

2.5 106 

• o 3 ? 7 

.>933 

. 3 3 08 

2.4974 

• o 4 0 2 

• 9 7 6 -j 

. 5 3 3 4 

2 . 4 o U 3 

• a 4 6 7 

» 9 b b 6 

. 3 3 4 9 

2.4602 

• d 1? u 7 

.9359 

. 33b 1 

2.4279 

• d > 3> 0 

• d4 1 2 

. 3 4 6 4 

2 . 13 6 5 

. a 7 9 9 

.714b 

. 3 4 79 

1 .8154 

• g c 3 6 

• 3 b 9 3 

. 339b 

1.4967 

• 3 o 2 4 

• 4 o 4 b 

.3216 

1 . 1 dOO 

• O 1 Co 

. ^407 

. 294 1 

• 8 6 5 4 

. 7 4 7 0 

• lei 1 7 6 

. 257 1 

. 5527 

o 

•*\ 

o 

. 

. 0952 

.2104 

.2416 

. 5 3 4 4 

-. 0267 

.1539 

-. 06 76 

. 5 9 08 

- • 1476 

. 087 1 

- . 3753 

.2211 

- . 2 0 7 6 

• 0 4 8 o 

-. 52 72 

.1*35 

- • 2 3 6 4 

. 0298 

- • b 0 0 4 

. 0 7 5 7 

- . 2 b 5 0 

.010b 

- • b 7 2 9 

. 02 7 0 

-. 29 3 2 

- . 0089 

-.7448 

-.0226 

-.3212 

- . 0289 

- . b 1 59 

-. 07 3 3 

-.3409 

- . U 4 9 3 

-.8862 

-.1251 

-.3763 

-.0701 

- . 9 55 7 

- . 17 7 9 

-.4033 

-.0913 

-1.0244 

-.2319 

-. 4 300 

-.1130 

-1.0922 

-. 287 0 

-.4564 

-.1352 

-1.1591 

-.3434 

-. 4 59 7 

-.1379 

-1.1675 

-.3503 

-. 504 1 

-.1772 

-1.2604 

- . 4 500 

-. 6 2 2 1 

-.2914 

- 1 . 5 o 00 

-. 7 4 01 

-. 7 4 06 

-.4217 

-1.8809 

-1.0711 

-.6596 

- . 5700 

-2.1833 

- 1 . 44 78 

-. 9 2 98 

-.6668 

-2.3615 

-1. 6 9 3 6 

-. 9 3 75 

- • 677b 

-2 . 38 1 3 

-1. 7209 

-.9415 

-.6641 

-2.3912 

-1. 7 3 7 6 

-. 94 4 0 

- . 6906 

-2 . 3976 

- 1 . 7 54 1 

-. 9 4 4 9 

-.6965 

-2.4000 

- 1 . 76 90 


a Coord i nates on L-plane blade sections defined by the intersection 
of the blade with planes tangent to cylindrical surfaces. 




Figure 39---Typical L-Section. 
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Figure 40. --Blade Average Foil Stress vs Radius. 
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The midspan damper stress distributions are shown in fig. 42. The damper 
was made as thin as possible for aerodynamic reasons and a slight taper with 
thickness increasing from the overhang to the blade junction was required to 
minimize the stresses. The maximum stress at 100 percent speed is 80 ks i 

(552 x 10^ N/cm^) , which is less than the allowable 102 ks i (704 x 10^ N/cm^) . 


The stacking axis displacements and blade untwist are shown in figs. 43 
and 44, respectively. The tilt optimization of the blade at the aerodynamic 
design speed has reduced these displacements to the small magnitudes shown in 
the figures. 


Blade v i brat ion . --The blade interference diagram for the first three modes 
is shown in fig. 45. A frequency band is shown, based on the effective lock-up 
range that can be obtained from the midspan damper. The midspan dampers were 
designed to lock up at 50 percent speed. The midspan becomes effective prior 
to a two-per-revolution (2 E) interference with the first flexural mode of a 
nonsupported blade. As shown in fig. 45, no blade resonances are expected at 
design speed. The first flexure mode occurs between 3E and 4E, the first tor- 
sion mode between 4E and 5E and the second bending mode between 6E and 8E. 
Resonances will occur at speeds below design speed, but are not expected to 
result in high stresses. The effect of blade dimensional tolerance on blade 
frequency is minor (± 1.25 percent). 


Fig. 46 is a modified Goodman diagram showing the allowable vibratory 
stress as a function of steady (mean) stress. The foil midchord at the midspan 
damper has been chosen as the critical point for vibratory stress. The allow- 

2 2 

able alternating vibratory stress (30 ksi or 207 x 10 N/cm ) is shown in 
fig. 46. 


Mode shapes for zero speed and 14 100 rpm (1476.6 rad/ s) are shown in 
fig. 47. The second flexural mode is a complex shape and cannot be considered 
a true second flex. Nodal lines are drawn from the results of displacement 
vectors obtained from the vibration program. The measured mode shapes based on 
a bench test are presented for comparison for the zero-speed case in the appro- 
priate figures. Both the measured and calculated results were with the blade 
dovetail and midspan damper clamped. Good correlation is shown between the 
measured and calculated results. 
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Equivalent stress ks i (N/cm ) 

N = 14 100 r pm ( 1476. 6 rad/s) 
Midspan damper radius = 12.0 in. 



Figure ^2. — Blade Midspan 
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Blade flutter . --Flutter is a self-excited blade vibration that occurs 
primarily in the torsional mode or coupled bending torsional mode. 


The effect of the midspan damper location on the torsional flutter was 
initially investigated by considering the blade to be split into two parts 
built in at the damper. One part is the portion extending radially inward 
from the damper and the other part extends radially outward from the damper. 
The flutter parameter versus midspan damper radial position for each of these 
independent blade parts and also for the complete blade, is shown in fig. 48. 
The optimum position of the midspan damper is 12.4 in. (31.5 cm) as can be seen 
from this figure. If -the midspan damper could be considered to lock the foil 
with 100 percent rigidity (cantilevered outboard from the damper), a flutter 
parameter of 1.6 could be obtained with a single damper located at the opti- 
mized radius. The achievable rigidity of the damper as shown by the cross- 
hatched region in fig. 48 can produce a flutter parameter between 1.1 and 1.2, 
where the flutter parameter is defined as ^C/V'. A flutter parameter of 1.6 
is considered conservatively safe and blades with a parameter as low as 1.0 
have operated successfully. 


cu = blade frequency, rad/s 
C = chord, ft (or cm) 

V 1 = relative velocity, ft/s (or m/s) 

The damper was positioned at a 12.0 in. (30.48 cm) radius giving a 
flutter parameter in torsion between 1.1 and 1.2 (fig. 47) at 100 percent 
of normal operating speed. This final position was chosen as a compromise 
between anticipated flutter stability and midspan damper Mach number. Posi- 
tioning the midspan damper out beyond 12.0 in. (30.48 cm) would have placed 
it at a radius at which supersonic through-flow was expected and would have 
compromised the achievement of this design flow condition. 

Fig. 49 shows the relationship of C/V versus radius. As can be seen, 
it is essentially a constant that makes the flutter parameter independent 
of the radial reference station used. The tip section was used as the 
reference station in the calculations. 

Rotor Disk 

Disk stresses . — The tangential and radial stresses versus radius are 
shown in fig. 50 * The bore and web stresses are within the allowable values 
given in the design criteria. The disk average tangential stress is 64 100 
2 2 

psi (421 x 10 N/cm ) at the overspeed point of 14 100 rpm (1476.6 rad/s). 


Ill 



Frequency, Hz 


Midspan damper location- 


• Optimum midspan damper location 


Frequency of part of blade 
outboard of midspan damper 
(perfectly rigid midspan) 
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Figure 48* — Blade First Torsional Frequency and Flutter 
Parameter vs Midspan Damper Location. 
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Figure 49. --Blade Chord Divided by Relative Velocity 


Stress, N/cm 








Burst and shedding speeds . --The minimum shedding speeds for the foil 
and dovetail and minimum disk burst speeds were determined for the worst 
tolerances and maximum temperature conditions. These shedding speeds were: 


B 1 ade 

foi 1 

18 

Blade 

shank 

22 

Disk 

shank 

18 

Disk 

burs t 

18 


**00 rpm (1927 rad/s) 
900 rpm (2400 rad/s) 
500 rpm (1938 rad/s) 
100 rpm (1895 rad/s) 


The required minimum burst and shedding speed is 17 600 rpm (1843.1 
rad/s), which is 125 percent of the mechanical design speed. The disk, 
which is the most critical element, exceeds the minimum required shedding 
speed by 3.0 percent. 


Radial growth . --The radial growths of the hub and tip leading and 
trailing edges are summarized in table 5* This table presents the centri- 
fugal, thermal, and total growth separately at the mechanical design speed. 
As expected, the maximum radial qrowth of 0.0432 in. (0.11 cm) occurs at 
the blade tip leading edge. 


TABLE 5 

HUB AND TIP RADIAL GROWTH 



Hub 

Tip 


i n . 

cm 

i n . 

cm 

Centri fugal growth , 
N = 14 100 rpm 





Lead i ng edge 

0.016 

0.0406 

0.034 

0.0864 

Trailing edge 

.016 

.0406 

. 

.028 

.0711 

Thermal growth, disk 
temperature = 150°F (33^°K) 
blade temperature = 250 F 
(395°K) 





Lead i ng edge 

.0025 

6.35 x 10“ 3 

.0092 

.0234 

Tra i 1 i ng edge 

.0028 

7.11 x 10" 3 

.0083 

.0211 

Total growth 





Leading edge 

.0185 

.0470 

.0432 

.1097 

Trailing edge 

.0188 

.0478 

.0363 

.0922 



Blade/Disk Attachment 


Attachment stresses . — The dovetail attachment stresses in the blade to 
disk have been summarized in table 6 . Neck, tang, and combined stresses 
are shown. The combined fillet stress is a combination of tang bending 
and neck tension excluding stress concentration effects, while the maximum 
fillet stress which effects the low cycle fatigue life of the attachment 
includes the stress concentration factor. All the calculated stresses 
show a positive margin of safety based on the design criteria given pre- 
viously. The disk, which is the weakest link as a result of using existing 
tooling, has a 0.03 margin of safety on the maximum fillet stress. The 
allowable stress levels are based on a 10 000 -starts requirement and can be 
deviated from for this test rig application. 

Stator Vane 

Geometric def in it ion . --The geometric summary of the final stator vane 
design is presented in tabular form in table 7 - 

Steady- stress .--The maximum airfoil steady-state stress is 6000 ps i (41.4 
2 2 

x 10 N/cm ) and is induced by gas bending loads. This stress is well below 
the al lowable. 


Vi brat ion . — The calculated natural frequencies of the fan stator are shown 
in fig. 51- The only exciting force will come from the 40 fan rotor blade 
wakes. Whenever the 40-exc i tat ions-per- revol ut ion line crosses a natural fre- 
quency, a possible resonance condition with an amplification of the vibratory 
stress level exists. Examination of the excitation diagram reveals possible 
resonances throughout the operating speed range; however no resonance problems 
are expected. The first four natural frequencies will not be excited 
above 5000 rpm (522 rad/s) and can be passed through rapidly during rig 
running without causing any problems. The allowable vibratory stress 

2 2 

based on a modified Goodman diagram is 34 ks i (234 x 10 N/cm ) as shown 
in fig. 52 . 

Stator vane flutter . --The frequency parameter, \ = coC/V', was calculated 
for torsional flutter. The calculated value of 2.75 is well above the con- 
servative minimum of 1 . 6 , indicating stable operation. 

Spool Critical Speed 

The critical speed analysis of the fan spool and facility drive system 
shows that no resonance problems will occur within the rig operating speed range. 
The first critical speed of the system with synchronous precession (whirl ratio 
= 1 ) is 16 370 rpm (1714.3 rad/s), which provides a 16 -percent margin over the 
maximum speed. The second critical speed is calculated to be 23 050 rpm 
(2413.8 rad/s). 


116 



TABLE 6 


BLADE-TO-DISK DOVETAIL ATTACHMENT STRESSES 



Blade 

Disk 





Margin of 




Margin of 


A1 lowab le 

Nomina 1 

Peak 

Safety 

Al lowable 

Nomi nal 

Peak 

safety 

Neck tension or 

87.4 ? 

36.4 

46.3 ? 

0.88 

87.4 ? 

61 .0 , 

76.3 0 

0.15 

tang bending 

(603 x 10 ) 

(251 x 10 ) 

(319 x 10*) 


(603 X 10 ) 

(421 x 10 ) 

(526 x 10 ) 


Comb i ned fillet 

121.3 , 

69-1 , 

86.5 ? 

.40 

121.3 , 

94.7 ? 

114.5 

.06 


(838 x 10 ) 

(476 x 10 ) 

(596 x 10*) 


(838 x 10*) 

(790 x 10*) 

(790 x 10 2 ) 


Max imum f i 1 let 

181.4 

0 

127.2 

.43 

181.4 

0 

177.0 , 

.03 

(low cycle fatigue) 

(1252 x 10 ) 

0 

(877 x 10*) 


(1252 x 10*) 

0 

(122 x 10- 5 ) 


Beari ng 

121.3 , 

47.7 ? 

71.8 

• 69 

121.3 , 

47.7 2 

71 .8 , 

• 69 


(838 x 10*) 

(329 x 10 ) 

(495 x 10*) 


(838 x 10*) 

(329 X 10 ) 

(495 x 10*) 


Shear 

50.7 , 

23.1 , 

27.0 , 

.83 

50.7 , 

21.9 , 

25.6 , 

• 93 


(350 x 10*) 

(159 x 10*) 

(186 x 10 ) 


(350 x 10*) 

(151 x 10' ) 

1 

(177 x 10*) 



NOTES: 

(1) All stresses in ksi and (N/cni ) 

(2) All stresses calculated at 100 percent speed 14 100 rpm (1^76.6 rad/s) 

(3) Margin of safety is based upon the larger of (1.2 X nominal) or (calculated peak) 

Margin of safety = (allowable stress/calculated stress) - 1 

(4) Allowable stresses include: 

5~percent reduction for surface finish and leading rate 
5-percent reduction for broach angle 



















































TABLE- 7 


STATOR GEOMETRIC SUMMARY 

[L-Plane Blade Sections, defined by the intersection of the blade and planes tangent to cylindrical surfaces 
Blade Height^ 6.076 in. (0.15^2m) Number of vanes = ^5 Airfoil Section 

Hub/tip ratio = 0.560 Aspect ratio = 3-10 Double circular arc 


D ^ A i no * ^ • 

13.499 

12.891 

11.980 

11.373 

10.765 

10. 158 

9. 550 

8.639 

8.031 

InoIQI Us 

m 

0.3429 

0.3274 

0. 3043 
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APPENDIX A 


THE CHARACTERISTIC COMPAT IB I L 1 TY EQUATIONS 


The purpose of this appendix is to derive the characteristic compatibility 
equations which apply to the flow in the outboard regions of the rotor. The 
derivation is divided into the following sections: 

(1) General discussion 

(2) Assumptions 

(3) Governing differential equations 

(4) Governing differential relationshi ps 

(5) The compa ti bi 1 i ty equations in P -3 form 

(6) The compatibility equations in § + , I form 

(7) Derivatives of the unit vector \ 


GENERAL DISCUSSION 

Consider the steady, supersonic, planar flow of a gas in some shockless 
region R. At each point P in R, two directions (referred to as the positive 
and negative characteristic directions) may be defined. The positive charac- 
teristic direction is oriented at an angle of +M- and the negative characteristic 
direction at an angle of -|*i with respect to the velocity vector at P, The angle 
M- is the Mach angle, defined as 

s in H m (A-1) 


A curve passing through point P, and tangent everywhere along its length 
to the local positive characteristic direction, is referred to as the positive 
character! Stic through P. Similarly, a curve through P everywhere tangent to 
the negative characteristic direction is called the negative cha rac ter i s t i c 
through P. The positive and negative cha rac ter i s t i c curves through P are 
denoted by C* and C~, respectively. If the subscript p is not specified, 

then C + and C refer to arbitrary positive and negative characteristics. 


Consider again the above flow field, with the additional constraints that 
the field be adiabatic and have constant entropy and total enthalpy, and that 
the gas be an inviscid, thermally perfect gas of constant composition. The 
theory of characteri sti cs as applied to the mathematical model of this flow 
field indicates that the field possesses a rather surprising property. This 


property is that the quantities and 5 , defined as = V | - ■ ■ — and § = 
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are invariant on the positive and negative characteristics, respectively, 
independent of how these cha racteri s ti cs weave and wind through the field. 

This property may be stated in differential form in the following manner 

d§ + = 0 on C + (A- 2 ) 

d£ = 0 on C (A-3) 

or more compactly as 

d§~ = 0 on C“ (A-4) 

Eq. (A-4) is referred to as the characteristic compatibility equation for this 
particular flow field. 

The extreme simplicity of eq . (A-4) allows flow problems which satisfy 
the above field constraints to be analyzed in a simple, though tedious, 
manner by the so-called "method of characteristics." The technique is, in 
fact, simple enough to be done by hand. 

Many flow fields which violate some of the above constraints can still be 
solved using the method of cha racteri st i cs ; however, the compatibility 
equations and consequently the solution technique become more complex. In 
general, the additional complexity is sufficient to make hand calculations 
impractical, and recourse must be made to high-speed digital computers. 


ASSUMPTIONS 

In deriving the compa t i b i 1 i ty equations which apply to the stream tube 
flows in the outboard regions of the rotor, the following assumptions are made 

i 

(1) The fluid is constrained to flow within an axisymmetric stream-tube 
of specified geometry. 

(2) The fluid is an inviscid, thermally perfect gas with constant 
composition and zero thermal conductivity. 

(3) The flow field is steady and fully supersonic relative to a 
coordinate system which is fixed to the rotor. 

(4) Electric, magnetic, gravitational, and thermal radiation effects 
are neglected. 

Assumption (1) above is particularly significant in that it reduces the 
fully three-dimensional problem to a quas i -three-dimensiona 1 problem. As a 
consequence, effects associated with the blade geometry on other stream tubes 
will not be accounted for except indirectly, insofar as they affect the 
local stream-tube geometry via the axisymmetric flow analysis. 
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GOVERNING DIFFERENTIAL EQUATIONS 


Within the accuracy of the above assumptions, the flow in the outboard 
regions of a rotor is governed by the following differential equations 
relative to a coordinate system fixed to the rotor. 

Conti nu i ty : 

■t V’ - VP + V • V' = 0 (A-5) 

p 

As presented by Wu (ref. 9) the momentum equation is written as follows: 

(V’ * V ) V' + 2u> A V’ - U> 2 7 + = o 

P 

Also, from ref. 9, eq.(lAa) for steady isentropic flow, 

V I = V’ A (V A V') (h-7) 

where I is the enthalpy function: 

I = H - u)(V Q r) 

Taking the dot product of V' on both sides of eq. (A- 7 ) yields 

V' -VI =0 ( A “8) 

Owing to the fact that V' * V' a (Va V') = 0 

4 From the equation of state and the conditions of isentropy, the 
following is obtained: 


P 

P 



(A-9) 


where the subscript, o, refers to an arbitrary reference state. 


For the purposes of this analysis, the system of eqs. (A*-5) through 
■(A- 9 ) , is more conveniently handled in a less conventional form developed 
as fol lows : 


From eq. (A -9), 


l 

Y 




(A- 10) 
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Defining a pressure function P as 



and taking the gradient of eq. (A-10) yields 

=vp 

p 

Substitution into eq. (A-5) yields 

V 1 • VP + V • V 1 = 0 


From eq . (A-1 1 ) 


\ 

or 


VP 


1 £E 

V P 


but 


therefore 



S£ 

P 



IE 

P 


££ = C 2 VP 


where C is the sonic velocity. Hence, eq. (A-6) becomes 

(V 1 • V) V 1 + 2u) A V' - u) 2 r + C 2 VP = 0 


(A-1 1 ) 


(A- 12) 


(A- 13) 


GOVERNING DIFFERENTIAL RELATIONSHIPS 

For convenience, eqs. (A-12) and (A- 13) are converted into a system of 

differential relationships relative to the unit vector direction i , i , and 
_ _ _ sn _ 

i.., where i is parallel to the velocity at each point, i is normal to i 
N s _ n _ s 

and lies in the axisymmetric stream surface, and i^ is normal to both i and 

i (see fig. A-l). Eqs. (A-12) and (A- 13) are resolved into these intrinsic 
components as follows: 

S i nee y i _ j y i 
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Figure A- I. -- Coordinate System - Derivation of Characteristic 
Compatibility Equations. 



and 


„„ - dp — dp - dp 

1 s ds + 1 n dn + 1 N dN 


we obtain 


— dp 

V' -VP = V gjf 


(A- 14) 


Also, 


V- V' = 


i. • §7 <i„V') + I. • sr (1.V-) + i, 


dn 


SN <i s V ' ) 


. avv 4 . S + , . S) 

ds ' n dn 'n dN ' 


(A- 15) 


Egs. (A-14) and (A- 1 5) substituted into eq. (A- 12) yield 

...dp dy 1 — d* s — ^ 

V ds + ds + V \ ' n ' dn 1 N ' dN J 


(A- 16) 


The first two terms of eq. (A- 13) are as follows: 

(v* • V)v> = V’ §J (T v) 


- v ’ (', Dr + v ' 57 s ) 


(A- 17) 


di 


= I 


V !¥■ ♦ V 2 ^ 


d s 


ds 


2u>AV' = 2 


1 1 1 k i 

s n N 


cjd u) u) , 

n N 


V' 0 0 

= 2v <r n « N - r N « n ) 


Accordingly, eq . (A- 13) becomes 


* dj 

d\/ 1 2 q — — 

LV' jdr + V' TT + 2V* (i W.. - i.,0) ) 
s ds ds v n N N n' 


7 ID 2 


r + c 2 (t If + 7 |£ + 7 |e) = 0 

\ s ds n d N ' n %H/ 


(A- 18) 


(A- 19) 
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In eqs. ( A- 1 6 ) and (A-19) , s, n, and N refer to displacements in the i s , 

i ^ , and i^ directions, respectively. Eqs. (A- 1 6 ) and (A-19) do not represent 

a set of differential equations with s, n, and N as independent variables. 

Since the metric coefficients have not been included in the derivations, these 
three variables do not define a coordinate system. More correctly, eqs. ( A- 1 6) 
and (A-19) should be viewed as a system of differential relationships which is 
valid at each point in the field. 


The expressions for i and derivatives of i g are presented in terms of 
the unit vectors i s> i n » and i^ under the headings of Derivatives of the Unit 
Vector i at the end of this appendix. In terms of these relationships, eq. 
(A- 16 ) reduces to 


*'f^Mt+ sin V° ,g + cosetg|.o 


)■ 


and eq. (A-19) may be rewritten as 

/ 

7 s ^'§ 7 '- w 2 r cos 3 s in a + C 2 || )+ V 


3 i 
3s 


+ i (2V ,<U sin ot + ^ r sin 3 sin & + C 2 


(A-20) 


- ( a 2 2 dp\ 

+ i N ^- 2V 1 & sin 3 cos a - u> r cos a + C $^1 = 


(A-21) 


Substituting the expression for the derivative of i in terms of the unit vectors 

into eq. (A— 17) and resolving the resulting equation into its components yields: 

/ 




cos 3 s in Of + C 


2 dp 

Fs = 


0 


(A- 22) 


- 1 ^ ( d3 + sin Of sin 3 \ 
\3s r / 


+ 2V’u ) sin Of + w 2 r sin 3 s in a + C 2 « 0 


V 2 (cos e |f - iinLLc°S^\ . 


2V* uj sin 3 cos ot - u) 2 r cos a + C 2 ||j 


(A- 23) 


0 

(A- 24) 


Since the flow is constrained to flow in a specified ax isymmetr ic 
stream-tube eq. (A- 24) is extraneous, serving simply to define the pressure 
gradient necessary to support this stream tube. 
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THE COMPATIBILITY EQUATIONS IN P-P FORM 

Multiplying eq. (A-20) by V 1 and subtracting it from eq. (A-22) yields, 
(after suitable algebraic manipulation of the resulting equation and 
eq. (A- 23 ) ) : 


M - 1 


dp 

3? 


BP 


= - K, 


(A- 25) 


1 Bp Bp 

- 2 ^ + SI = - K 2 (A- 26) 

M 


where 


s i n oi o . U 1 . o Ba 

K. = cos P [ 1 + — 9 + cos P r-; 

1 r l v ,2^ BN 


= sin a 
K 2 r 


in P 




Noting that by definition 


s in M = ^ 


cos U = 


Vm 2 - 1 


Eqs . (A- 25) and (A- 26 ) can be rewritten as 


s in 2M» dp .. dP . 

2 C0S “ d? + s ,n ** = " s ,n ^1 

S - n ^ sin M* + cos 4 ^r- = - cos I^K 
2 B n ds 2 


Adding and subtracting eqs. (A- 3 1 ) from (A- 3 2 ) yields 


(A- 27) 
(A- 28 ) 

(A- 29 ) 
(A-30) 

(A-31) 

(A-32) 


s in 2 ji 


cos 


dP 

U dl 1 


s in Mt 


rr\ibos H It ± 




ds 


s in l-i 


BP\ 

dn/ 


(s in 


1 *K 1 ± 


cos 


kk 2 ) <*-33) 


The parenthetical expressions on the left hand side of eq. (A — 3 3 ) represent 
direction derivatives in the positive and negative characteristic directions. 
Hence, in terms of the notation of eq. (A-4) , eq. (A-33) becomes 


sin 2|x 
2 


dP ± dP = 


- K, dl on 
4 


+ 

C" 


(A-34) 
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where dl is a differentia] displacement in the appropriate characteristic 
direction and is given as 

sin (P ± M.)^1 + ± cos p (A-35) 

Eqs (A-34) and A-35) represent the characteristic compatibility equations in terms 
of the pressure variable P and the flow angle P. For development of the blade 
sections, it is more convenient to express the compatibility equations in 
terms of the variables 5 and The conversion from the P - P to the - § 

form is accomplished below. 


K _ cos P da + sin ff 
k M dN _ r 


THE COMPATIBILITY EQUATIONS IN § + , §" FORM 
The pressure function P, eq. (A- 11), can be expressed as 

P = - log I s — — ) 

V * e \P tr p o / 


2 

Y 


log,. 


(*)• '••• (i) 


where p is the local value of relative stagnation pressure, 
tr 

of P is 


The differential 


dP= Y 


1 VW + 1 dp tr 


Y P 


tr 


tr 


In the f i rst ternl 
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In the second term 


Therefore , 


‘ 


JL 

y Y ” 


• fit 

_/_tr 


p 

l t 

/ 



X 0 : 

r 


■ 1o 9e Po 

_ _x_ 

Y-1 

( ,0 3 e T tr - 

,0 9 e 

1 d p 

1 tr 

1 

d T 


tr 


Y P„ 
tr 

y-1 

T tr 


dP 

1 

diH 2 ) ... 

1 

1 

: ,+i h 

" Y-1 

f v, the 

Prandt 1 

-Meyer angle 

is 

dv = 

n/m 2 " 




dT* 


tr 


2M 


(l ♦ ^ M 2 ) 


(A-36) 


(A-37) 


Noting that 


s in 2M- • .. ^M 2 - 1 

— 1 - = s in M> cos p> = — 

Z M 


Eq. (A-36) yields 


sin-lit dP = . 


- 1 d(M^) sin 2M* r 
2H 2 (l + ^ M 2 ) 2 ( Y '» T tr 


Accordingly, eq. (A-34) becomes 


„± ,, . s in 2U> dT tr „± 

dv ± dP = dl + 777-^77 y on c 

/ t r 


S ince 

eq. (A- 38 ) becomes 


dv ± dP = 2d§ _ 


dT 


K ± . s in 2M» 
dl k 2 (Y - 1 ) T dr dl 


dr (on n 


(A-38) 


(A- 39) 
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The relative stagnation temperature T is 


Therefore , 


Also, 


T tr = T o + 2C 


dT 2 ,, 2 

tr a) r U 


/ '■ A 

( r - r o) 


dr C C r 
P P 


c = (Y - 1)C T 
tr p tr 


Substitution of eqs. (A-40) and (A-4I) into (A-39) yields 


(A- 40) 


(A-41) 


2d$~ 


r 2 
± sin 24 U_ dr 

K 4 2 r _2 dl 


L 


tr 


dl (on C ) 


(A-42) 


The following substitutions are made into (eq. A-35) and into the 
second term of eq. (A-42) 

da = I dh 
dN h dm 

2 2 2 
V ' = M C 


1 (^f) ± = cos (S ± (.) (on C*) 


and' 


sin 2P cos p 
2 “ M 

Finally, eq. (A-42) is written 

dS* = e 1 dl (on (T) 

where 



(A- 43) 


(A-44) 
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DERIVATIVES OF THE UNIT, VECTOR i 

s 

Consider the intrinsic unit vectors, i^, i^, and i^ originating from the 

same point as the cylindrical coordinate unit vectors, i ^ , i^, and i^, 

fig. A-2. • The geometrical re] a t ionsh i ps between the two sets of unit 
vectors are developed as follows. 

From the terminus of the i vector, a normal is drawn to the i , i 

s _ _ z r _ 

plane (component 3 ) followed by a normal drawn in the i^, i ^ plane to the i^ 

vector (component 2 ) and concluded by a segment on the i vector to its 

origin (component 1). (The foregoing components are numbered in descending 
order because the imaginary path followed during this construction is opposite 
in direction to the vector sum def in i ng 7^ . ) Tracing from the origin along 

the foregoing orthogonal segments to the i terminus defines the magnitude and 

sense of the following components. 

Component 3 : 1 q sin P 

Component 2: i cos 3 sin o' 

Component 1 : i cos P cos ot 

Accordingly 

I = i cos P cos a + i cos P sin ot + i’q sin 3 (A-45) 

Similarly, the i unit vector is resolved into i , i , and ig components 
by means of ( 3 ) a normal from the i^ terminus to the ig, i plane, ( 2 ) a 
normal drawn in the ig, i plane to the ig vector, and (l) a segment on the 
ig vector to its origin. As before, the magnitude and sense of the components 
a re : 

Component 3 : -i sin 3 cos 01 

Component 2: -i sin3 sincv 

Component 1 : ig cos 3 

and the vector sum yields 

i = i 0 cos 3 - L sin 3 sin or + i sin P cos a (A-46) 

n B r z 

The 7.. unit vector lies in the i , i plane; its components are 
N 2 r 

i = - i sin oi + i cos ot (A-47) 

n z r 
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The foregoing reasoning applied to the resolution of i , i^, and ig into 
their intrinsic components yields 


1 r 


lA = 


i = 
2 


Taking the der 


cos P sin o' - i sinPsina+L. cos & 
s n N 


s in P + i cos P 
s n 


cos P cos <y - i sin P cos O' - i sin Q? 
s n N 


(A-48) 
(A-49) 
(A- 50) 


vative of i eq . (A-45) with respect to any arbitrary 


variable X (noting that the unit vector, i , is invariant) yields 


S — a • Scy T . q ^P 

SjT “ • ' z cos P sin a gj - i z sin 3 cos a ^ 

t . d . 33 — q Ba 

-i sin P sin O' jr + i cos P cos Q ' -^r 

r °X r °X 


3 i 




— 3p ^ 1 r ~ 1 q 

+ ig cos P ^ + 5 — cos 3 sin a + sin 3 


(A-51) 


and 


di 


da 


= cos 3 (- i 2 sin a + i r cos a) ^ 



+ (- i ^ sin 3 cos ot - i ^ s in 3 s in » + ig cos 3)^ 


3 i 


di 


e 


r a u 

+ cos 3 s in ot + s in 3 


(A— 52) 


Substitution of eqs. (A-46) and (A-^+7) into eq . (A- 50 yields the 

fol lowing : 


Si 

Sx 1 


cos 3 


3a — 33 d * ^ '0 . 

Zx + 5x + sr cos e 5ln “ + 5T s,n e 


(A- 53) 
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Since the variable X operates on i and ig through the polar coordinate, 

0 (i.e., changes in i or i Q comprise rotations about the 2 axis), the follow- 

r w 

i ng are obtained (fig. A-3) : 


Si Si 

r _ r S9_ 

sx se sx 

_ - se 
" '0 sx 


(A- 54) 


5l e 5l e se 

sx se sx 


-r se 

r SX 


(A- 55) 


Substitution of eqs. (A-25b) and (A-25c) into eqs. (A-27a) and (A-27b) yields 

= (T sin P + 7 cos P) || (A ’ 56) 

3X s n ' Sx 

sT, 


and 


'0 




= ~ (i cos P sin a - i sin P s in Qf + i„ cos ol) ■r— 
s n N Sx 


SX s n 

From eqs. (A- 53) , (A- 56 ) , and (A-57) 

Si 


(A-57) 


SX 


- /se . se\ 
■ (sx + sm “ ax) 


+ i„(cos 0 Sa - sin 0 cos a S6 ) 


(A- 58) 


SX 


SX 


The genera] variable X is now identified successively as s, n, and N. 

The partial derivatives, S6., S6_, and S0 are evaluated from geometrical 

Ss Sn SN 

considerations. For an incremental translation As along the s-axis, a transla- 
tion in 6 resu Its: 

As sin 0 = rA6 

Accordi ngl y, 


S6 sin p (A- 59) 

Ss ~ r 


135 





Similarly, for a translation along the n-axis, 

An cos 3 = rA0 

99 _ cos B (A-60) 

9n r 


However, the N-axis lies in the meridional (e, r) plane. A translation along 
the N-axis yields no change in 9: 


99 _ Q (A- 6 1 ) 

9N " U 


Finally, eqs. (A - 58 through (A-61) yield 


9 

9"s 


's -r /93 , s in O' s in 
s = ‘n\9s + r / 


+ 7 r Sin 2 3 

+ ' N C0S P 9i - ~ 


3 cos cy I 
r / 


(A- 62) 


9 i 


_s _ t /93 sin g cos 3 

' “ 'n 9 n + 7 




cos 3 


3a sin 3 cos 3 cos 
9 n " r 




(A-63) 


9 i 
5 

9N 


' 93 + • 

n 9n 1 N 


9a 

9N 


(A- 64) 
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APPENDIX B 

TRAILING EDGE MIXING CALCULATION 


The equations used to calculate the trailing edge blockage and the 
inviscid core flow conditions are derived by relating the core plus boundary 
layer flow to the uniform downstream flow. 

Consider a vanishingly thin control volume enclosing the trailing edge 
region of two adjacent blades as shown in fig. B-1. The fluid density, pres- 
sure, velocity, and flow angle at the discharge plane of the control volume 

(denoted by p , p , yt , and 0 1 , respectively) are assumed to be uniform and 
e e 0 e 

equal to the cor respond i ng section discharge values calculated by the axi sym- 
metric program. The corresponding flow quantities at the inlet of the control 
volume (denoted by p, p, V 1 , and |3 T , respectively) are assumed to be nonuniform. 
The pressure acting on the trailing edge of the blade is assumed to be uniform 
and is denoted by p^. 

The conservation of mass, meridional momentum, and tangential momentum 
applied to the control volume results in the following three equations: 


i 

\ 



pV* cos 3 T dy 


p V' cos 
e e 


0* Y 
e e 


(B-1) 


9Ph 


Y 

(VM + At f 


pV* cos^ 


p dy = 


gp Y 

3r e e 


+ p V« 
e e 


cos 


P’Y 

e e 


(B-2) 


where 



P dy 



pV* cos P' s i n P' dy 


p V' cos P' sin S' Y 
e e e K e e 


(B-3) 
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If the flow in the inviscid core region is assumed to be uniform, and the 
corresponding flow properties are denoted by the subscript c, then eqs. (B-1), 
(B-2) , and (B-3) can be rewritten as 


1 

Y 



pV*cos (3 


:*c 


cos 


y 


dy = F 


(B-4) 


Y 

c 


1 

Y 


o 


pV* 2 cos 2 B' 
P c V, c cos2 P c 


dy 



(B-5) 


Y 


c 


1 

Y 


o 


pV ,2 sin P* cos j dy 
P c V’ 2 sin B' c cos &' c / 


(B-6) 


I 

\ 


where 


F i = 1 - 


p v Y q 
e me e 

P V Y 

c me c 


(B-7) 



(B-8) 


F 


3 


p v V* Y 

e me 6e e 

P V o Y r 
c me Be c 


(B-9) 


F 


P 


P Y 

e e 


P c Y c - Pb< Y e - V 

p V 2 Y 
c me c 


(B-1 0) 


and where V and V’ D are the meridional and tangential components of V’. in the 
m B 

core region, the integrands of the above equations are zero, whereas in the 
boundary layer region, the integrands are nonzero. In order to evaluate these 
integrals in terms of boundary layer parameters, assume that the flow angle in 
the boundary layer is the same as in the core, and change the variable of inte- 
gration from dy to dn where 
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, dn 

dy = 


(B-tl) 


and the upper limit of the integral becomes Y c = Y c cos (3\ 


The integrals in eq. (B-4) and the integrals in eqs. (B-5) and (B-6) 
may then be written respectively as , 

11 (' ' Mf) dn ( 


(B— 1 2) 




(B-13) 


Algebraic expansion of the integrand of 12 yields 


cos 3' 


^) dn+ /« (' '^i) dn (B ' ,4) 


If 6^ and 9* are used to denote the sum of the displacement and momentum 
thicknesses, respectively, of the pressure and suction surfaces at the trailing 
edge, then it follows that the first integral is 


cos P' 


(B-15) 


and the second integral is 


12 = — '—zr ( 6 * + 0 *) 
cos 9' ' ’ 

c 


(B-16) 


Substituting eqs. (B-15) and (B-16) into eqs. (B-4), (B-5), and (B-6) and 
noting the previous assumption that the flow direction is the same for the core 
and boundary layer flow yields 


Y cos P' F 1 
c c 


(B-17) 


6 * + 8 * 

Y cos p' 2 " p 
c c r 


(B-18) 
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6 * + 9 * 


(B-19) 


Y 

c 


cos K 


F 


3 


Finally, simultaneous solution of the preceding three equations with 
appropriate substitutions for the F terms of eqs. (B-7) , ( B -8) , and (B-9) yields 


( B— 20) 


where 



H 


6 * 
9 * 


defines the boundary layer form factor. 


Eqs. (B- 1 7) , ( B— 18) and (B-19) can be handled numerically in the follow- 
ing manner: 

(t) The relative total temperature and' total pressure of the core flow 
are determined from the flow conditions upstream of the control 
volume. In cases where the section is assumed to be operating in 
the unstarted mode, the core relative total pressure is reduced 
from the ideal value by an amount consistent with the calculated bow 
shock loss. 


(2) A discharge "shape factor" H is specified by the designer. If a 
shape factor of 1.5 is assumed to correspond to a flat-plate boundary 
layer, and a shape factor of 2.5 is assumed to correspond to a boun- 
dary layer at the point of incipient separation, then the choice of 

a discharge shape factor between these two values would appear 
reasonable. In the present rotor design, a value of 2.0 was used. 

(3) A value for the overall trailing edge blockage X (the sum of the boun- 
dary layer displacement thickness and the trailing edge metal block- 
age) is assumed. 

(4) With the information from steps 1, 2 and 3, the continuity equation 
in the form 


P 


V 

c me 


0 - M 


p 


v 

e me 


plus eq. (B-20) is used to calculate the core flow conditions and, 
consequently, the values of F^ , F^, F^, and Fp in eqs. (B-7) through 

(B-10). 
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(5) The shape factor is calculated using eqs. (B-17) and (B-19). If the 
value is not equal to the specified value, a new value of blockage 
factor is assumed, and steps 4 and 5 are repeated until agreement is 
achi eved . 

At present, no convergence problems have been encountered in applying the 
trailing edge blockage calculation. 

The flow conditions leaving the control volume for streamline No. 4 as 
calculated from the axisymmetric flow analysis are: 

M =1.094 
e 

3 ^ = 60.24 deg 

V = 624.9 ft/s (190.8 m/s) 

Vg e = 1092.7 ft/s (333.7 m/s) 
p e = 2491 lbf/ft 2 (1.193 X 10 5 N/m 2 ) 

P = 5282 lbf/ft 2 (2.529 X 10 5 N/m 2 ) 

P e = 0.0848 lbm/ft^ (1.358 Kgm/m^) 

T te = 682. 5°R (379. 2°K) 

The core fluid properties (assuming no upstream shock loss) are 

p tc = 5529 lbf/ft 2 (2.647 x 10 5 N/m 2 ) 

T = 682. 5°R (379. 2°K) 

The discharge pitch and the core pitch (trailing edge tail radius = 0.0073 in.) 
a re 


Y g = 2.103 in. (5-342 cm) 

Y c = 2.089 in. ( 5-306 cm) 

With the trailing edge back pressure assumed to be equal to the core static 
pressure and the trailing edge shape factor specified as 2 . 0 , the trailing 
edge blockage and core absolute swirl velocity are 
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X = 6.18 percent 


V(abs) * - 28U5 ft/s (“86.0 m/s) 

The corresponding core flow conditions are 

M c = 1.143 
8 C = 59-35 deg 


V = 664.5 ft/s (202.9 m/s) 

me 

V = 1121.3 ft/s (342.4 m/s) 

P c = 0.085 lbm/ft^ (1.362 Kgm/m^) 

The corresponding boundary layer parameters are: 

6* = 0.0512 in. (0.130 cm) 

9* = 0.0256 in. (0.0650 cm) 


The discharge absolute swirl velocity calculated by the axi symmetric flow 
analysis program is 


V(abs) * - 3 ' 0 -' ft/s I " 34 ' 7 m/s) 

As a consequence, the absolute total temperature of the core flow is lower 
than the corresponding temperature of the downstream flow field. The differ- 
ence is accounted for in terms of the higher total temperature associated with 
boundary layer fluid. 

Similarly, the core static pressure is lower than the static pressure of 
the downstream flow, the difference being the result of the momentum exchange 
involved in the mixing process. 
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APPENDIX C 

LEADING EDGE SHOCK CALCULATI ONS --STARTED MODE 


Once the strength of the leading edge shock is set, either by calculation 
for the outboard section, or by design input for the central sections, the 
problem remains of calculating the shape of the shock in the blade passage and 
the point where the shock intersects the suction surface. The techniques used 
to determine the shock shape and the intersection point are described in the 
following discussion. 


TWO-DIMENSIONAL ISOLATED WEDGE 

Consider a two-dimensional wedge with wedge angle 6 and circular leading 
edge nose radius r in a supersonic flow field at zero incidence (see fig. C-1). 
If the wedge is assumed to extend to downstream infinity, then theoretically 
the supersonic upstream Mach number can only be sustained if 6 is less than 
some maximum attachment angle 6^, which is a function of the specified upstream 

Mach number. Assuming that 6 is less than 6 , then a shockwave will stand 

m 

forward of the finite leading edge of the wedge as shown in the figure. At 
points far removed from the wedge, the shock angle is asymptotic to the ideal 
shock angle for a sharp wedge of the same included angle. In the region of 
the leading edge, the shock gradually becomes stronger until it becomes a nor- 
mal shock on the wedge centerline. Following the technique described in ref. 

4, the shock is assumed to have a hyperbolic shape described by the equation 

_2 _2 2 

X - Y cot o m = 1 (C-1) 

where o m is the shock angle far from the wedge, and X and Y are the X and Y 
coordinates normalized to X Q (fig. C-1). The local shock angle a is dY/dX; 
accordingly differentiating eq. (C-1) yields 

2 2 

cot a = (1 - 1 ) cot a (C-2) 

x 2 

At some point on the shock, the discharge Mach number from the shock is just 
sonic. This point will be called the sonic point. The sonic line dividing 
the subsonic flow in the leading edge region from the supersonic flow far 
from the wedge begins at the sonic point and ends somewhere on the wedge sur- 
face. Specifically, the sonic line is assumed to terminate on the wedge sur- 
face at the point m where the wedge surface angle corresponds to 6^ (ref. 4). 

Obviously, this point lies on the leading edge radius at 


X = X - r sin 5 (C-3) 

me m 
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(C-4) 


Y = r cos 6 
m m 


The symbol r represents the 
addition, the sonic line is 
and oriented at an angle 6^ 

location of the termination 


blade leading edge radius normalized by X . J_n 

assumed to be straight, of normalized length 

+ 90 where 6 is the average of 6 and 6 . The 
a 3 s m 

point in terms of H is 


X 

m 



Y 

m 



+ H sin 6 


- j l cos 6 


a 


a 


(C-5) 

(C-6) 


For two-dimens i ona 1 flow, 
continuity as follows: 


the value of J is approximately related to Y^ by 



where the subscript 00 denotes the freestream flow before the shock, the area 
ratio is based on the upstream Mach number, and AS is the average entropy rise 
of the fluid passing through the shock between Y^ and the wedge centerline. 

The latter quantity can be determined by (1) taking a number of equally spaced 
points between Y^ and the centerline, (2) using eq. (C-2) to calculate the 

associated shock angle (and consequently the entropy rise at the point), and 
(3) averaging the entropy rise over all such points. 

Combining eqs. (C-4) , (C-6), and (C- 7) 


Y 

s 


cos 6 

m 



AS 

R 


cos 6 

a 


(C-8) 


Eq. (C-6) may be used to calculate r. For a specified leading edge radius, 
the value of x (and consequently the physical dimensions of all other norma- 
1 i zed quant i ties) is 


x D = z (C-9) 

r 

Eq. (C-3) yields 


= X 

+ r sin 

6 


m 


m 


= X 

+ t 5 i n 

6 + r s in 6 

(C-10) 
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The shock standoff distance (from the center of the leading edge radius) 


i s 




D 


(C-11) 


THREE-DIMENSIONAL CASE 

In determining the location of the shock in the quasi-three-dimensional 
section, the shock is first treated as a simple two-dimensional shock in the 
manner described above. The flow angle into the shock is then assumed to vary 
linearly across the blade passage from the value at the wedge to the value at 
the shock intersection point on the suction surface. The inclination of the 
shock in the field is then assumed to be equal to the local flow angle minus 
the shock angle calculated for the simple two-dimensional shock a similar dis- 
tance from the wedge centerline. In addition, the standoff distance of the 
shock is assumed to be the same as in the two-dimensional case. With the stand- 
off point and shock inclination angle determined in the above manner, the shock 
location is determined by 


m 


sh 


m 

so 



0 

so 



(C-12) 


Since the suction surface angle at the shock intersection point is not 
known a priori, the above technique implies an iterative solution. Rather than 

iterate for the suction surface, an approximate technique was used in the 
section design program. 


LEADING EDGE SHOCK SHAPE--STREAMLI NE NO. 4 


The Mach number into the streamline No. 
magnitude of the wedge angle are, 


4 


leading edge wedge and the 


M 1 = 1.502 


= 4.34 deg 
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The correspondi ng downstream Mach number and shock angle far from the 
leading edge radius are 

M oo - 1.351 
a oo - 46.91 deg 

The maximum deflection angle, sonic deflection angle, and sonic shock 
angle correspondi ng to the inlet Mach number are 

6 = 12.17 deg 

m 

6 s = 11.75 deg 
= 62.24 deg 

The corresponding values of and from eqs. (C-1) and (C-2) are 

X = 1 .210 
s 

Y s = 0.728 

The average entropy rise to the sonic point, normalized to R is 

AS _ 0.0524 
R 

The cor respondi ng value of & and r are 

1 = 0.651 
r = 0.0927 

The actual leading edge radius and the value for X q are 

r = 0.0073 in. 

X = 0.079 

o 

The value for X and the shock standoff distance (from the center of the 
c v 

leading edge radius) are 
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X c = 1 .360 

X = 0.028 
so 


LEADING EDGE SHOCK BLOCKAGE 


In Appendix D, a correction for leading edge shock loss is incorporated 
into the overall leading edge blocking calculation for development of suction 
surface. The technique used to calculate the magnitude of this correction is 
disoussed in the following paragraphs. 


Consider a streamtube of width dy^ at a distance above the wedge cen- 
terline upstream of the shock. The shock angle, and hence the entropy rise of 
the fluid in the streamtube, can be calculated from eqs. (C-1) and (C-2). 


Downstream of the shock, the streamtube is located at a normal distance y 
from the wedge and has the width dy. Since the streamtube flow area is propor- 
tional to streamtube width, and the shock pressure recovery, in terms of 
entropy rise, is 


J5t2 

ptl 


AS 

R 


(C-13) 


the relationship between streamtube widths upstream and downstream of the shock 
is as fol lows : 



(C-1 4) 


where the subscripts 1 and 2 refer to conditions before and after the shock, 
(A/A*)^ and (A/A*) are the A*/A values corresponding to the upstream and down- 
stream Mach numbers, respectively, and the exponential term accounts for the 
change in A* across the shock. Eq. (C-l4) is written as 

dy = T dy 1 (C-1 5) 
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where F represents the term in brackets in eq. (C-l4). The normal distance 
of any arbitrary streamtube above the wedge surface can then be determined 

by integrating eq . (C-15) from the wedge centerline to Yj a> i.e., 



(C- 16) 


The blockage associated with the shock losses in the vicinity of the leading 

edge is taken to be the difference, normalized by the leading edge radius, 

between V for the real shock wave and Y (denoted Y .) for the idealized shock 
a a a 1 

wave attached to an infinitely sharp wedge. Denoting this quantity as V yields 


V 



or 




(c-17) 


Since the suction surface wedge angle is zero, F. = 1 • Extending the upper 
limit of integrations to infinity yields, finally 


00 

Y = 7/ (T - 1) dy t 

Jo 


(C— 1 8) 


LEADING EDGE SHOCK BLOCKAGE--STREAMLI NE NO. b 


Setting the wedge angle to zero and using the same inlet Mach number and 
leading edge radius, the calculated shock parameters for streamline No. b are 
as fol lows 

M 1 = 1.502 
6 = 0 

M 2 . = 1.502 
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= 41 .74 


CT 2i 

6 = 12.17 

m 

6 $ = 11.75 

a = 62.24 
s 

x s = 1.133 

Y s = 0.474 

1 = 0.425 
7 = 0.060 

x = 0.121 
o 

x = 0.150 
c 

/ 

Ax = 0.029 

It may be noted that despite the differences in the intermediate numbers, 
the shock standoff distance calculated for the upper and lower surfaces 
(6=0 deg and 6 = 4.34 deg, respecti vely) are approximately the same. This at 
least provides some justification for treating the actual wedge of 4.34 deg at 
an angle of attack of 2.17 deg as a superposition of an 8.68-deg wedge and a 
0-deg wedge, each at a 0-deg angle of attack. 

The normalized shock blockage for the above shock as calculated by a 50- 
step numerical integration to the upper limit 5 Y g is 

Y s = 0.332 

The corresponding incidence angle correction to streamline no. 4 is 

i f = 0.17 deg 
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APPENDIX D 


CASCADE START MARGIN 

TWO-DIMENSIONAL CASCADES 

Experimental studies of high-solidity, two-dimensional cascades with 
supersonic inlet flow indicate that at high back pressures, such cascades, 
in general, operate with a strong normal shock on the suction surface forming 
a bow wave ahead of the leading edge of the adjacent blade. As the back pres- 
sure is reduced with inlet Mach number held constant, the incidence angle 
decreases and the suction surface shock moves aft. Continued reduction gener- 
ally results in one of two possible occurrences: either the cascade chokes or 

the shock reaches the leading edge of the adjacent blade. In the latter case, 
further reduction of cascade back pressure causes the suction surface shock to 
"jump" past the effective throat of the cascade, and to be replaced by an 
oblique shock system in the leading edge region. When this phenomena occurs, 
the cascade is said to be operating in the "started mode." Prior to this 
occurrence, the cascade is operating in the "unstarted mode." Once the cascade 
starts, the incidence angle no longer changes with changes in cascade back 
pressure. The incidence angle at which start occurs (hereafter referred to as 
the "start incidence") is, in general, only a function of the inlet Mach 
number for a particular cascade. Further reductions in back pressure tend to 
draw the trailing edge normal shock out of the trailing edge region, and for 
sections which are similar in shape to the outboard rotor sections, this will 
result in fully supersonic operation of the cascade. 

START MARGIN DEFINITION 

In order to ensure the starting of a cascade operated in the manner 
described above, the cascade passages should be sufficiently open to allow 
the cascade to reach the start incidence without choking. The start margin 
is defined as the percentage by which the metal passage area exceeds the mini- 
mum allowable area. For calculating the start margin, the operating conditions 
closest to choke are assumed to occur with the cascade at start incidence with 
the unstarted inlet shock structure. This condition need not be the critical 
condition for cascades operating at high Mach numbers with low flow angles. The 
minimum allowable passage area in this case is simply the sonic area for condi- 
tions downstream of the inlet normal shock. In general, cascade experience 
indicates that start margins of 3 to 4 percent are adequate to insure starting. 
The distribution of that margin for rotor sections 1 to 8 i s presented in fig. 
D-1 . 


In order to calculate the start margin for a particular cascade, it is 
necessary to determine the starting incidence angle and the Mach number into 
the passage inlet normal shock. An approximate expression for the start 
incidence angle and shock Mach number for the simple case of a two-dimensional 
cascade with a straight suction surface is derived below. The application of 
this approximate expression to the design of the rotor outboard sections is 
discussed later in this appendix. 
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STARTING INCIDENCE 


Consider the flow field for a started cascade with an external wave pattern 
as iVlustrated in fig. D-2. Let the upstream Mach number and flow angle be 
denoted as and $ 9 respectively, and the corresponding quantities on the 

suction surface downstream of the wedges as M g and (3^ The start incidence 

relative to the suction surface i is by definition 

s 7 

i s = e«-e s (D-i) 

for started flow. The corresponding suction surface Mach number (and hence 
the Mach number into the passage inlet in the unstarted mode) can be calculated 
approximate ly if the suction side bow wave losses are neglected 

= v » + » s (0-2) 

where v is the Prandtl -Meyer expansion angle. 

The upstream flow area per blade per unit span is given by 

A oo “ Y COS Poo (D-3) 

where Y is the cascade pitch. The passage inlet flow area per unit span is 
given by 


Aj = Y cos P s - t 


(0-4) 


where t is the suction surface metal blockage. The change in flow area 
normalized to the inlet area is therefore 


A I - 


AA 

A ro 


Y cos (Pa, - i g ) - t 
Y cos Pa, 


(0-5) 


For small incidence angles such that cos i s « 1 and sin i « this area 
change is approximately 


AA , a v 

— = (t^ P J 


Y cos P, 


(0-6) 


Eq. ( D-6) then relates, for a given cascade, the required starting area 
margin to the starting incidence with respect to the blade suction surface. 
From continuity requirements, the change in flow capability (induced by the 
change in Mach number) between the upstream and suction surfaces must equal 
the change in available flow area. It can be shown that for small values of 
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Lead i ng edge 
ob 1 i que shock 



cond i t ions 



Figure D-2r-Leading Edge Region of a Single Two-Dimensional Cascade 
with Straight Suction Surface in the Started Mode. 


i and a weak oblique shock from the leading edge, the change in flow capacity 
is approximated by 




(D-7) 


Eq. (D-7) then defines the approximate, requl red area margin in terms of 
the upstream Mach number for the fixed cascade. Equating eqs. (D-6) and (D-7) 
yields 


t 

Y cos Po, 

S tan - \/m£ - t 

^co oo 


(D-8) 


THE OUTBOARD ROTOR SECTIONS 

At the design point, the design intent is to have Thp outboard sections 
operating in the started mode, the upper central section operating in the started 
mode with a strong, second shock internally, and the lower central sections 
operating in the unstarted mode. Due to the exp^ctptipn of a continuous leading 
edge shock surface, a gradual transition from started sections near the tip to 
unstarted sections near the hub is anticipated. There 1$ np operating point 
analogous to this transition region in the two-dimensional cascade. The 
extent of this transition region, the manner. in which it occurs, and the nature 
of the flow within it are questions that require further analytical and experi- 
mental investigations. 

In addition to the structure of the leading edge wave pattern at specified 
speed and back pressure, the behavior of this wave pattern at a specified 
speed and varying back pressure is open to investigation. In particular, to 
determine if the transition region moves smoothly up and dov/n the blade with 
varying back pressure, or "jumps" into and out of the rotor at certain critical 
values of back pressure with a consequent abrupt starting or unstarting of the 
tip section. The latter case would most probably be observable as a hysteresis 
phenomena in the fan operating map. Although this behavior has been looked 
for in many transonic and supersonic compressors, hysteresis of performance 
has not been observed. However, for entirely supersonic rotors there have 
been occasions where the sound level has made abrupt changes up and down. 

Independent of whether smooth or abrupt transitions occur, a question 
exists as to whether the rotor sections require more or less start margin 
than the corresponding two-dimensiona 1 cascades. This question is significant 
in that the start margin is a major parameter in the design of the outboard 
sections. In the present rotor design, no benefit was assumed to accrue 
relative to start margin due to the three-dimensional nature of the rotor 
flow field. 
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The techniques used to calculate the start incidence and the corresponding 
start margin for the outboard section are discussed below. 

SUCTION SURFACE DEVELOPMENT 


The suction surface of the outboard rotor -secti on as far aft as the first 
captured Mach line is developed by means of a direct analogy with the simplified 
two-dimensional cascade discussed previously. The corresponding simplified 
two-dimensional cascade quantities in the analogy are listed below. 


(1) Straight suction surface aft of the leading edge wedge 

(2) Flow conditions specified at upstream infinity 

(3) Leading edge wedge blockage 

(4) Uniform Mach number and flow conditions on the straight suction sur- 
face after correction for leading edge wedge blockage 


Analogous 
fo 1 lowi ng 

0 ) 

U) 

(3)' 


three-dimensional section items in the analogy are given in the 
list. 

Free streamline suction surface aft of the leading edge radius 
Calculated axi symmet ! ri c inlet conditions 
Leading edge radius and shock loss blockages 


(4) Varying Mach number and flow conditions along the calculated free 
streamline after correction for leading edge radius and shock loss 
b lockages 


In setting the finite leading edge radius equivalent to the leading edge 
wedge blockage, it is assumed that the dividing streamline is parallel to the 
free streaml ine- and, hypothetically, passes through the leading edge center. 
Consistent with this assumption, the suction surface thickness associated with 
the leading edge radius is set equal to the leading edge radius. The "leading 
edge shock blockage" refers to effective blockage associated with the detached 
position of the leading edge shock in the immediate vicinity of the leading 
edge radius. The technique used to calculate this blockage is presented in 
Appendix C. 


The free streamline suction surface incidence is calculated as the sum 
of two incidence angles i and i^.. The former accounts for the leading edge 

radius blockage and is calculated by replacing t by r. in eq. (D-8). The 

, L • t • 

latter accounts for, the shock. blockage and is calculated as 


(D-9) 


The technique used to calculate Y is discussed in Appendix C. 
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The total incidence correction is used to calculate the free streamline 
inlet Mach number. 

In summary, a free streamline suction surface is used in an attempt to 
reduce the wave disturbances upstream of the section, with the free streamline 
flow condition specified in a manner which is (at least to a first approxima- 
tion) consistent with started flow at the design point. 

APPLICATION TO STREAMLINE NO. 4 

For streamline no. 4, the inlet Mach number and flow angle are 

Mj = 1.479 

3 1 = 62.85 deg . 

and the leading edge nose radius and leading edge pitch are 

r = 0.0073 in. (0.01854 cm) 

Y = 2.13 in. (5*41 cm) 

The incidence correction i calculated from eq. (D-8) is 0.50; the shock 

s ^ 

blockage correction per eq. (D-9) for streamline no. 4 is 

i f - 0.17 

and hence, the total incidence is 

i + i ^ = 0.67 deg 

The values of Mach number and flow angle used in the free streamline 
calculations are 


M, = 1.502 

1 corr 


6, = 62.17 deg 

1 corr 3 


The shock blockage correction was introduced into the blade calculation 
procedure via the incidence factor term discussed in the text. The actual 
value used in the calculation was 0.18 deg rather than the more correct value 
of 0.17 deg. 
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APPENDIX E 


C 

D 


DCA 

f 



Ft 

y 

9 

h 

H 

i 

1 1' 
i 

s 

j 


NOMENCLATURE 


Chord, in. (cm) 
Diffusion parameter 


rotor 


V ' /V ' 
2 ' 1 


r 2 V 6 2 ~ r l V 6 1 

(r 1 + r 2 ) v ct 


stator = 1 - V^/V^ + 


r 2 v e 2 


V'ej 


< r 2 + r 3 ) V 2 a 


Double c i rcular arc 

Ratio of suction surface to total boundary 
Layer displacement thickness 
Suction surface transition factor 


Pressure surface transition factor 


Ultimate tensile stress 


Tensile yield stress 

2 2 

Gravitational acceleration, ft/s (m/s ) 

Radial height of stream tube, consistent units 
Boundary layer form factor, 6 /0 
Blade incidence angle, deg 

Section incidence correction due to inlet shock loss 

i 

Section incidence correction due to leading edge thickness 
Arbitrary index 
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t 

m 

M 

MCA 

N 

P 

P 

a 

Nb 

r 

R 

S 

T 

t 

Sn/C 

U 

V 
W 

Y 
Z 

a 

P 


Sonic line between leading edge radius and the 

Meridional distance in radial-axial plane, in. 

Mach number 

Multiple circular arc 

Rotational speed, rpm (rad/s) 

2 

Total pressure, psia (N/cm ) 

2 

Static pressure, psia (N/cm ) 


bow wave 
(cm) 


-AS 

Loss parameter e ° 


Pj/P 

<V T 2> 


2 

y/y - 1 


Blade number 
Radius, in. (cm) 

Gas constant, ft- 1 b^/1 b |T| -°R (m-N/kg-°K) 
Entropy 

Temperature, °R (°K) 

Blade thickness, in. (cm) 

Maximum thickness to chord ratio 


Wheel tip speed, ft/s (m/s) 

Absolute velocity, ft/s (m/s) 

Airflow, lb/s (kg/s) 

Blade pitch, in. (m) 

Axial distance, in. (m) 

Streamline in radial-axial plane, deg 
Flow angle measured from meridional, deg 
Stagger angle, deg 
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Inlet suction surface angle, deg 


0 


Iss 


Y 

6 



6 ° 

* 

6 

## 

6 

e 

n 


i* 

V 

S' 


Ratio of specific heats, C /C 
r ’ p v 

Ratio of total pressure to NASA standard sea level pressure 
of 14.696 psia (10.133 N/cm 2 ) 

Leading edge shock strength, deg 

Trailing edge shock strength, deg 

Pressure surface kink angle, deg 

Trailing edge wedge angle, deg 


Deviation angle with respect to blade metal mean camberline, deg 
Displacement thickness, deg 

Deviation angle with respect to envelope meanline, deg 
Function characteristic equation 
Adiabatic efficiency 


rotor = 


(p 2 /p 1 ) y " 1/y -i 

(T/rp - 1 


stage = 


(p 3 /p 1 ) y ' 1/y -1 
oy^) - i 


Parameter modifying the subsonic section surface curvature 
d istr i but ion 

Mach angle 

Prandtl-Meyer expansion angle, deg 

v + 8 

Characteristic parameter, — - — 

v - 8 

Characteristic parameter, — r — 
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p Fluid density, lb/ft^ (kg/m^) 

a Blade sol idity 

= conical chord/blade spacing (based on average radius) 
0 Streaml ine slope angle, deg 

to Angular velocity, rad/s 

to Total pressure loss coefficient 



and a^ = upstream total acoustic velocity 


stator 



P 


P 


1 

2 


to Profile loss coefficient 

P 

w Shock loss coefficient 

s 

Subscripts 

1 Rotor leading-edge station 

2 Rotor tra i 1 i ng-edge , stator leading-edge station 

3 Stator tra i 1 i ng-edge station 

00 Freestream just upstream of cascade 

or ss Suction surface 

p Pressure surface 

fs free stream 

T Tip 

MCL mean camber line 


M 


Meridional component 
LE Leading edge 

TE Trailing edge 

ex exit 

Q Tangential component 

z Axial component 

r Radial component 

Superscripts 
* Critical condition 

** Blade envelope parameter 


t 


Relative to rotating part 
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